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EXECUTIVE  SUMMARY 


A  three-dimensional  mathematical  model  of  a  seat,  (x;ciipant(s),  and  restraint  system  has  been 
developed  for  use  in  aircraft  crashworthiness  analysis.  Programs  SOM-LA  (Seat/Occupant  Mcxlel 
-  Light  Aircraft)  and  SOM-TA  (Seat/Occupant  Model  -  Transport  Aircraft)  combine  a  lurnptui 
parameter  model  of  aircraft  occupants  with  a  finite  element  model  of  tlie  seat  structure.  SOM-LA 
models  a  single  occupant,  whereas  SOM-TA  has  the  capability  to  model  up  to  three  passengers. 
The  intent  of  these  programs  is  to  aid  In  evaluation  of  the  performance  of  aircraft  seats  and  restraint 
systems  in  crash  environments.  Because  the  programs  have  been  wntten  for  use  primarily  by 
engineers  concerned  w'ith  the  design  and  analysis  of  seats  and  restraint  systems,  an  effort  has  been 
made  to  iiunimize  the  input  that  is  required  to  describe  the  occupants.  Characteristics  of  two 
standtird  occupants,  one  dummy  and  one.  human,  ai^e  included  within  the  program,  and  an  option  is 
provided  to  simulate  other  occupants  by  providing  additional  input  data.  The  structural  model 
includes  beam  elements  and  has  a  maximum  capacity  of  approximately  450  degrees  cf  freedom,  as 
determined  by  array  dinrsensions  within  the  programs.  'Fhe  beam  elements  can  accommodate  large 
plastic  deformations  and  include  the  capability  for  cross  section  reduction  due  to  local  instabilities. 
Four  different  cases  are  describcxl,  and  a  listing  of  input  data  is  provided  for  each.  These  examples 
are  (1)  a  simple  three-passeiiger  airline  seat  model  with  three  occupants,  (2)  a  single-ot:cupant 
general  aviation  seat  with  a  more  complex  structural  configuration  than  that  of  the  first  exarriple,  (3) 
an  energy-absorbing  helicopter  seat,  and  (4)  a  case  in  which  two  seat  rows  are  modeled,  in  order  to 
demonstrate  the  effects  of  passenger  impact  on  the  seat  backs  in  front  of  them. 

A  line-by-line  description  of  input  data  is  provided  in  an  appendix.  Another  appendix  includes 
examples  of  input  data  for  nonstandard  occupants,  seveTal  cushion  materials,  and  a  number  of 
structural  alloys.  Program  organization  is  described  in  detail,  as  are  the  funcbons  of  all 
subroutines.  A  complete  set  of  output  data  for  one  of  the  examples  is  also  included. 

Details  of  the  mathematical  models  and  solution  algorithms  for  the  SOM-LA  program  were 
reported  in  DOT/r^'AA/CT-8^/33-I  (March  ^983)  and  for  the  SOM-TA  program,  in  DOT/FAA/CT- 
86/25-1  ( ^  ugust  1986).  The  program  documentation  was  originally  presented  as  a  second  volume 
of  each  of  the  above  reports  (bOT/FA.A/Cr-82/33-n  and  DOr/FAA/CT-86/25-II).  Due  to 
program  nuxlifications,  those  separate  manuals  (the  second  volumes)  are  now  superseded  by  this 
single  document. 


1^0  imfisiL'Cim 


Proj^fanis  SOM-TA  (Seat/Occupant  Model  -  Transport  Aircraft)  and  SOM-LA  (Seat/Occupant 
-  Light  Aircraft)  combine  a  lumped  parameter  model  of  ;  craft  occupants  with  a  finite 
clement  mwlel  of  the  seat  siructure.  SOM-LA  models  a  single  occup.  nt,  whereas  SOM-TA  has  the 
capability  to  model  up  to  three  passengers.  'The  intent  of  these  programs  is  to  aid  in  evaluation  of 
the  performance  of  aircraft  seat  and  restraint  systems  in  crash  environments.  Because  the 
programs  have  been  written  for  use  primarily  by  engineers  concerned  witli  the  design  and  analysis 
of  seats  and  restraint  systems,  an  effort  has  been  made  to  minimize  the  input  that  is  required  to 
describe  the  occupants.  The  program  allows  simulation  of  one,  two,  or  three  passengers,  of  the 
same  or  different  sizes.  Characteristics  of  two  standard  occupants,  one  dummy  and  one  human, 
are  included  within  the  program,  and  an  option  is  provided  to  simulate  other  occupants  by 
providing  additional  input  data.  ITie  structural  model  includes  beam  elements  and  has  a  maximum 
capacity  of  approximately  450  degrees  of  freedom,  as  determined  by  array  dimensions  within  the 
program.  The  beam  elements  can  accommodate  large  plastic  deformations  and  include  the 
capability  for  cross  section  reduction  due  to  lo  a!  instabilities.  As  an  option  to  reduce  both 
modeling  complexity  and  cccution  costs  for  cases  where  only  the  restiaint  system  or  cabin 
( :  infiguration  is  of  ccincern,  or  for  cases  where  the  (I.;:,:ai;s  of  the  seat  de.sign  may  not  yet  be  known, 
a  rigid  scat  mode),  in  which  scat  pan  and  back  planes  defned  by  ifiput  are  maintained  in  fixed 
positions  in  the  aircraft,  is  avaiiable.  Details  of  the  inrithematic. !  models,  validation  programs,  and 
solution  procedures  are  contained  in  Reference  1  for  .SOM-i.A  and  Reference  2  for  SOM-TA. 
Input  instructions  for  the  vSOM  LA  and  SOM-TA  computer  programs  were  originally  preseiited  in 
References  3  and  4,  respectively,  which  this  re|X)rt  now  replaces. 

The  following  sections  of  this  report  present  insumetions  necessary  for  the  use  of  Programs  SOM- 
TA  and  SOM-LA,  and  infomtarion  to  enable  the  user  to  operate  fJie  program  most  efficiently. 

Sections  2.0  and  3.0  describe  program  input  and  output,  respectively,  including  options  available 
ro  the  user.  Section  4.0  outlines  an  efficient  prtK'edure  for  development  of  a  mathematical  mcxlel. 
Section  5.1  then  provides  detailed  descriptions  of  .sample  input  cases.  Appendix  A  defines  all 
input  variables,  line  by  line.  Appendix  B  provides  exampies  of  material  properties  and  occupant 
characteristics  required  as  input  data.  Appendix  C  describes  program  organization  and  the 
functions  of  aft  subroutines.  Appendi.:  D  di.splays  the  complete  set  of  output  data  for  the  e.xample 
descrilied  in  Sectio.i  5.1  and  Appendix  A. 


2,0  PROGRAM  INPUT  DATA 

Input  data  are  read  in  the  following  six  bloc  ks: 

1 .  Simulation  and  output  control  information 

2.  Cushion  properties 

3.  Restraint  system  description 

4.  Crash  conditions 

.  Occupant  description 

6.  Seat  design  information. 

A, II  input  data,  except  those  pertaining  to  the  seat  (Bltx'k  6),  are  read  by  subroutine  INPT;  the  seat 
data  are  read  by  subroutines  SRATIN,  BGF'OM.  CABIN,  and  RF.ADIN. 

The  cixtrdinate  system  that  is  fixed  to  the  aircraft  at  the  fl(K>r  has  the  following  positive  directions: 

X  Forward 
Y  Ixft 
/  -  Upward 

The  basic  input  data  deck  consists  of  a  minimum  of  26  lines  of  data  for  execution  of  Program 
SO.M  LA/S()M-'lA  with  one  passenger.  These  are  described  ifi  detail  in  Appendix  A.  f'he  basic 
case  makes  use  of  a  rigid  seat  model,  specified  by  NSFA  T  -  •  0  on  i,ine  3  Mtxleling  an  actual  seat 
with  the  fume  element  analysis  would  re>iuire  a  number  of  additional  lines,  beginning  with  Line 
27.  Rec|uesting  the  storage  of  plot  data  on  externa!  files,  unit  14  for  the  (H,cupani  and  unit  20  for 
the  sciit,  by  setting  NOPT  I'  >  0  on  Line  4  or  N,SPLT  >  0  on  Line  27,  requires  additional  lines  to 
ilescribe  the  plots.  .Mixleling  luinstandard  iKcupantts)  requims  an  .idditional  12  data  lines  for  each 
occupant,  after  Line  22.  If  the  setit  n.w  m  front  of  the  passenger(s)  is  to  bt*  mixleled,  six  atiditionai 
lines  of  data  must  tx‘  added  .it  the  end  of  the  input  tleck 

The  tollowing  sections  ol  this  chaptei  preseni  des^:np^lous  ot  each  ol  die  six  input  d:ita  blocks, 
mciikiing  niou*  detailec.!  delinitunis  of  the  atxive  'options,  !  me  bv  line  dc-scriplions  of  mimi  data 
!or  ,m  ex.imple  are  presented  in  .A[)|h  iuIi\  .A 

2  i  s IMLLLV ntjN, ,\NI).ULIl*L!T.idJNIKL)L iNTGKMATlGN 

2  1  i  .Systems  of  l.l'iits.  ITe  N1  'NTI  p.uana  tei  on  Line  |vniuis  die  user  to  s()ecitv  cither  die  SI 
Ol  Tnglish  system  ert  units  lor  both  in[iut  .iiui  output  il.iia  bnglisli  unus  are  ptx'senierl  ihioueliout 
l!ie  !n[ru!  instructions  in  this  re|>on  aiut  .ue  userl  in  the  s.imple  m|nil  cases,  in  the  SI  system  ot 
units,  .id  lei'gdis  are  s'lecitied  m  meteis,  m.isses  m  kilogiams,  and  lorces  oi  weigfils  in  newtons. 

2  I  2  Scdl .ULSicUti  I  be  NSI'.'X  I  n.u.mieier  on  Tiik'  '  ,dluw  ■,  du'  iisei  to  seleci  eii'u  i  .i  iieul  mmi 
iiusiel  Ol  .I  biiiie  eleiiumi  se.it  nxxlel  I  he  ripKl  se.ti  iiukIcI  (.  onsists  ol  i\',(i  pl.aies  ili.il  lepicseni  the 
MUiI  irtn  ,iiui  se.it  b.u  k  (he  posiiumsid  I'u  sc  jdaiies  ate  spevdi.e.i  tn  du  \  and  /  .ooidm.ile^  o! 
liicn  in'ei ''ta.  iH'ii  i.i  l.ilei.'l  line)  .uni  iwo  .ingles  w  ha  h  s|Ha  d\  llum  posmons  lel.mo-  lo  lun  i.'onuil 
.a  ii  i  er !  K  .i !  j  u.i  nr  s  i  e  sjx--  (i  \  r  l  \ .  !  he  length  oi  the  se.ii  p.in  ,iiid  the  heiu  li!  o!  the  sea,  I  I  si,  k  .n  a  u  -  ed 
io  date!  I'l  !i !  le  i (te  I n till s  ol  f ne  sin  i.iees  w  iliuii  hu  h  it.e  ^e.l^  p.iii  .m.l  hsu  k  ean  \  I oia  i' s  lo  i!ie 
esviip.inl  (lishioiis  ol  aqait  spciilitd  diiek  nesses,  ,ne  ou  hided,  on  lop  ot  ihe  se.n  .ind  se.ii 

i  sii  k  s Li :  I  a.v  (' 


The  type  of  seat  (siinglc-,  double-,  or  triple-occupant)  is  specified  on  Line  3,  along  with 
identification  of  the  positions  that  are  occupied. 

2.1.3  Occupant  Derrf.es  of  Freedom.  'Die  NDIM  parameter  on  Line  3  permits  selection  of  cither 
two-  (NDIM  =  2)  or  three-dimensional  (NDIM  -  3)  occupant  response,  'fhe  three-dimensional 
occupant  model  consists  of  12  rigid  segments,  illustrated  in  Figure  1,  with  rotational  springs  and 
clamj^rs  at  the  joints.  Each  of  the  torso  joints  possesses  three  rotational  degrees  of  freedom,  or,  in 
other  words ,  is  a  ball-and-socket  type  joint.  Because  of  the  hinge-type  motion  at  elbow  and  knee 
joints,  the  position  of  a  forearm  or  lower  leg  relative  to  an  uppter  arm  or  thigh,  respectively,  is 
described  by  one  additional  angular  coordinate.  In  total,  this  occupant  mcxiel  possesses  29  degrees 
of  freedom. 

An  alternative  occupant  model,  which  is  restricted  to  plane  motion,  is  specified  by  NDIM  -  2  on 
Line  3;  it  consists  of  1 1  segments,  as  shown  in  Figure  2.  Beam  elements  in  the  torso  and  neck  arc 
capable  of  flexural  and  axial  deformation.  Although  restricted  to  two-dimensional  response,  this 
occupant  option  does  permit  more  direct  evaluation  of  accident  severity  by  output  of  forces  and 
moments  in  the  spine  and  neck.  Restraint  system  forces  on  the  ellipsoid^  contact  surfaces  are 
computed  three-dimensionally,  but  only  the  X-  and  Z-components  are  used.  Therefore,  the  plane- 
motion  model  should  be  reserved  for  cases  in  which  both  the  impact  conditions  and  the  restraint 
system  are  symmetrical  with  respect  to  the  X-Z  plane. 

2.1.4  Output  Control  Data.  Ten  blocks  of  program  output  can  be  selected  on  Line  4.  Tlie  data 
include  time  histories  of  the  following  variables,  which  are  stored  during  solution  at  predetermined 
print  intervals: 

1 .  Occupant  segment  positions  (X,  Y,  Z,  pitch  and  roll)* 

2.  Occupant  segment  velocities  (X,  Y,  and  Z) 

3 .  Occupant  segment  accelerations  (x,  y,  z,  and  resultants)* 

4.  Restraint  system  loads 

5 .  Cushion  loads 

6 .  Aiiciaft  displacement,  velocity,  and  acceleration 

7 .  Injury  criteria,  including  spinal  forces  and  moments 

o .  Details  of  contact  between  (lie  occupants  and  the  scat  or  interior  .surfaces  in  front  of  them 

9.  Scat  strticture  nodal  displacements  and  forces 

1 0.  Scat  structure  element  stresses 


*lJpper  case  X,  Y,  Z  refer  tc  r'crtiaS  ur  at  re,  a  It  fixe  tl  ciKmlitiatc  s/sicm;  lower  m-A'  y,  /.  refei  to  'A-gtiiciii  lixe;i 
coondittates. 


®  SEFGMENT  MASS  CENTER 
O  JOINT  . ® 


Figure  1.  Twelve-segment  (genera!  three -dimensutnal)  tx;ciipant  mcxlel. 


Printer  plots  are  provided  for  occupant  segment  accelerations,  restraint  system  loads,  and  ciisiiion 
loads,  d’he  option  t)f  two  different  filters  is  also  provided  for  llie  (K-cnpant  segment  accelerations 
and  cushion  loads.  The  particular  occupant  for  vdiich  out[nil  data  arc  displayctl  is  spccd'ird  in 
Line  4. 

if  plots  are  reipiesied  for  the  iK-ciij-  'lu  and/or  seal  on  Line  4,  then  additional  lines  mus.  Ix"  incUaicd 
to  specify  plot  times  (up  to  ciglit)  and  viewing  angles.  As  explained  in  the  line  by  line  inpiu  data 
descrijnions,  if  plots  are  axjuestcd,  the  job  control  language  must  define  externa!  files  14  a'lii  20  to 
lie  saved  for  postprocessing. 

I’lograui  output  data  arc  described  tieifieT  in  LIupptc!  ..f 


2,1  5  Solution  Control  Data  The  occupant  model  utUizes  an  Adams -Moulton  prcdictoi  corrector 
solution  pr(x:cdurc  with  a  variable  step  size.  Data  on  Line  .S  control  the  .step  si/e  and  error  bounds 
for  the  solution,  a;,  well  as  the  duration  of  the  simulation, 

2.'^  SrxXjIs'DARY  ’MPACT/SHAT  BACK.  (DNTACr 

li  coiitaci  hct’vvcen  (l!C  fiasscnpcrs  and  the  scat  in  froni:  of  them  is  to  Ik."  snmilated  by  SOM  I'A,  ’lie 
i  iph!  [dane  sui  faces  illissuated  in  lapu'c  must  lx-  delt.'a'd  by  input  (d  liK:aiions  and  dmKawioi.’s  f  !' 
ilic  scat  back,  fiay  tal»k\  aitd  ,irm  rests  a!id  of  loice  denectioi!  liHH.Iioris  for  these  sui  laces.  1  >nnni,! 
excL  ution  ot  tfu'  jar.tyuani,  tiie  d-i.^ance  between  eac  h  «,;!  these  suiKice>  aiid  ilte  2ti  Indy  K.citac! 
ciiilisoids  K  cal(  aLued.  A  distance  less  tiian  zero  indicates  jx-netranon  of  a  suiiacc;  a  contact  force 
is  coni|ai!eti  fioin  tins  [jcneu  inon  and  appiici.!  at  the  appropriate  ixunis  on  the  ho(„lv  aini  (he  seal 


ActiKil  iransport  aircraft  scats  have  hticks'  that  arc  hitigcti  to  rotate  ff)rward  if  pusticd  I'n'rn  !)l1  did 
'i’lic  SOM  -'l’A  program  pcmiiis  tiie  scat  back  to  rotate  forv/ard  ;ilx)ut  a  transverse  hinge  axis  at  the 
base  of  the  back.  Should  an  (Kveupant  strike  any  of  the  seat  back  surfaces.,  the  moir.en.  of  the 
impact  force  v,'ith  respect  t(t  the  hinge  axis  is  coinpiitcd  as  the  product  of  the  nonrial  component  of 
file  force  multiplied  by  the  distaticc  from  me  axis  to  the  contact  point.  I'he  appliei.l  moment  is 
coniitarekl  widi  the  resisting  moment,  and,  if  the  net  ntoraent  is  greater  than  zero,  seat  back  angular 

acceleration  a  is  calcuhUtjd  according  to 

EM  -  .M,..,  1y-  (1) 

where,  --  the  moment  applied  to  the  seat  back  by  the  cxfcuptuit, 

=  the  resisting  moment  of  the  seat  back,  and 
iy  -  the  moment  of  inertia  of  the  seat  back  about  a  transverse  (Y-) 

axis  at  the  hinge  point. 

The  resisting  niornent  depends  on  the  cutrent  angular  displacement  of  the  seat  back  from  its 
initial  position  and  is  calculated  by  interpolation  in  a  table  of  input  moments  and  displacements 
which  pnxiuce  the  function  shown  in  P'igure  4.  Loading  along  the  initial  slope  up  to  point  2  is 
elastic,  so  that  if  a  slight  "bump"  were  to  occur,  causing  an  angular  displacement  less  than 
DDBO(2),  the  seat  back  would  return  to  its  initial  position  after  removal  of  the  load.  Once 
DDBO(2)  has  been  exceeded,  the  resisting  moment  remains  cor.stard  at  the  "breakover"  moment 
FFBO(2)  up  to  a  displacement  DDBO(?).  Beyond  DDBO(?),  the  resisting  moment  increases 
rapidly,  simulating  a  mechanical  stop. 

Use  of  the  seat  back  contact  feature?  is  specified  by  input  of  IOUT{4)  >  0  on  Line  4.  Si.x 
■idditioria!,  input  lines  ate  then  required  follow'ing  Line  26.  The  input  format  is  described  in 
Appendix  A, 3. 

2.3  CUSHION  PROPERTSFS 

Seat  cushion  forces  applied  to  the  oecupiutt  iiKxlel  are  calculated  from  cushion  dellections  Pased  ()ri 
an  exponential  mlationsh’p: 

F  =  C(e”'^-!)  (2) 

l..incs  8  through  10  require  input  cf  the  C  and  B  c(x:fficiefits  for  this  equation,  along  with  dam[)ing 
coefficients  and  thicknesses.  The  force  denecl ion  lehuiotiship  for  the  seat  cushio't  itlso  includes 
c;)mp!iai'ce  of  the  occujiant  butUX'ks.  Thew-'dre,  the  mlatiouship  ior  ati  (K,'cupant  s.ttitig  directly  on 
a  liard  seat  pan  would  be  tlie  fcuco  tiefleciion  curve  lor  the  (K'cupaul  hiitux.'ks.  Several  sarn|de 
force  delleclion  curves  with  their  appropriate  coefficients  are  I'lnivKied  in  Afiixmdix  B. 


2.4  RFSTRAINT  SY.SI'Ff 


Several  iest.raiiu  system  eonfigurations  are  avadable  in  .SOYI  'I'A/SOM  1,,'\;  lap  k'it  only,  lap  with 
(liagona!  sb.oulder  beit  over  eillier  shoulder,  and  tlouble  shoulder  belt  vvitli  or  'A  ithmit  a,  la|i  bell 
tiedown  siraj).  .As  specified  on  l  .ine  3,  Irotli  the  lap  belt  itiui  shouhlei  lairne.  s  can  he  alta.clicd  ti' 
c ntn  r  ih  '  airlr.-nne  or  tl'C  seal. 


riic  ion  e  def!cc:;(.’ii  i  iiaractenst.c,  ol  the  restiaint  svsieiu  wi'lihine  are  [irovnliai  iw  nijaii  of  lafiie:.. 
oi  iot..'es  iiiid  sli  lln,^.  l'i()peili..s  C't  repre'seiuaiive  restraint  welihang  i>'pes  are  nu  lueied  in 
■Appendix  F 


o 


Moment,  M 


Figure  3.  Seat  back  contact  surfaces. 


Fij'ua’  4.  .Scat  back  resisting  fiionicnt  liutction. 


For  a  seat  with  a  shoulder  harness  in  which  an  inertia  reel  is  mounted  on  the  seat  b<tck  and  a  length 
of  inertia  reel  strap  is  passed  along  the  seat  back  to  a  slot  above  the  txxupant  shoulders,  the 
XTRAL  parameter  on  Line  14  defines  the  length  of  the  shoulder  strap  behind  the  seat  back. 

2.5  IMPACT  CONDITIONS 

Six  components  of  the  acceleration  of  the  aircraft  ctxtrdinate  system  are  provided  on  NIMFF  Lines 
beginning  with  21  A:  time,  X,  Y,  Z,  j'aw,  pitch,  and  roil.  Acceleration  components  are  directed  in 
the  aircraft-fixed  coordinate  system. 

2.6  OCCUPANT  DFSCRIPTION 

The  IMAN  parameter  on  Line  3  identifies  the  type  of  occupant  being  simulated.  Initial  positions 
for  each  passenger  are  specified  by  data  on  Line  22.  Data  for  the  standard  occupants,  a  50th- 
percentile  U.S.  male  and  a  50th-percentile  (Part  572)  anthropomorphic  dummy,  are  included 
within  the  program.  For  nonstandard  occupants,  additional  data  may  be  provided  (for  each 
occupant)  following  Line  22  to  define  segment  lengths,  center  of  mass  locations,  weights, 
moments  of  inertia,  contact  surface  radii,  properties  for  the  spine  and  neck,  and  compliances  for 
the  chest  and  abdomen  under  restraint  system  loading.  Examples  of  these  data  are  included  in 
Appendix  B. 

2.7  SEAT  DESIGN  INFORMATION 

2.7.1  Rigid  Seat  Option.  For  cases  where  the  details  of  seat  response  are  not  important  or  not 
worth  the  greater  execution  costs  that  would  be  incurred  by  the  use  of  the  finite  element  structural 
model,  a  rigid  seat  option  is  provided.  Plane  surfaces  representing  the  seat  pan  and  seat  back 
support  the  cushions  and  remain  fixed  in  the  aircraft  coordinate  system,  e.xcept  where  the  energy¬ 
absorbing  option  is  used. 

2.7.2  Simplified  Energy-Absorbing  Seat  Option.  If  the  SEA'I  M  parameter  on  Line  24  is  greater 
than  0,  a  simplified,  two-degree-of  freedom  seat  model  is  used.  Intended  for  use  in  simulation  of 
a  guided  energy-absorbing  seat,  this  model  permits  the  stroking  of  a  rigid  seat  bucket  in  a 
prescribed  direction.  Because  elastic  bending  of  the  supporting  frame  has  been  observed  in  testing 
of  such  seats  and  may  influence  occupant  response,  the  second  degree  of  freedom  is  added  to 
simulate  rotational  elasticity  of  the  frame. 

Although  the  finite  element  analysis  can  provide  a  complete  evaluation  of  a  seat's  crashworthy 
performance,  the  simple  stroking  seat  nuxlel  can  prove  useful  in  other  aspects  of  seat  design.  For 
example,  the  two-degree-of-freedom  iiKxlel  can  aid  in  economically  estimating  the  optimum  energy 
absorber  limit  load  for  protection  of  occupants  of  various  size,  as  well  as  in  evaluating  alternative 
restraint  system  configurations. 

Input  data  for  this  seat  nxxiei  include  the  weight  of  the  movable  part  of  the  seat,  the  direction  along 
which  it  will  stroke,  the  mass  moment  of  inertia  with  respect  to  a  iateral  axis,  force  del  lection 
characteristics,  and  unloading  sIoik's 

2  '  Fimle  Flcmcni  Structural  Analysis.  ’I’he  I’lniie  element  seal  miHlel  coniameil  in  Ih'ogram 

.")(  >M  LA/SOM  i'A  uses  lieain  elements,  flie  beam  elemeiils  can  acconiiiUHlaic  large,  plastic 
clctormaiions  ami  liK'ali/cd  buckling  of  elcniciits  with  hollow  cross  sections,  i'hc  program  has  a 
ca[iacity  lor  /5  luxlcs  atul  degrees  ol  trccdom.  However,  a  more  severe  restnetion  is  phiced 
on  ll;e  si/e,  N.  ol  the  master  stilTness  matrix,  given  by: 

N  Mi:g  I  MUD  M2  *  Mi  O  M(4  >  I )/.’  (M 


S 


where  MUD  is  the  length  of  the  maximum  upper  diagonal  of  the  banded  stiffness  matrix  given  by: 

MUD-6*(J+ 1)- 1  (4) 

MEQ  equals  the  total  number  of  degrees  of  freedom  and  J  equals  the  maximum  difference  between 
node  numbers  across  elements  in  the  model,  as  illustrated  in  Section  4.0. 

As  determined  by  array  dimensions  in  Program  SOM-LA/SOM-TA,  the  quantity  N  is  limited  to 
16,000. 

The  finite  element  seat  model  uses  an  unconditionally  stable  solution  algorithm.  However,  stability 
does  not  necessarily  imply  convergence  to  the  correct  solution,  and  solution  accuracy  will  depend 
on  the  size  of  the  time  step,  a  smdler  time  step  yielding  more  accurate  results.  Because  the  seat 
step  size  is  governed  by  that  for  the  occupant  model,  reducing  DMAX  and  DMIN  on  Line  5  will 
produce  a  more  accurate  solution.  However,  little  improvement  can  be  expected  in  reducing  the 
seat  step  size  below  that  normally  required  for  stability  in  the  occupant  solution. 

Material  properties,  irxluding  a  three-slope  approximation  to  the  stress-strain  curve,  are  provided 
on  Lines  33-35,  which  must  be  repeated  for  each  material  used.  (The  number  of  mateiials  is 
specified  as  NUMAT  on  Line  27.)  To  assist  in  input  of  material  properties,  summaries  of  input 
data  for  metals  typically  used  in  seat  frames  are  presented  in  Appendix  B. 

Beam  cross-sections  can  be  either  open  or  closed,  but  a  plastic  problem  requires  a  closed  cross- 
section  to  generate  all  the  terms  required  by  the  tangent  stiffness  matrix.  If  plastic  deformation  of 
an  open  "I"  is  anticipated,  the  cross-section  can  be  modeled  as  a  closed  "box"  beam,  which  is 
equivalent  for  one  bending  direction,  provided  that  the  erroneous  properties  for  other  bending 
directions  can  be  tolerated. 

The  Nl/MDS  parameter  on  Line  27  specifies  the  number  of  nodes  that  are  attached  to  the  aircraft 
structure.  Then,  floor  attachment  conditions  are  specified  on  Line  43,  one  of  which  must  be 
inserted  for  each  of  the  NUMDS  nodes.  Element  cross-sections  are  described  by  data  on  Lines  36 
and  37,  which  must  be  repeated  for  each  cross-section,  the  number  of  which  is  specified  by 
NSECT  on  Line  27. 

Nodal  coordinates  are  provided  on  Line  38,  which  is  repeated  for  each  node  in  the  model 
(NUM.NP  on  Line  27)  and  for  each  beam  pointer  node  (NCC)RD  on  Line  27).  As  illustrated  in 
Appendix  A,  the  pointer  node  is  required  to  specify  the  initial  orientation  of  the  y-axis  of  a  beam 
cross-section.  A  real  node  can  be  used  as  a  pointer  ritxle.  or  (NCORD)  additional  nodes  can  lie 
added  solely  to  serve  as  pointers. 

Element  data  are  provided  on  Line  39,  which  must  be  repeated  for  each  element  (NIJMEL  on 
Line  27).  Data  for  each  element  inch  de  identification  of  its  end  ntxles,  the  pointer  ntxle  that  is 
used  to  orient  the  cross  section,  the  cross  section,  the  material,  and  end  release  conditions. 


3.0  £RQGEAM.QjUmJXQAIA 

Output  data  are  available  from  the  following  four  sources; 


1 .  Printer  (unit  6) 

2.  Occupant  position  plots  (unit  14) 


3 .  Seat  structure  plots  (unit  20) 

4.  Plots  of  other  data  (unit  26) 


which  are  described  further  in  the  following  sections. 

3.1  PRINTED  QU'fPUT 

Printed  data  can  be  selected  from  the  ten  blocks  listed  in  Section  2,1.4.  The  interval  at  which  these 
data  are  printed  is  selected  in  subroutine  INPT,  based  on  the  total  solution  time.  The  interval  is 
sized  to  provide  a  maximum  of  5 1  lines  (appmximately  one  page)  for  each  vaiiable.  For  examnle, 
a  solution  time  between  0.100  and  0.150  sec  results  in  a  print  interval  of  0.003  sec,  a  solution  time 
between  0.250  and  0.300  sec,  an  interval  of  0.006  sec,  etc. 

Accelerations,  severity  indices,  vertebral  forces  and  moments,  and  restraint  system  forces  are 
printed  in  tabular  and  graphical  formats.  Other  data  are  provided  in  tabular  form  only. 
Acceleration  output  data  are  compute<l  each  0.001  sec,  equivalent  to  a  1  KHz  sampling  rate.  Input 
Line  4  provides  the  option  of  applying  a  Class  180  (3(30  Hz)  or  Class  60  (100  H/)  liter  to  the  data 
prior  to  their  printing. 


If  specified  in  input  Line  4,  data  for  up  to  eight  plots  of  occupant  posi'ion  can  be  stored  on  external 
file  14.  The  times  for  these  plots  are  defined  on  input  Line  7,  along  with  viewing  angles,  which 
itfe  illustrated  in  Figure  5.  The  right-side  view  of  Figure  6  was  obtained  using  an  angle  of  zero 
degrees.  The  front  view  of  Figure  7  was  obtained  using  an  angle  of  90  degrees. 

The  IOUT(4)  parameter  on  Line  4  can  be  used  to  draw  the  image  of  the  seat  in  front  of  tliat  being 
mi>deled.  A  value  of  0  causes  no  seat  to  be  drawn.  A  value  of  i  or  2,  respectively,  prtxiuces  a 
plot  of  the  forward  seat  in  its  undefoimed  or  deformed  position  and  uses  subroutine  IMPAC'F  for 
prediction  of  contact  with  the  forward  seat.  Ilie  seat  image  is  spaced  according  to  the  SPII  O! 
piu'ameter  (Lane  49). 

The  job  control  language  used  in  executing  SOM-LA/SOM-TA  must  define  external  file  14  as  a 
permanent  file  to  be  saved.  The  (Kcupant  plotting  program  can  then  be  executed  using  this  same 
permanent  file  as  input. 

3.3  SEAT  STRUCTURE PLQ'fS 


Just  as  described  in  Section  3.2  for  (xrcupant  position,  data  for  plots  of  the  seat  structuri'  can  be 
requested  on  I.ine  27,  As  shown  in  Figure  8,  n(Hie;s  are  indicated  and  numbered.  'I’hc  viewer 
position  for  the  seat  structure  is  defined  by  both  elevation  and  a/.imudi  angle  ,  '•  a.nd 
rcsi-'cctively,  as  shown  in  1-igure  9.  The  view  of  I'igure  H  w'as  obtained  with  0  20  degrees  and  o 

--  4,3  degrees 


Figuix'  5.  Definition  of  plot  viewing  angle. 


The  job  control  language  must  save  external  file  20  for  subsequent  use  as  input  to  the  seat  plotting 
program. 

3.4  ADDITIONAL  DATA  FOR  PLOITING 

Although  the  printer  plots  of  accelerations  and  forces  me  probably  satisfactory  ouq^ut  for  most 
purposes,  therc  may  be  cases  where  plots  with  a  higher  level  of  rc.>olutjon  are  desired.  Also,  pen- 
drawn  plots  may  be  required  for  use  in  rept)rts.  To  meet  these  netids,  32  variables  are  written  on 
externa!  file  26  at  0.00! -sec  intervals.  The  data  .are  written  in  eiuier  FI0.3  or  F10..'i  format  and  tu-e 
arranged  as  illustrated  in  Figure  10. 


PROGRAM  SOM-TA 


PROGRAM  SOM-TA 
TRANSPORT  ASRCRAF 
TIME  =  0.0000  SEC. 


SEAT 


S' 


0  =  Azimuth  angle  in  X-Y  plane  in  degrees  (-1 80'’  <  e  <  +180”) 
(j)  -  Elevation  angle  in  degrees  (-90”  <  (j)  <  +90”) 


Field 

Forma) 

Viirial^k 

1 

FI  0.5 

Time  (sec) 

2 

F10.3 

Aircraft  X-accci  (G) 

3 

FI  0.5 

Aircraft  Z-accel  (G) 

4 

FI  0.5 

Aircraft  res.  accel  (G) 

5 

FI  0.5 

Aircraft  res.  vei.  (ft/sec) 

6 

F10.5 

Aircraft  res.  displ.  (in.) 

7 

F10.5 

Seat  X-accel  (G) 

8 

FI  0.5 

Seat  Z-accel  (G) 

9 

FI  0.5 

Head  x-accel  (G) 

10 

FI  0.5 

1  lead  z-accel  (G ) 

11 

FI  0.5 

Head  res.  accel  (G) 

12 

FI  0.5 

Chest  x-accel  (G) 

13 

FiO.5 

Chest  z-  accel  (G ) 

14 

FI  0.5 

Chest  res.  accel  (G) 

15 

FI  0.5 

DRl 

16 

F10.3 

Seat  cushion  force  (lb) 

17 

F10  5 

Pelvis  x-accel  (G) 

18 

F10.5 

Pelvis  z-accel  (G) 

19 

FI  0.5 

Pelvis  res.  accel  (G) 

20 

FI  0.3 

Lumbiir  tixial  load  (lb) 

21 

F10  3 

[..umbar  y-moment  (in.-Ib) 

22 

I- 10. 3 

Neck  axial  load  (lb) 

23 

FI  0.3 

Neck  y  moment  (in. -lb) 

24 

i!0.3 

Pack  cushion  force  (lb) 

25 

1  10.3 

Right  lap  bell  force  (lb) 

26 

1  10.3 

I  .eft  lap  bel(  force  (lb) 

27 

1  10.3 

Right  siioulder  ixdt  force  (Ih) 

28 

FlO  3 

1  eft  shoulder  Ix'll  force  (lb)'* 

29 

110  3 

l  .nergy  ab.sorber  force  (!b) 

30 

[•10.5 

Seat  displacement  (in. ) 

31 

F  lO  3 

l  iHilresi  .X  force  (Ihi 

12 

110  1 

i  (H)lR‘sl  /.  foree  (lb) 

‘Rcplai'c*!  by  scii!  annul.if  ilisjilai  ciiu'ni  tiM*  I'lKTpv  il'surhiiij;  M'as  iiunk'l  I'a  ah  NSi' A  I  0  a.iul  M  \  IM  •  0) 


f'lj’iire  1!)  t )al.i  tonnal  for  cxHTiial  li!'-  -f: 


If. 


4.0 


This  chapter  is  intended  to  guide  program  users  through  an  efficient  process  of  preparing  input 
data.  The  recommended  procedure  is  summarized  in  Table  1.  It  is  suggested  that,  if  time  permits, 
one  or  more  of  the  sample;  cases  described  in  detail  in  Chapter  5  be  run  initially  in  order  to  be 
certain  that  the  program  runs  properly  on  a  particular  computer  system.  Storage  of  plot  data  on 
permanent  files  and  subsequent  access  of  these  files  using  the  related  occupant  and  seat  plot 
programs  should  be  attempted  first  with  the  samiple  cases  to  assure  that  the  plotting  programs  are 
compatible  with  the  computer  system  and  that  the  job  control  language  is  structured  properly. 


_  TABLE  1 .  SUMMARY  OF  INPUT  DATA _ _ _ 

1 .  On  sketch  of  seat,  locate  aircraft  floor  and  coordinate  system. 

2 .  IXMcate  restraint  system  anchor  points. 

3 .  Ixxate  footrest  and/or  heel  position  (at  Z  =  0). 

4.  Estimate  initial  position  angles  foi  occupant  upper  torso,  head,  and  arms. 

5 .  ftepare  input  data  for  and  run  rigid  seat  case  fo*"  short  time. 

6 .  Plot  occupant  initial  position  and  check  whether  it  appears  reasonable. 

7 .  Add  seat  structure  input  after  Line  26. 

8 .  Run  short  case  with  civmplete  input  data. 

9 .  Check  plot  of  seat  structure  at  initial  lime. 

1 0.  Run  complete  ca.se. 


The  essential  starting  point  for  any  simulation  case  is  a  sketch  of  the  seat  of  interest,  o.n  which  the 
aircraft  floor  and  aircraft  coordinate  axes  can  be  Ux’ated.  On  this  sketch,  the  restraint  system 
anchor  points,  which  can  l)e  fixed  to  either  the  seat  or  the  aircraft  structure,  can  be  located,  as  can 
the  position  of  a  fexttrest  or  }iedal,  if  applicable.  The  required  seat  design  data  for  a  rigid  scat  case, 

i. e.,  the  l(x:ations  and  dimensions  of  the  seat  pan  and  back,  can  then  be  detenriined.  Both  ihe  scat 
cushion  and  back  cushion  are  assumed  to  be  plane  surfaces  parallel  to  the  scat  pan  and  back 
surfaces.  Using  an  average  cushion  thickness  in  the  area  of  confaei  between  trie  occupttnt  and  the 
seat,  the  cushion  surfaces  can  now  be  added  to  the  sketch. 

Initial  angular  positions  of  the  torso,  head,  and  ami  scgment.s'  tire  required,  dhe  torso  segnienl.s 
can  be  tissumed  parallel  to,  or  one  or  two  rlegrees  forward  of,  the  seat  back,  dhe  position  of  flic 
head,  in  a  normal  seated  position,  wouhi  range  from  vertical  to  several  degrees  froward  of  vertictil, 
as  illustrated  in  the  Section  5.0  sample  cases. 

!n  order  to  be  ceiltiin  that  the  occupant  mitia!  posiuon  is  rea.sonabic  for  the  configmatioi!  being 
studied,  it  is  wise  to  run  a  short  rigid  seat  case  prior  to  adding  the  seat  structure  iiq-iut.  Sttiriing 
vv.tli  ti  s.mali  viuue  of  fintil  time,  Tig  on  [,ine  .'s,  such  as  0.010  see,  tiie  posilioti  ttnd 

accelerations  of  the  <K:eup;int  segments  and  the  exterrtal  forces  ettn  be  chcckra!  prior  to  initiatiiv:  a 
longer,  pei'ha[)s  more  expensive  ease.  The  plot  data  saved  on  unit  11  by  SOM  'FA  can  then  F.e 

ii. put  to  the  (K.a  u[nint  [doi  progr.tm  and  the  initial  posthon  ol  all  the  (Heu[)ant  segments  reviewed. 


If  the  occupant's  initial  position,  as  calculated  by  subroutine  INITIL,  is  geometrically  impossible, 
the  program  will  be  stopped  and  informative  messages  printed.  An  example  of  this  type  of  error, 
commonly  encountered  in  initially  running  a  case,  is  in  attempting  to  locate  the  heel  position 
beyond  the  reach  of  the  legs.  If  the  input  parameters  yield  a  geometrically  feasible  initial  position 
and  NOPLT  >  0  on  Line  4  and  TOPLT  =  0  on  Line  7,  then  data  for  a  plot  of  the  initial  position  will 
be  stored. 

Once  the  desired  initial  position  has  been  achieved  for  the  occupant,  the  input  data  for  the  finite 
element  seat  structure  model  caw  be  added.  The  NSEAT  parameter  on  Line  3  should  then  be 
changed  from  0  to  1  to  signify  modeling  a  nonrigid  seat.  Once  again,  prior  to  running  a  complete 
simulation,  a  case  with  a  small  TF  should  probably  be  run  in  order  to  check  the  seat  structure  plot 
at  the  initial  time. 

When  it  has  been  determined  that  both  the  occupant  irdtial  position  and  the  seat  structure 
configuration  are  as  desired,  a  complete  simulation  case  can  be  run. 


5.0  MMaJE^I^LAIIQIiJCASES 


This  section  contains  descriptions  of  input  data  for  sample  simulation  cases.  A  simple  airline  seat 
model  with  three  occupants  is  used  to  illustrate  program  capabilities,  including  the  preparation  of 
input  data  and  output  data  available  from  the  program. 

A  complete  line-by-line  input  data  listing  is  discus.sed  in  Section  5.1.  Tlic  .contents, of  Appendix  A 
show  explicit  input  requirements  for  this  case.  Particular  features  of  different  cases  are  discussed 
in  Sex;tions  5,2  -  5.4. 


5.1 


LmSiLLl: 


Descriptive  t'tles. 


>  LISTING  IN 


LiaSLS.:  NDIM  =  2  for  plane-motion  simulation;  IMAN  =  1  for  a  standard  50th- 

percentile  (Part  572)  dummy;  NSEAT  =  1  for  use  of  finite  element  seat  model; 
IRSYS  0  for  lap-beit-only  restraint  system  configuration;  lEUKL  =  0  (used 
to  identify  the  nature  of  shoulder  belt  attachment  to  the  lap  belt  if  IRSYS  >  0); 
HdBLT  =  1  for  lap  belt  attachment  to  seat;  IS)TNS  -  0  (iis^  to  specify  whether 
shoulder  harness  is  attached  to  aircraft  or  seat  if  IRSYS  >  0);  NIMFF  =  4  for 
input  of  tour  points  in  rime  and  acceleration  (Lines  21A-21D)  specifying  a 
trapezoidal  pulse  shape;  NUNTT  -=  I  for  Engli.sh  units;  NDCC  =  3  fur  three 
axupants  to  be  simulated;  OYTE  =  3  for  a  thre;sr-occupajU  seat  type;  ISEAT(l) 
-  1,  ISE,AT(2)  “  i  and  !SEAT(3)  =  1  to  specify  that  all  seat  positions  are 
occupied. 


Line.4:  K)UT(l)  =  I  and  iOUT(2)  ~  1  request  segment  position  and  velocity  data; 

IOUT(3)  ~  2  requests  occupant  segment  acceleration,  filtered  with  class  180 
digital  filter,  IOUT(4)  0  calls  for  no  secondary  impact  sinmlation;  IOUT(5) 
=  1  requests  restraint  system  forces;  I01JT(6)  --  1  requests  spinal  loads  and 
injury  criteria;  IOUT(2)  =  I  requests  seat  external  forces  fillered  with  a  class  60 
digital  filter;  iOUT(8)  -  1,  IOUT(9)  -  1,  and  lOUT(lC)  -  1  request  sear 
structure  deflection,  support  reactions,  and  beam  loads  and  stresses, 
respectively;  NOPIT  --  8  for  eight  occupant-position  plots;  IPASS  =  2  to 
sjjecify  that  die  output  data  M'ill  be  stored  and  printed  for  die  center  occupant. 

Line  5:  '11  =  ()  and  TF  =  0. 17 5  sec;  DMAX  and  ! )MIN  set  to  ().(XX)3  sec  for  fixed- time 

step  integration;  integration  error  boutids,  EUR,  and  ELR,  set  ai  10  and  0,1 
percent,  respectively;  initial  time  step  size,  DTI,  is  same  as  fixed  step  size  (If 
DMAX  is  not  ti:ie  same  as  litMlN,  D'H  should  ixr  set  equal  to  DMIN). 


CKP'I'IN  0.025  for  interval  (in  seconds)  at  which  restart  data  are  to  be 
written  on  unit  25. 


UofiSiZAcZIi: 

Liod.ii: 


CX:cupant  plot  data  are  to  be  writien  on  unit  14  at  the  rimes  specified  on  these 
eight  iine.s,  (Eight  lines  are  retjuired  by  input  of  NOP?  T  =  8  on  lone  4  )  ‘"’lot 
viewing  angles  (Figiiie  5)  are  all  0  deg,  indicating  right-side  views. 

Conil-ined  seal,  bottom  custiion  ami  oceupaii!  luitlocks  forcc-deflectioi! 
rcl.Mtionshtj)  of  the  form  F  --  760  (eZOft  j)  sliowu  in  f  igure  1  !:  d.m.pifig 
luseiflcieru  of  2.4i)  tt  zero  load.  A  cushion  tlnckncss  of  f  5  .u.  w.u.  useti. 


Figure  11.  Exponential  approximation  of  force-deflectioii  characierisucs  for 
dummy  and  1.5  in.  thick  cushion. 


Idr.;  . :  Seat  back  cushion  properties  wid:  same  load-deflection  cuive  as  bottom 

cushion. 

Line  10:  Headrest  cushion  properties  with  same  load-deflection  curve  as  seat  back  and 

scat  bottom  cushions. 


Line  11:  Forces  and  strains  lor  2.0-ji>.  nylon  webbing. 

Line  12A:  L.ap  belt  anchor  point  ctxjrdinates  for  pa.;sengei  No.  1 . 


yn£J21i: 

LmLj,2£^ 

Liri£..0: 


I  .ap  beb  anchor  point  coordinates  for  passenger  h'o.  2. 

Lap  belt  anchor  point  ewtruinates  for  passenger  No.  3. 

Sbo'ilder  belt  force-dcneciion  properties  (blank  not  used  in  i!iis  case  because 
LR.SVS  =0  on  Line  3). 


iJnes  14A  -  i4C: 


Shoulder  belt  anchor  pe.int,  Rlj'Kl..,  and  XTRaI.  fo"  occispa'ns 
rcspecd''ely  (blank  •  fiOt  iisc(i  in  tics  t;ase). 


anej 


20 


titlkNl 


Line  15:  Lap  belt  tiedown  strap  properties  (blank  -  not  used  in  this  case). 

Lines  16A  16C:  Tiedown  strap  anchor  point  c<x>rdinates  for  occupants  1,  2,  and  3  (blank  -  not 
used  in  this  case). 


line  17:  Restraint  system  damping  coefficients  and  slack  all  zero. 

Line  18:  COEFFS  =0.18  and  COEF'FR  =  0.25,  friction  coefficients  for  seat  cushion 

and  floor,  respectively;  XFR  and  ANGFR  =  0.0  (used  only  if  a  footiest  is  to 
be  modeled). 

Line  19:  Initial  position  of  "aircraft"  coordinate  system  in  inertial  system  is  assumed  to 

be  (0.0,  0.0,  0.0)  with  yaw,  pitch,  and  roll  also  zero. 


Line  20:  Initial  velocity  of  "aircraft"  coordinate  system,  30  ft/sec  in  X-direction. 

Lines  21 A-21D:  Table  of  times  and  floor  accelerations  (X-component  only,  in  this  case).  Four 

lines  are  included  for  die  trapezoidal  pulse  shown  in  Figure  12,  as  requiied  by 
setting  NIMPT  =  4  on  Line  3. 


l-igurc  1 2.  .Anproxunation  to  sied  accderaiion. 


1  » 


Line.  22-1: 


Line  22-2: 
Um22zl: 
Line  23: 


Ussi24z26: 
Line  27: 


Initial  position  for  passenger  No.  1.  GAM(1,1)  and  GAM(2,1)  -  -16  degrees; 
GAM(3,1)  =  7  degrees;  GAM(4,1)  -  -16  degrees  for  upper  arms;  GAM(5,1) 
-  60  degrees  at  elbow;  heels  at  32  in.  (these  coordinates  are  illustrated  in 
Figure  A-5);  YPASS(l)  =  -20  in.  for  Y-coordinate  of  mid-plane  for  nassenger 
No.  1. 

Initial  position  for  passenger  2.  Same  as  line  22-1  except  YPASS(2)  =  0  in. 

Initial  position  for  passenger  3.  Same  as  line  22-1  except  YPASS(3)  =  20  in. 

The  coordinate  system  was  located  under  the  rear  edge  of  the  seat  pan,  so  that 
XSEAT  =  10  in.;  the  height  of  the  seat  pan  is  12  in.;  seat  pan  and  back  angles 
are  8  and  16  degrees,  respectively.  Seat  pan  length  and  width  are  15.15  in. 
and  20  in.  respectively.  The  height  of  the  seat  back  is  39  in.  (These 
coordinates  arc  illustrated  in  Figure  A-8). 

Blank  due  to  use  of  finite  element  seat  model  (NSEAT  =  1  on  Line  3). 

The  finite  element  seat  model  is  illustrated  in  Figure  13.  NUMNP  =  20  (real 
nodes);  NUMEL  -  27;  NUMAT  =  2  for  two  materials;  NUMDS  =  4 
(restrained  nodes  on  the  floor);  NCORD  =  2  (beam  pointer  nodes,  numbered 
21  and  22);  and  NSECT  =  2  for  two  beam  element  cross  sections;  NSPLT  =  8 
for  eight  seat  plots. 


1,2, . 22  Node  numbers 

(D.C?}.  .  @  Beam  element  numbers 


I6gurc  !  3.  Sample  case  three  passenger  seat  finite  element  rntxlel. 


7  7 


LimL2S: 


KNTRL(  l)  =  5  indicates  that  up  to  5  iterations  will  be  used  at  each  time  step. 
KNTRL(2)  =  5  indicates  that  5  load  increments  will  be  used  to  enforce  the 
floor  warping. 


Lines  29A-29H:  Seat  plot  data  are  to  be  written  on  unit  20  at  the  times  sp'^cified  on  these  eight 
lines.  (Eight  lines  are  required  by  input  of  NSPLT  =  U  on  Line  27.)  All  eight 
plots  are  to  be  made  with  aii  elevation  angle  of  +20  deg  and  an  azirnutli  angle  of 
+  45  deg,  i.e.,  viewed  from  left-front  quaiter. 

Line  30:  Nodal  output  data  aie  requested  for  nodes  1  through  20,  indicated  in  Figure  13. 

Line  31:  Beam  loads  and  stress  output  data  are  requested  for  elements  1  through  27, 

indicated  on  Figure  13. 

Line  32:  Seat  structure  output  at  intervals  of  0,025  sec. 


Lines  33-35:  There  are  two  groups  of  material  properties  corresponding  to  NLIMAT  =  2. 

Material  type  No.  1  is  2024-T4  aluminum;  material  No.  2  is  4130  steel. 

Lines  36-37:  There  are  two  groups  of  cross-section  properties  (NSECT  =  2),  shown  in 

Figure  14.  The  circular  tubing  cross-section,  defined  first,  is  approximated  by 
eight  segments;  the  square  cross-section  is  made  up  of  four  segments.  The 
orientation  of  the  cross  section  is  specified  by  the  beam  pointer  node,  which 
locates  the  beam  y-axis  (The  beam  coordinate  system  is  illustrated  in  Figure  A- 
13). 


-0.593,-0,593  -0  593,0.593 


Note:  All  Dimensions  in  inches. 

f  igure  14.  lilemen!  cross-sc  eiion  mixlels  used  for  sea!  struelure  l>eam  eknneiils 


LineJS 


LineJg: 


Line  40: 


Line  41: 


Line  42: 


I.ines43-44: 


Twenty-two  lines  of  luxlal  coordinate  data  corresponding  to  20  real  nodes 
(NNODE  =  20)  and  2  pointer  ntKles  (NUMNP  -  2).  Newie  point  ItK'ations  and 
number  are  shown  in  Figure  13. 

Twenty-seven  lines  of  lK;am  clement  dat.-j  corresponding  to  NIJMHL  “  27. 
Beam  eletix'-nt  connectivity  and  numbers  are  shown  in  Figure  13. 

Seat  pan  cushion  load  is  distributed  on  nodes  (5,  6,  13,  14)  for  the  first  (right) 
occupant,  (6,  7,  14,  15)  for  the  second  (middle)  c»ccupant,  and  (7,  8,  15,  16) 
for  the  third  (left)  occupant.  The  order  in  which  the  seat  pan  nodes  are 
specified  is  shown  in  Figure  A- 14. 

Back  cushion  load  is  distributed  on  nodes  (5,  6,  17,  18)  for  the  first  (right) 
occupant,  (6,  7,  18,  19)  for  the  second  (middle)  occupant,  and  (7,  8,  19,  20) 
for  the  third  (left)  occupant.  The  order  in  which  the  seat  back  nodes  are 
specified  is  shown  in  Figure  A- 14. 

Lap  belt  loads  are  applied  at  nodes  (9,  10)  for  the  first  (right)  occupant,  (10, 
11)  for  the  second  (middle)  occupant,  and  (11,  12)  for  the  third  (left)  occupant. 
Three  lines  are  used  for  the  three  passengers. 

Five  lines,  which  specify  the  constraint  conditions  at  the  bottom  of  the  seat-leg 
elements.  For  nodes  1,  2,  3,  and  4  all  forces  and  only  moments  in  Z-direction 
can  be  supported.  In  addition,  node  2  is  moved  down  in  Z-direction  by  0.5 
in,,  corresponding  to  user-specified  (pitch)  floor  warping  condition. 


A  complete  lis’  ,ng  of  the  input  data  for  this  san:iplc  case  is  presented  in  Figure  15.  Examples  of 
plots  genex-atec  .  this  case  are  presented  as  Figures  16-18. 
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PROGFiAM  SOM-TA 
TRANSPORT  AIRCRAFT  SEAT 
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PROGRAM  SOM-TA 
TRANSPOF^T  AIRCRAFT  SEAT 


TIME  =  0.1750  SEC. 


5.2  SAMPIJi  CASH  NO.  2:  PRQDUC?  ION  GENERAL  AVIATION  SHAT 

Fnont  and  side  views  of  the  prociuction  general  aviation  seat  treated  lucre  as  an  example  are  shown 
in  I'igure  1 9.  Note  that  the  ctxtrdinate  system  has  been  placed  on  the  floor  at  the  centerline  of  the 
seat  far  enough  aft  that  all  points  on  the  seat  will  be  positive.  Although  this  is  not  a  requirement, 
such  location  of  tlie  coordinate  system  does  facilitate  preparation  of  input  data. 

Because  of  the  nonsymmetric  restraint  system,  a  three-dimensional  occupant  simulation  is 
requested  by  NDIM  ==  3  on  Line  3.  As  illustrated  in  Figure  20,  the  seat  structure  was  modeled 
using  28  nodes  and  36  beam  elements.  The  seat  siruciure  is  fabricated  of  6061-T6  aluminum 
alloy.  The  cross  section  of  all  beam  elenncnts  is  illustrated  isi  Figure  21  along  with  the  rectangular 
approximation  utilized  in  the  model.  A  listing  of  input  data  is  presented  as  Figure  22.  Because 
neither  the  lap  belt  nor  the  shoulder  harness  is  attached  to  the  seat,  the  restraint  system  nodes  are 
not  required.  The  seat  has  one  fore-and-aft  adjustment  locking  pin  at  the  left-front  track 
connection.  Tire  track  connections  are  assui.ied  to  consnain  nodes  1 , 2,  and  15  against  translation 
in  the  Y  and  Z  direction  but  leave  them  free  in  the  other  directions.  Node  16,  on  the  other  hand, 
where  the  adjustment  locking  pin  is  located,  is  constrained  in  all  directions  except  Z  rotations. 
Some  judgment  is  required  as  to  the  ability  of  the  adjustment  lockir.g  pin  to  resist  rotations  about 
the  X  or  Y  axes. 
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Figure  22  (contd).  Input  daui  listing,  case  no.  2. 


A  production  energy-absorbing  helicopter  seat  was  tested  at  CAML  The  test  configuration  is 
illustrated  in  Figure  23.  A  complete  listing  of  input  data  is  presented  as  Figure  24. 

Headrest  properties  are  provided  on  Line  10.  In  addition  to  lap  belt  and  shoulder  belt  properties 
and  locations  on  Lines  11-14,  the  lap  belt  tiedown  strap  of  the  five-point  restraint  system  is 
described  on  Lines  15  and  16.  The  webbing  is  a  low-deflection  polyester  type,  whose  load- 
elongation  propertie;;  are  illustrated  in  Figure  B-5.  The  five-point  restraint  system  is  indicated  by 
IRSYS  =  4  on  Line  3. 

As  shown  in  Figure  23,  the  seat  was  rotated  on  the  horizontal  sled  in  order  to  simulate  a  near¬ 
vertical  impact.  The  input  acceleration  is  input  in  both  X-  and  Z-components,  on  Lines  21A-21H. 
The  pitch  of  -73  degrees  is  entered  on  Line  19. 

Energy  absorber  data  are  entered  on  Line  25.  The  load-stroke  characteristics  for  the  seat  are 
illustrated  in  Figure  25.  The  guide  tubes  shown  in  Figure  23  are  oriented  4  degrees  from  the  Z- 
axis,  and  this  angle  is  input  on  Line  24,  along  with  the  movable  seat  weight  of  60.6  lb.  It  is  this 
nonzero  sear  weight  that  causes  the  stroking  seat  model  to  be  used 

Line  24  includes  the  unloading  slope  of  4308  Ib/iri.  and  the  damping  coefficient  of  0.55  Ib-sec/in., 
which  was  determined  b}^  matching  the  measured  energy  absorber  force-time  history,  A  moment 
of  inertia  of  148  !b-in.-sec™  with  respect  to  the  aircraft  coordinate  system  was  estimated  for  the 
seat.  Rotational  stiffness  parameters  on  Line  26  were  estimated  from  static  tests  of  the  seat. 
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Configuration  for  dynamic  test  of  energy-absorbing  helicopter  seat 
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Figure  24  (contd).  Licting  of  input  data,  case  no.  3 
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Figure  25.  Energy  absorber  load-stroke  characteristics 


5  4  SAMPIJ'  CASE  NO.  4:  SEAT  BA CK  CON  I  ACI~ 

'I'his  example  simulates  one  of  a  serie.s  of  tests  conducted  at  the  FAA  (hs'il  Acroinedical  Institute 
(CAMI),  in  which  two  nnvs  of  seats  were  uistalled  on  the  deceleration  sled.  This  lest  condition 
vva.s  intended  primarily  to  observe  the  lunuence  v)f  aft  siratcd  passengers  on  seat  loading  (Ref.  5). 
nte  first  test  conditions  simulated  are  iho.se  of  test  AS7  (MO,  for  which  the  seats  were  installeil  at  a 
row  pitch  of  .10  in.  (SPl’l'ClI  ■=  30.0  on  l.ine  49),  and  the  .scat  back  movenicni  .in  the  forward  scat 
was  unrestricted.  The  "breakover"  moment,  which  would  resist  the  scat  back's  rotation,  was  set  in 
the  simulation  at  (lO  lb- ft,  the  approximate  lower  end  of  tlie  range  measured  by  ("AMI  tliiriiig  llieii 
test  program.  The  sled  deceleration  that  was  measurcil  in  the  test  was  digiti/eci  md  aiiplievi  as  input 
io  the  nuxlel;  the  digitized  pulse  is  shown  l\i  figum  2(i.  A  senes  of  (K'cigiant  pUits  sho's'ing  the 
scat  tiuck  rotation  is  incluitcrl  as  lagnic  37,  and  a  coiujiletc  listing  in  figiiic  3S. 
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Figure  26.  Test  sled  acceleration,  test  A87040. 
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5'iUiiTe  2H  I  isiiiig  ol  irpiif  (Lila,  v  as'‘  ncj.  I 


i 


0 . 1800 

0,00 

-13,0 

-13.0 

10.0 

-18.0 

AB.O 

-1,3. 0 

-13.0 

10.0 

-18.0 

BD 

--T  :■) .  0 

-13.0 

10.0 

-18. 0 

4  8.0 

2.2 

11.1 

5.00 

13.0 

20.0 

14 . 0 

16.0 

10.0 

22.0 

3.0 

7  60.0 

0.68 

3.0 

760.0 

1  .  0 

2.4 

7  60.0 

1.0 

2.4 

10.0 

8C.0 

30.0 

720.0 

720.0 

3500  0 

0.0349 

1.274 

2».0 
?H.  0 
2  8,0 
18.0 


0.0 


1  .  S2  3 


-18 . 0 
0.0 
18 . 0 
4  6.0 


16.0 


00021  1  .W 
G002:'000 
00022010 
00022020 
00023000 
00024  OOoj 
000250001 
00026000! 
0004500o| 
ooo46ond 


00047000| 


0004»OOC 


00050^  00 


Figure  28  (contd).  Listing  of  input  data,  case  no.  4. 
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APPENDIX  A 


INPUT  DATA  REQUIREMENTS 

in  tliis  Appendix,  a  line-  by-line  description  of  the  input  data  I’equired  by  Prograrxi  SOM-i  xA/SOi^'l- 
TA  for  example  no.  1,  described  in  section  5. ) ,  is  presented.  As  desenbed  in  sccuon  2,  there  are 
also  a  nivraber  (,»f  optional  lines  of  data,  such  as  Lines  22A  through  22L,  which  are  u.sed  only  for  a 
nonstandard  occupant.  These  are  included  follo  wing  the  Line  22  input.  beginnitAg  on  ^  A-32. 

Lines  24  tlirough  26  for  tlie  energy-absorbinp  seat  oteion  are  not  used  in  example  no.  1  and  are, 
therefore,  blank.  An  exio  ptle  of  th-riv  :  se  c.tn  be  found  .0,  exainple  no.  .1,  de.^'''ribed  in  sectiou  5.3. 
Data  for  the  lin'te  element  .f’ct  model  l>egins  with  Line  27,  which  follows  page  .A  -57. 

Input  data  for  rnodeh-rg  Cv)  i-act  b-:vwfcen  i.c-'upaptr  a.'n*  the  seat  backs  in.  front  of  L  .;;rn,  used  in 
example  tfc.  4,  begin  ■  n  page  A-Si,  foUcwiing,  the  input  h  -exanpie  no.  V 


i 


LINES  1  AND  2:  Case  Ideritirication 
DES£SiEfiQI;i:  I’iile  of  case  on  two  lines. 


. 1 

2| 

!T234567§W 

T2l45C7S50j 

$^^5678^1 234567595  IE4567|5^5^() 


3 


I  NANIEI  1 

NAME2 

THREE-PASSENGER  TRANSPOR'F  AIRCRAFT  SEAT 

SAMPLE  CASli  NO.  1  ! 

mw 

R^RMAT 

mmiENTS 

NAMEl 

7A10 

Alphanumeric  title  of  case  to  be  centeied  at  top  of  printed  output  and 
plots. 

NAME2 

7AI0 

Se^'-ond  line  of  alphanumeric  data. 

LINE  3:  Case  Control  Parameters 

DES CRIPTIQN :  Type  of  case,  type  of  seat,  and  occupant  locations. 

EQRMAIi^mjaMlPLE: 


r _ _  .1 

HHHK 

IIIIIIMIIIIIIllllllflllll^^ 

HHIIillHIIC 

tWtg&ffiriagil 

mmm 

pDlM 

MAN 

NSEAT 

IRSYS 

IBUKl 

ILBLT 

ISWJS 

NMPT 

NUNir 

NOCC 

iriTPE 

_ 

i(>s*va 

mg 

iiig 

1 

. . i 

i 

1 

hkhh 

mm- 

— 

_ 1 

1 

1 

mil 

FORMAT 

CQNTlWrS 

NDIM 

15 

Defiiiition  of  occupant  degrees  of  freedom 

NDIM  =  2  :  Two-dinM;nsional  (plane  motion)  simulation 
NDIM  =  3  :  Three-dimensional  sirr.ulation  (default). 

MAN 

15 

Identification  of  occupant 

IMAN  =  0 :  Standard  50tli-pei-cenale  male  human 

IMAN  =  1 :  Standard  50th-percentile  (Part  572)  dummy 
IMi\N  =  2  :  Nonstandard  human 

IMAN  =  3  :  Nonstandard  dummy. 

NSEAT 

15 

Seat  model 

NSEAT  =  0 :  Rigid  seat  model 

NSEAT  "  1 :  Finite  element  seat  model. 

IRSYS 

15 

Restraint  system  configuration 

IRSYS  =  0:  Up  belt  only 

IRSYS  =  I  :  Diagonal  she  Ider  belt  over  right  shoulder 
IRSYS  -  2  :  Diagonal  shoulder  belt  over  left  shoulder 
IRSYS  =  3  :  Di^uble  shoulder  belt. 

IBUKL 

15 

Buckle  connection  type  (see  Figure  A-1). 

LL3LT 

15 

Up  belt  attachment 

Il.BL'r  =  0  :  Attached  to  aiifiame 

IIBLT  -  1  :  Attached  to  seat. 

!SHNS 

15 

Shoulder  himiess  attachment 

ISHNS  =  0  ;  Attached  to  aiiframe 

IS  HNS  =  1  :  Attached  to  seat. 

NIMPT 

15 

Number  of  noints  in  table  of  aircraft  acceleration  vs.  time 
(uete'Tiiines  number  of  L  ine  21  inputs  recjuired,  a  maximum  of  40). 

NUNiT 

System  of  units 

NUNIT^-^O;  Si  units 

NUNl'I  -  1  :  English  uni's. 

N(x:c 

1.5 

Niinitcr  of  CKauipara!'  to  tx',  imdi'lcd. 

A 


OTPE 


ISEAT 


15  Seat  type  in  terms  of  occupant  positions 

ITYPE  =  I  :  Single  seat 
ITYPE  =  2  :  Two-passenger  seat 
ITYPE  =  3  :  Three-passenger  seat. 

315  Loc:ations  of  occupants,  specifying  whether  a  given  seat  position  is 

occupied  (1)  or  empty  (0) 

ISEAT(1):  Right-most  position 

ISEAT(2) :  Center  position  for  three-passenger  seat,  left 

position  for  two-passenger  seat 
ISEAT(3) :  Left  position  for  thrce-passenger  seat. 


A  I 


1  =  Shoulder  belt  fixed  to  buckle 

2  =  Shoulder  belt  and  one  side  of  lap 

belt  are  one  length  of  webbing 

3  =  Shoulder  belt  and  lap  belt  attached 

to  fixed  point 

NOTE:  BUKL  parameter  is  defined  on  lines  14A,  14B,  and  14C. 


Figure  A  !.  Types  of  buckle  connections  specified  by  IBUKL  on  line  3. 


LINE  4:  Output  Selection 

DESCRlFnON:  Definition  of  output  data  to  be  stored  for  printing  and  number  of  plots. 

EmMALMJD^AMELE: 


1 

4 

5 

7 

vmmm 

1 

I 

1 

mmm 

1 

i 

lOUTf]) 

!OUT(31 

IOUT(4) 

IOLnr(5|OUT(6) 

[OUT(7) 

[OUT(8) 

!OliT(9) 

“ 

Kx;r 

(70)- 

NtmT 

ITRMX 

IPASS 

Hi 

1 

_ 1 

2 

(] 

_ L.  1 

2 

1 

1 

Ml 

_ i 

■^'1 

_ 1 

HELD 

lOUT 


RIRMAT 

1015 


CONreNTS 

Vector  of  O's,  I's,  2's 
printed  (1 , 2,  or  3)  or 
lOUT(l) 
IOUT(2) 
IOirr(3) 
IOUT(4) 
iomx5) 
10171X6) 
10177(7) 
10171(8) 
lOUT(9) 
lOlLr(lO) 


and  3's  indicating  which  output  data  are  to  be 
not  printed  (0) 

Occupant  segment  position 
Occupant  segment  velocity 
Occupant  segment  acc6leration(l) 

Secondary  impact  prediction(2) 

Restraint  system  forces 
Injury  criteria 

Seat  external  loads  (cushions,  floor)  (1) 

Seat  structure  deflections(3) 

Seat  structure  support  reactions 
Stresses  in  seat  structure  beam  elements(4). 


NOPL;r 

15 

Number  of  requested  occupant  position  plots  (up  to  20). 
(Determines  number  of  Line  7  inputs  to  be  included.) 

mmx 

15 

Number  of  iterations  in  initially  seating  occupant(s).  (Default  ~  5.) 

IPASS 

15 

Identification  of  occupant  for  which  output  data  (position,  veIcKity, 
acceleration,  belt  loads,  etc.)  are  stored  and  printed 

IPASS  =  1  :  Right-most  passenger 

IPASS  =  2  ;  Center  position  for  three-passenger  seat  or 
left  position  for  two-passenger  seat 
IPASS  =  3  ;  Left  position  for  three-passenger  seat. 

The  example  specifics  a  three-passenger  seat  which  is  fully  occupied.  Output  data  are  stored  and 
printed  for  the  center  passenger. 


(1)  For  IOU'r(3)  and  IOUT(7),  an  input  value  of  I  results  in  unfiltercd  outpiil.  A  value  oi  2  or  t  results  in 
application  of  a  class  180  (3(K)  Hz)  or  class  (:>()  (UK)  H.',)  filter,  rcs{>cilively. 

( 2 )  10UT(4)  --  0  :  No  seciintiary  iiiifKict  prediction  and  the  forward  seal  is  not  plolted. 

IOU'r{4)  ■:  !  :  .Suhroutine  IMPAf'T  is  calleii  htr  jirecliction  of  contact  with  the  scat  bach,  and  (;cni|)ajil 

plots  sliow  die  forward  .seal  in  its  nndeforinezi  position, 

IO'.)T(4)  “  2  :  .Subroutine  IMl’ACT  is  called  for  [trcdiction  ot  contact  wilh  inc  .seat  tc.ick  ,  and  oi  cii();!ni 

plots  sltow  the  forward  .scat  deforniation  as  dial  vtf  the  scat  iH  ing  itUKlckat, 

(  a  )  II  10U'i'(8)  -  i ,  data  on  I  .inc  30  wili  Ire  u.scd  ti  select  the  tUHies  for  stress  OLit()nt  dcf'lcciion 

(  4  )  If  !()U'I'(  HI)  1  _  on  i  UK  3!  will  Ih'  used  to  .sclcci  the  fvjani  elcntcnis  lor  siress  output , 


A  0. 


LINE  5:  Simulation  Control  Data 

DESCRIITION:  Parameters  for  control  of  solution  duration,  step  size,  and  enor  bounds. 

FORMAT  AND  EXAMPIJi: 


r  I 

--  5 

3 

,  - 1 

5L„  6 

7 

1 

1 

1 

mmm 

1 

1 

'  TI  I 

TF 

DMAX 

DMIN 

.  EI.R  _  ^ 

Dn  1 

o.iSd 

■H&n 

HHuISiS 

mMSSSm 

FIELD 

FORMAT 

CONTENTS 

Tl 

FIO.O 

Initial  solution  time  in  seconds.  Nonnally  taken  as  0. 

TF 

FIO.O 

Final  solution  time  in  seconds. 

DMAX 

FIO.O 

Maximum  step  size.  A  value  of  0.001  sec  has  been  used 
successfully. 

DMIN 

FIO.O 

Minimum  step  size.  A  value  as  large  as  0.001  sec  has  beer,  used 
successfully,  but  the  use  of  very  stiff  restraint  system  webbing  or 
seat  cushions  may  require  a  smaller  value,  such  as  0.00001.  Tne 
solution  can  be  accomplished  with  a  fixed  step  size  by  setting  DMIN 
=  DMAX. 

EUR 

FIO.O 

Maximum  bound  on  error  between  predictor  and  corrector.  A  value 
of  0.05  to  0.10  is  suggested,  corresponding  to  a  range  of  5  to  10 
percent  If  the  error  in  any  variable  is  larger  tlian  this  value,  the  step 
size  is  halved,  maintaining  solution  act:uracy. 

ELR 

FIO.O 

I-X)wer  bound  on  error  between  pretlic  or  and  corrector.  A  value  of 
0.001  is  suggested,  corresponding  to  0.1  percent.  If  the  error  in  all 
variables  is  smaller  than  this  value,  the  step  size  is  doubled, 
preventing  the  computer  execution  cost  from  becoming  excessively 
high. 

Note:  Becau.se  doubling  the  step  .size  multiplies  the  tmneation  error 

in  the  Adams-Moulton  integrator  by  a  factor  of  2^,  ELR  sliould  be 
chosen  less  then  ElJR/32  if  the  advantages  of  doubling  arc  not  to  be 
short-lived. 

DD 

FIO.O 

Initial  step  size,  norniaily  set  equal  to  DMIN. 

A  7 


LINE  6:  Reiitait  Data  Inter\'a! 


DESCRIPTION :  Tiire  interval  at  which  data  are  to  be  written  on  unit  25  for  potential  use  in 

subsequently  restarting  solution. 


FORMAT  A  ND  EXAMPLE: 


]. 

2 

u, 

4 

5\  6 

7 

1234567890 

1234567890 

mmm 

1234567890! 1234567890 

1234567890 

CKPTIN 

0.025 

FIELD  FORMAT  CONTENTS 

CKPTIN  FI 0.0  Time  interval  in  seconds. 


A 


LINE  7:  Occupant  Plot  Times  and  Viewing  Angles  (number  of  lines  required  =  NOPLT  on  Line  4) 


Times  when  occupant  plot  data  are  to  be  stored  on  unit  14,  which  must  be 
saved  as  a  permanent  file  for  subsequent  plotting.  Viewing  angles 
corresponding  to  times  are  measured  in  degrees  in  the  horizontal  plane,  as 
illustrated  in  Figure  A-2.  An  angle  of  0  degrees  results  in  a  right-side  view; 
90  degrees,  a  front  view;  and  180  degrees,  a  left-side  view. 


FORMAT  AND  EXAMPLE: 


"—-IT—  2 

3 

4 

5 

6"  7 

te4567890 

1234567890 

1234567890 

"1234567890  1234567890 

TOPLT  !  ANGVU 

0.1751  O.Q 

EIEm  FORMAT  CONTENTS 


TOPLT  FI 0.0  Plot  f  me  (sec). 

ANGVU  FI 0.0  Occupant  viewing  angles  (deg). 


Viewing  angle  =  a 

Viewing 

F-’^ositlon 


i'igiin.:  A  2.  Dcimiiion  nt  iu'i'upaiii  pic?  virwui!.’,  anglf 


LINE  8:  Combined  Sea:  .shion  and  Occupant  Buttocks  Properties 

DESCRIPTION:  Force- deflection  characteristics  and  damping  for  seat  cushion  and  buttocks 

combined;  thickness  for  seat  bottom  cushion,  llie  force,  F,  is  computed 
from  total  deflection,  6,  according  to  F  =  C(e®5  -  1). 

FQRJ^IAT  AND  EXAMPT£: 


1  1 

5 

3 

5 

6 

7 

!TIJ456780O 

“12^4567890 

1234567890 

1234557590 

1234567890 

1^5575901 

1234567890 

\  CSC 

b5c“ 

DPS^ 

THSCE 

_  . .  1 

1  76^ 

0.50 

1.50 

„  .  _  1 

FIELD 

FORMAL’ 

CONTENTS 

FIO.O 

Ct  -'fl  lent  C  i 

1  ove  equation  (lb). 

BSC 

i-.O.O 

CoefTtcie  i:-.  B  i. 

ve  equation  (in.’^). 

DISC 

■.0 

Diu  ping  .-oe' 

i.‘f  '  at  zero  lead  (Ih  "  in.) 

7HSCE 

FIO.O 

Unl.'>ade  .hi 

•  ‘.s  if  cu'bk  nfe  nitttxrks  (in.) 

LINE  9:  Back  Cusliion  Properties 


DESCRIPnON:  Force-deflection  characteristics,  damping,  and  thickness  for  back  cusiiion. 

Tliese  characteristics  should  be  measured  using  an  indentei  with  the  fomi  of 
the  occupant  torso.  If  a  dummy  torso  is  used  in  measurement,  the 
deflection  should  be  based  on  the  ‘best  accelerometer  location.  The  force, 
F,  is  computed  from  cushion  deflection,  5,  according  to  F  =  C(e^^  -  1). 

FORMAT  AND  RXAlVfPI  E: 


1 

2 

4 

5 

6 

7 

T2545§7S90 

T^156:TO 

1234567890 

1234567896 

1234567890 

1234567^ 

12^4567890 

msmm 

J3BC 

DPBC 

THBCE 

760. 

0,50 

■HKEE 

EEI.Q 

e:2rmm 

CONTENTS 

CBC 

FIO.O 

Coefficient  C  in  above  equation  (lb). 

BBC 

FI  0.0 

Coefficient  B  in  above  equation  (in.'^). 

DPBC 

FIO.O 

Damping  coefficient  at  z£:ro  load  (Ib-secAn.) 

THBCE 

FIO.O 

Unloaded  thickness  of  cushion  in  center  of  seat  back  (in.). 

A  r.’ 


LINE  10:  Headrest  Cushion  Properties 


DESCRIPTION:  Force- deflection  characterisdes,  damping,  and  thickness  for  headrest 

cushion.  The  measuretnent  should  be  made  using  a  headfonn,  'Fhe  force, 
F,  is  computed  :*rom  cushion  deflection,  f>,  according  to  F  "  Cfe®*' -  1).  Jf 
CliK  "  0  (or  blaitk),  the  headrest  is  omitted  from  the  scat  configuratton. 

FORMiVT  .\ND  EX.VMPLE: 


—J 

"  - ^ 

- 5 

- T"  5 

6 

- - - 

1154567890 

155^675^ 

:i5456'.'f^9& 

— mR~'^ 

BHP  :  DPHIl 

ITUlRfe  ! 

..  _ ] 

"7^ 

a.q 

1 

mui 

E2RMAI 

CQNIRKTS, 

CHR 

FiO.O 

Coefficient  C  in  above  e/iuation  (lb). 

BHR 

FIO.O 

Coefficient  B  in  above  ec^uaticn  (in.’S- 

DPfiR 

FIO.O 

Damping  coefficient  (Ib-secAn.)- 

THHRE 

F’0.0 

Unloaaed  thickness  of  cushion  behind  head  (in,). 

line  11:  Lap  Belt  Pniperties 


LESCRlPnON:  Tables  of  forces  and  deflections  define  an  approxirjiation  to  force-deflection' 

curve  by  tliree  linear  segments,  as  illustrav.J  in  !'igiir>:  A-3.  Tbc  force  ij'nd 
deflection  at  point  1  are  assumed  to  be  zero. 

E^EMAILAmMAMELE. 


1 

— ? 

3 

.  "I 

. 5 

^  7) 

■11545571590 

1^34567890 

T5545e'7g90 

12345671590 

n'534567’590|T234j67S9r)i 

FFlB(2r 

fflb(3) 

fTT,B(4) 

DDLB(f)'^ 

DDLB(.';)  !  1 

1500: 

i  "  ^^.56. 

“TOita 

0:T(T4I5 

O-TbilT  1 

FIEIX) 

FORMAT 

CONTENl’S 

FFLB 

3F10,0 

Forces  (lb). 

DDLB 

3F10.0 

Strains  corresponding  to  forces  (’ 

1 


Force 


Figure  A-.'F  Forcc-deriectio'i  modr!  forrcs'ruint  system  wcbftitig. 
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LINFi  .i2A:  Lap  Belt  Anchor  Points  atid  Footrest,  Passengei'  No.  1 

DESCRIPnON :  Coordinates  of  right  and  left  lap  belt  anchor  points  iu”  aircraft  coordinate 

system. 


FORMAT  AND  EXAMPLE: 


i_ ij 

■■HK 

r"  3"  ■■  1 

6 

L _ 2i 

mmm 

1234.567890  1234567890 

12?456789() 

4734.5078901 

1234567890j 

i  XLB(l,t) 

YLB(1,1) 

\  OB  (2,1) 

ZLBQJ) 

j 

12.51  -30.q  i5.()|  ,12.5 

-lO.G 

_ I 

EIELQ 

imMAi 

XLB(1,1) 

YLB(1,1) 

ZLB(I,1) 

3F10.0 

Coordinates  of  right-hand  lap  bel*  anchot  point  in  aircraft  coordinate 
sy.stern  (in.) 

XLB(2,1) 

YLa(2,l) 

ZLB(2,1) 

3F’0.0 

Gxirdinates  of  left-hand  lap  belt  ar'chor  point  in  aircraft  coordinate 
system  (in.). 

LItffi,  jlLV  Lap  Bell  Arichor  Feints  and  Footrest,  Passenger  No.  2' 

DESCRiFn'ON :  G-omior'scs  of  righ‘.  and  lef.  'ap  belt  anchor  poinr;  in  aircraft  c<M>rdinate 

svst.’m 


4i  5 

7 

11 i  I5^h7;7B^lB45^/S5o 

T254?57S^  i5W<57^0 

15335^®’^ 

rXLB(Tiri  Yl^lF^)  f  ZLBt’l.f) 

XLB(2,^.) 

YLB0,2) 

zui(I^Tl 

r  il3r  "HB.q 

"~n3 

1(53 

I5.q 

„ 

aiiU?. 

K>S^:<AI 

.CQiCL£m:s 

XLB(i,2) 

3F10.0 

L'oordivnates  of  right-har.o  lap  belt  anchor  point  in  aircraft  coordinate 

YLB(?,2) 

ZL3(1.2) 

system  (in.). 

XLB(2,2) 

3F10.() 

Coordinaicr  of  left-hand  lap  belt  anchor  point  in  aircraft  coordinate 

YLB(2,2j 

ZLB(2.2) 

system  (in.X 

InrludiMi  (^n!v  i!  '  !  Ylh-  >  I  lasd  ISF^ATu’;  I  <>'!  l.inr 


A .  I  7 


LINE  12C:  Lap  Belt  Anchor  Points  and  Footrest,  Passenger  No.  3* 

DESCRIPTION:  Coordinates  of  right  and  left  lap  belt  anchor  points  in  aircrai't  c(K)rdinalc 

system. 


FORMAT  AND  EXAMPLE: 


FIELD 

FORMAT 

CON'IENTS 

XLB(1,3) 

YLB(1,3) 

ZLB(1,3) 

3F10.0 

Coordinates  of  right-hand  lap  belt  anchor  point  in  aircraft  coordinate 
system  (in.). 

XLB(2,3) 

YLB(2,3) 

ZLB(2,3) 

3F10.0 

Coordinates  of  left-hand  lap  belt  anchor  point  in  aircraft  coordinate 
system  (in.). 

‘  Included  only  it  Id'Vp!-  i  jnd  I.SI  'A  ILL  I  u!i  I  .ine  1 


A  IH 


LINE  13:  Shoulder  Belt  Pro[)erties  (used  only  if  IRSYS  >  ())* 


DESCRIPTION:  I'ables  of  forces  and  deflections  define  an  approximatitjn  to  force-  deflection 

curv'e  by  three  linear  segments,  as  illustrateti  in  Idgurc  A-3.  The  force  and 
deflection  at  j)oint  1  arc  assumed  to  be  zero. 

mfil'IAXANDJiXAJdELE: 


r  I 

5 

7 

mmm 

mmm 

immm 

1 

I 

1 

msmm 

FFSH(2) 

FFSH(3) 

FFSH(4) 

DDSH(2) 

DDSHO) 

DDSH(4) 

FIELD 

EQBMAI 

mfflENIS 

FFSH 

3F10.0 

Forces  (lb). 

DDSH 

3F10.0 

Strain  corresponding  to  forces  (in./in.). 

*  Not  used  ill  sample  case,  thcrefoic  blank. 


A  i') 


LiNfiJiAx-LlLi-MC:  Sfiouldcr  Belt  Anchor  Points  (usct!  only  if  IRSYS  >  ())  > 
DESCRIFriON:  Coordinates  of  shoulder  l>clt  anchor  point  in  aircraft  ctKtrtlinate  system. 


FORMAT  AND  EXAMPLE: 


HHHHD 

3 

'I 

1  7 

1234.567890 

1234567890 

123456789011234567890 

XSH(l) 

YSH(I) 

ZSH(l) 

BUKL(l) 

XTRAL(l) 

XSH(2) 

YSH(2) 

BUKL(2) 

XTRAL(2) 

XSH(3) 

YSH(3) 

ZSH(3) 

BUKiA3) 

X'TRAL(3) 

-  - — 

field 

FORMAT 

CONTliNTS 

XSH(l) 

YSH(l) 

ZSH(l) 

3F10.0 

Coordinates  of  shoulder  belt  anchor  point  in  aircraft  coordinate 
system,  or  point  from  which  belt  passes  to  shoulder  in  a  straight 
line. 

BUKL 

FI  0.0 

Ixngth  of  lap  belt  webbing  attached  to  buckle,  as  illustrated  in 
Figure  A-L 

XTRAL 

FT  0.0 

Length  of  shoulder  strap  beyond  (XSH,  YSH,  ZSH)  if  strap  is  not 
in  straight  line  from  anchor  point  to  shoulder,  as  shown  in  Figure  A- 
4  (not  used  if  IRSY'^S  =  0). 

Lines  14B  and 

14C  repeat 

Line  14A  for  passengers  2  and  3;  they  are  included  only  if  NOCC>l. 

Not  used  lii  s'liiijiie  ; Mst.  ,  t.hervfore  Idank. 

A  M) 


LINE  15:  Tiedown  Strap  Propciries  (used  only  if  IRSYS  -  4)* 

DESCFJR''IQN:  Table  of  forces  and  deflections  define  an  approximation  to  idree -deflection 

curve  by  three  linear  segnienis,  as  illustrated  in  Figure  A  •3.  The  force  and 
deflection  at  point  1  are  assumed  to  be  zero. 

ECMMLMD£XAMPLE: 


1  .  ^ 

4 

r“  5 

6 

7 

i 

I 

1234567890 

1234567890 

1 234567890 

1234567890 

FFTD(2) 

FFfD(3) 

FFTD(4) 

DDTD(2) 

DDTi:)(4) 

J 

BBHB 

FIELD 

FDRMaXT 

CONTENTS 

FFTD 

3F10.0 

Forces  (lb). 

DDTD 

SFIO.O 

Strain  con'esponding  to  forces  (in./in.). 

# 


'  N('t  used  m  Aaiiinlc  vMS'*.  ilurw* 


HV 


LINE  16A.  i6B.  16C:  Tiedown  Strap  Anchor  Points  (used  only  if  IRSYS  =  4)* 

DESCRII^TON:  Ccwrdinates  of  lap  telt  tiedown  strap  anchor  points  in  aircraft  coordinate 

system. 


r 

iiiiiiiiiiiiiiliHiiiHQji^^ 

6 

7 

mmm 

mmm 

1 

i 

lI34557gW)IT2345678^{5 

"1234567890 

1234567890 

XTD(l) 

YTDO) 

ZTD(l) 

r  YTi5{2r~ 

trDi'ij 

xtd(5) 

YTD(3) 

ZTD^) 

EIHJ2  FORMAT  CONTENTS 

XTOT)  3F10.G  Coordinates  of  right-hand  lap  belt  anchor  point  in  aircraft  coordinate 

'VTD(l)  system  (in  ). 

ZrD(l) 

Lines  16B  and  16C  repeal  Line  i6A  for  passengers  2  and  3;  they  are  included  only  if  NOCC>l. 


'No:  u\c(!  in  sai  'Ir  c  ise,  ii'Crfirta'  'liank 


A 


•>  A 


LINE  17:  Additiofud  Belt  Projjeities 

DESCEIITIQN:  Damping  coefficient  and  belt  slack  for  lap  belt,  shoulder  belt(s),  and 

tiedown  strap. 

FORMA  T  AND  EXAMPLE: 


1 

“T  "  .3 

4 

5 

7 

[1234567850 

1234567890 

1234^890 

DPLB 

SLAB 

DPSH 

DPIT) 

SShhH 

■hBBh 

FIELD 

FORMAT 

CONTENTS 

DPLB 

FI  0.0 

Lap  belt  damping  coefficient  (lb-sec). 

SLAB 

FI  0.0 

Lap  belt  slack  (in,). 

DPSH 

FIO.O 

Shoulder  belt  damping  coefficient  (Ib  -sec). 

SLSH 

FIO.O 

Shoulder  belt  slack  (in.). 

DPIT) 

FIO.O 

Tiedown  stop  'damping  ci>effic:ent  (lb-sec). 

SLTD 

FIO.O 

Tiedown  strap  slack  (in.). 

LINE  18:  Other  seating  arid  restraint  data 


DESCRIPTION :  Friction  coefficients  and  footrest  location. 


FIELD 

FORMAT 

CONTENT'S 

COEIffiS 

FIO.O 

Seat  cushion  friction  coefficient. 

COEFFR 

FIO.O 

Floor-foot  friction  coefficient. 

XFR 

FIO.O 

X-coordinate  of  footrest  in  aircrafi  cixndinate  system,  at  intersectioii 
with  floor,  where  Z  -  0. 

ANGFR 

FIO.O 

Angle  between  footrest  and  floor  in  degrees. 

LINE  19:  Aircraft  Initial  Position 


DE^jCRIPIlON:  Components  of  aircraft  initial  position,  in  earth-fixed  coordinate  system, 

and  attitude. 


EmMAT  AND  EXAfdPLF: 


HEIX)  FORMAT  CQNT'ENT’S 


XA 

YA 

ZA 


3F10.U  Position  of  aircraft  coordinate  systein  in  inertial  system  (earth-fixed 

system  in  which  gravity  acts  in  the  -Z  direction)  (in.).  These  initial 
coordinates  are  normally  taken  tis  (0 , 0.,  0.)  unless  displacement 
iron'  a  specific  point  is  desired.  For  example,  if  the  simulation  is  to 
be  in' tinted  at  some  horizontal  distance  from  a  barrier,  such  as  bO 
in.,  in.  initial  position  could  be  specified  as  (-60.,  0.,  0.).  If  the 
simulation  is  (o  Iregin  10  in.  above  tr.e  ground  in  a  vertical  drop,  the 
initial  position  could  be  specified  as  (0.,  0.,  10.).  These  crxirdinates 
are  not  used  in  the  simulation  bur  only  in  output  of  aircraft  position. 


YAW 

PITCH 

ROU. 


31  10.0 


Initial  attitude  of  aircraft  relative  to  earth-fixed  system  (deg). 


LML2!]:  ALvcraft  Initial  Vet  xity 


D>^SCRIPTION :  Components  of  aircraft  initial  velocity,  in  aircraft  cooidinate  system, 

traits!  ation  and  rotation. 


\'X  3F10.0  Components  of  aircraft  initial  velocity  in  aircraft  coordinate  system 

VY  (ft/set;). 

VZ 


DYiVW  3FiO.O  Yaw,  pitch,  and  roll  rates  (rad/sec). 

DPITCH 

DPax 


Lih'E  21:  Aircraft  Acceleration 

DHSCRIITION:  The  time  v«i  lotion  of  tne  i-ix  cornpenents  of  the  acccleraticii  of  the  aircraft 

coordinate  system  is  approximated  by  up  to  40  points  in  acceleration  and 
time.  NIMPT  lines  must  be  included  (up  to  40). 

PQRMAT  AND  EXAMPLE:  (4  lines  for  NIMPT  -  4  on  line  .r) 


- 2 

5 

f"  7| 

fmmmi 

immm 

wmsm 

fl 

mmm 

\mmsm 

1  TA 

AX 

AY 

.AZ 

AVAW 

APiT 

ARQL  ! 

r  “TTs 

mmmm 

(Hi 

iKO 

1 — -  ■ '  — TS’tJ 

r 

-6.D 

1 

O.Oi 

: 

-GiTT) 

^  T)t1 

HHKSwl 

_ M _ ^ 

o.q  O.Oi 

HMKIlt 

r 

L- 

(TO 

HELD 

FGRMAT 

TA 

FI  0.0 

AX 

Pio.o 

AY 

FI  0.0 

/kZ 

FlOO 

AYAW 

FI  0.0 

APIT 

FI  0.0 

/tR-OL 

FIO.O 

C:QN"rENTS 
Time  (sec), 

X  -acceleration  (G). 
Y-accckration  (C). 
Z-acccIeraiion  (G). 

Yaw  iccelcratior.  (ratVsec/setO- 
Pitch  acceleration  (rad/scc/scc) 
Roll  acceleration  (rad/'seo’scc). 


ONE  22:  Occupant  iniii..,'  ‘‘osition,  Psssenger  No.  I 

■  Initial  pcshion  angles  and  heel  X- position,  as  illustrated  in  Figure  A-5.  The 
heels  are  assumed  to  begin  at  Z  =  0  The  torso  is  aligned  according  to 
GAM(1,1),  GAM(2,1).  and  GAM(3,1),  and  the  position  is  then  determined 
from  static  equilibrium,  allowing  for  compression  of  the  cushions.  Also, 
the  Y-r.oordinatc  of  the  occupant  plane  ’  symmetry  is  included.  (Line  22 
must  be  included  for  each  occupant,  three  in  this  example.) 

mKMAIAC!a2£2SLAMmi: 


rzzzi 

4 

- § 

- - r  -1 

uamtm 

GAM(2.1) 

GAMO.l) 

gf^arEim 

)mEEL(l)|  YPaW)! 

FORMAT 

—11  i  1..  ■■■— .1  -  i.i  1  .—I 

mNTEm 

GAM(I,1) 

5F10.0 

Vector  of  initial  position  angular  coordinates,  as  illustrated  in  Figure 
A-5  (deg). 

XHEELd) 

FIO.O 

X-coordinate  of  heels  in  aircraft  coordinate  system  (in.). 

YPASS(l) 

FIO.O 

Y-coordinatc  of  mid-plane  (plane  of  symnoetry)  for  occupant  (in.). 

LINE  22:  Occupant  Initial  Posiiion,  Passenger  No.  2  (included  only  if  NOCC>l) 


i 


D£SCRrPn(,>N:  Initial  ptxsition  angles  and  heel  X -position,  as  illustrated  in  Figure  A-5.  The 

heels  are  assumtid  to  begin  at  Z  -  0.  The  torso  is  aligned  according  to 
GAM(1,2),  GAM(2,2),  and  GAM(3,2),  and  the  position  is  then  determined 
from  static  equilibrium,  allowing  for  compression  of  the  cushions.  Also, 
the  Y'Coordinate  of  the  occaptint  plane  of  symmetry  is  included. 


1 


“TSIJ 


UMTCZir!  GAM(3,2) 


7l! 


■S:am(4,2) 


TO 


TaM(5,2) 


* - 6f 

1154557^5^ 


XHEEL(2) 


Tl345^789() 


YPASS(i) 

- cm 


FIELD 

FORMAT 

contents 

GAMa,2) 

5F10.0 

Vector  of  initial  position  angular  coordinate,  as  illustiated  in  Figure 
A-5  (deg). 

XfUiEL(2) 

FT  0.0 

X-coordinate  of  heels  in  aircraft  coordinate  system  (in.). 

YPASS(2) 

FIO.O 

Y-coo;dinate  of  mid-plane  (plane  of  symmetry)  for  occupant  (in.). 

LINE  22:  Occupant  Initial  Position,  Passenger  No.  3  (included  only  if  NOCC  =  3) 


DESCRIPTION:  Initial  position  angles  and  heel  X-^wsition,  as  illustrated  in  Figure  A-.5.  The 

heels  are  assumed  to  begin  at  Z  =  0.  The  torso  is  aligned  according  to 
GAM(1,3),  GAM(2,3),  and  GAM(3,3),  and  the  position  is  then  determined 
from  static  equilibrium,  allowing  for  conripression  of  the  cushions.  Also, 
tlie  Y-coondinate  of  the  occupant  plane  of  symmetry  is  included. 

FORMAT  AND  EXAMPLE: 


1 

.  3 

4 

5 

7 

1254567S90I 

1234567890 

mmmsi 

1234567890 

1234567890 

1234567890 

1234567890 

dAM(l,5) 

GAM(2,5) 

GAM(3,3) 

GAM(4,3) 

GAM(5,3) 

XHEEL(3) 

YPASS(3) 

-16.0 

wmmm 

32.0i 

FIELD 

FORMAT 

CONTENTS 

GAMa,3) 

5F10.0 

Vector  of  initial  position  angular  ctxjrdinate,  as  illustrated  in  Figure 
A-5  (deg). 

XHEEL(3) 

FIO.O 

X-coordinate  of  heels  in  aircraft  coordinate  system  (in.). 

YPASS(3) 

FIO.O 

Y-coordinate  of  mid-plane  (plane  of  symmetry)  for  occupant  (in.). 

A  >2 


Nonstandard  Occupant  Input  Data* 


If  nonstandard  occupants  are  requested  by  setting  IMAN  =  2  (human)  or  IMAN  =  3  (dunnn>)  on 
Line  3,  then  12  additional  lines  must  be  inserted  for  each  occupant.  The  format  for  these  12  line'., 
referred  to  as  22A  -  22L,  is  explained  on  the  following  15  pages.  If  IMAN  =  0  (standard  50th 
percentile  human)  or  IMAN  -  1  (standard  50tli~percentile  dummy),  skip  this  section  and  proccfxi  ro 
Line  23. 


‘  Fhcsc  Inscs  nuis;  hi-  p/  ivulcii  !or  each  occupan:  after  each  i.inc  .’2  ''['>ci.  il\ liig  ia'i'cspi'nilirai 
tvcu[iant  iiiiiia!  [Hcsiuot!  They  arc  not  inciiulcii  its  the  sample  c.'ise,  but  an  example  is  I'lrcvkiei!  in 
AjifKavdis  H. 


Alt 


LME_22A:  Scgiiient  Lengths  (otily  if  IMAN  =  2  or  3) 

DESCRIPTION:  lengths  of  the  spine  and  segments  3,  4,  5,  8,  and  9  as  described  in  Figure 

A-6.  llie  lengths  of  segments  6,  7,  10,  and  1 1  are  obtained  from  these  by 
symmetry  (m.). 

FORMAT  AND  EXAMPLE: 


1 

99B9B 

5 

6 

7 

T^56785(r 

mmMm 

1234567850 

T2^45'S7895 

SPL 

XL(3) 

XL(4) 

XL(5) 

XL(8) 

XL(9) 

r 

FIELD 

FORMAT 

CONTENTS 

SPL 

FI  0.0 

Spinal  length. 

XI.(3) 

FI  0.0 

Head  length. 

Xl.(4) 

FI  0.0 

Upper  ami  length. 

XL(5) 

FIDO 

I^wer  arm  length  -  elbow  to  mid-point  of  hand. 

XI..(8) 

FI  0.0 

Upper  leg  length. 

XI,  (9) 

FIO.O 

1  ower  leg  length  -  knee  to  ankle. 

A  34 


^  Center  of  Mass 
#  Joint 

:IID  Beam-Column  Element 


{■ii'iirt*  A  6. 


B(xU  segiiK-nl  diiiu’n.'  icvfi 


LINE  22B:  Segmeni  Center  of  Mass  fxication  (only  if  IMAN  =  2  or  3) 

DESCRIFFION:  Cent«-r  of  mass  locations  for  segments  1,  2,  3,  4,  5,  8,  and  9.  See  Figure 

A-6  for  datum  plane  description  (in.)- 


r 

p- - 2 

3 

5 

7 

■n;545^7890 

. 

1234567550 

12345(j7850 

1234567890 

RHO(l) 

RHO(2) 

HHO(3) 

RHO(4) 

RHO(5) 

RHO(8) 

RHO(9) 

RELD 

FORMAT 

'.20NTEr4TS 

RHO(l) 

PIO.O 

l^wer  torso  center  of  mass  vertical  distance  from  hip  pivot. 

RHO(2) 

FIO.O 

Upper  torso  center  of  mass  distance  from  base  of  neck. 

RHO(3) 

FI  0.0 

Head  center  of  mass  distance  from  base  of  neck. 

RHO(4) 

FIO.O 

Upper  ami  center  of  mass  distance  from  shoulder  pivot. 

RFIO(5) 

FIO.O 

Lower  arm  center  of  mass  distance  frons  elbow  pivot, 

RF?0(8) 

FIO.O 

Upper  leg  center  of  mass  distance  from  hip  pivot. 

RHO(9) 

FIO.O 

Lov;er  leg  center  of  mass  distance  from  knee  pivot. 

A  36 


LINE  22C:  Segment  Weight  (only  if  IMAN  =  2  or  3) 


Weights  of  segments  1 , 2,  3,  4,  5,  8,  and  9  (lb). 


r:  1 

~j - 3 

4 

7 

1 

i 

1 

1 

1 

''11545^78^ 

nmomi 

SW(4) 

isyfeM 

SW(8) 

IHHH 

FIELD 

R)Ri>/IAT 

CON’rENTS 

SW(1) 

FI  0.0 

Lo  Afcr  torso  weight. 

SW(2) 

FIO.O 

Upper  torso  weight. 

SW(3) 

FI  0.0 

Head/neck  weight. 

SW(4) 

FiO.O 

Upper  arm  weight. 

SW(5) 

FIO.O 

Lower  arm  weight. 

SW(8) 

FIO.O 

Upper  leg  weight. 

SW(9) 

rlO.O 

Lower  leg  weight. 

A  v; 


LINE  22D:  Segment  Moment  of  Inertia  with  Respect  to  Ixtca!  x-axis  (only  if  IMAN  =  2  or  3) 


DESCRIPTION :  Moments  of  inertia  with  respect  to  x-axis  for  segments  1,  2,  3,  4,  5,  8,  and 

9  (lb-in. -sec^). 

FORMAT  AND  EXAMPLE: 


1  _ 

"  "  ^ 

3 

4 

5 

6 

I 

I 

1 

T2345S7M) 

C1X(2) 

CIX(3) 

CIX(4) 

CIX(.5) 

CIX(8) 

CIX(9)  i 

iSfiSH 

! 

r  1 

FIELD 

FORMAT 

CONTENTS 

CIX(l) 

El  0,0 

Lx)wer  torso  x-axis  moment  of  inertia. 

CIX(2) 

FIO.O 

Upfxir  torso  x-axis  moment  of  inertia. 

CIX(3) 

FIO.O 

Head/neck  x-axis  moment  of  ineitia. 

CIX(4) 

FIO.O 

Uppter  ajTn  x-axis  moment  of  inertia. 

ClX(5) 

FIO.O 

Ix'wer  aiTT!  x-axis  moment  of  inertia. 

CIX(8) 

FIO.O 

Up;x;r  leg  x-axis  moment  of  inenia. 

CIX(9) 

FIO.O 

lx)wer  leg  x-axis  iiKwnen*  of  itsertia. 

LINE  22E:  Segment  Moment  of  Inertia  with  Respect  to  L-ocal  y-axis  (only  if  IMAN  =  2  or  3) 

DESCRimON:  Moments  of  inertia  with  respect  to  y-axis  for  .segments  1,  2,  3,  4,  5,  8,  and 

9  (lb-in. -sec^). 


FORMAT  AND  EXAMPLE: 


- 2 

3 

4 

5 

^  _  7J 

mmm 

mmmm 

1 

1 

1 

1 

s 

1 

1 

i 

CIY(l) 

nww 

BBS 

i 

1 _ _] 

1 

i _ 

FIELD 

FX3RMAT 

CONTENTS 

CIY(l) 

FiO.O 

Lower  torso  y-axis  moment  of  inertia. 

aY(2) 

FIO.O 

Upper  torso  y  -axis  moment  of  inertia. 

CIY(3) 

FIO.G 

Head/neck  y-axis  moment  of  inertia. 

C1Y(4) 

FIO.O 

Upper  ann  y-axis  moment  of  inertia. 

aY(5) 

FIO.O 

Lower  arm  y-axis  moment  of  inertia . 

CIY(8) 

FIO.O 

Up{>er  leg  y-axis  moment  of  inertia. 

C1Y(9) 

FIO.O 

lx)wer  leg  y  axis  moment  of  inertia. 

LINE  22F:  Segment  Moment  of  Inenia  with  Resjject  to  Local  z-axis  (only  if  IM  AN  -  2  or  3) 


DESCRlEnQN :  Moments  of  inertia  with  respect  t(>  z-axis  for  segments  1,  2,  3,  4,  5.  8,  and 

9  (lb-in. -sec^). 

FORMAT  AND  EXAMPLE: 


1 

3 

mniiMK 

7 

1234567^1)' 

immm 

1234567890 

mmm 

CIZ(l) 

az(2) 

CIZ(4) 

C!Z(5) 

CIZ(9) 

_ 

_ 

I3£LQ 

PORM/tT 

a)Mnii:n:s 

CIZ(l) 

FI  0.0 

Lower  torso  z-axis  moment  of  inertia 

CIZ(2) 

FI  0.0 

Upper  torso  z-axis  moment  of  inertia. 

CIZ(3) 

FI  0.0 

Head/'neck  z-axis  tnonienf  of  inertia. 

CIZ(4) 

FIO.O 

Uppt^r  ann  z-axis  moment  of  inertia. 

CIZ(5) 

FI  0.0 

Ix)wer  amt  z-axis  nxtrrtent  of  inertia. 

CIZ(8) 

FIO.O 

Upper  leg  z-axi.s  ittoment  of  inertia. 

CIZf9) 

FIO.O 

.'.xtwer  leg  /.  axis  moment  of  inenia. 

A  Ul 


LINE  22G:  Contact  Smface  Radii  (only  if  IMAN  -  2  or  3) 


DESCRIPTION: 


Radii  of  contact  surfaces  1,  2,  3,  4,  5,  8,  and  9  (in.).  (See  Figure  A-7  and 
Table  A-1  for  human  occupant.) 


.FORMAT  AND  EXAMPUI 


- 2 

.nnnngiiiiiig 

:  ”1 

1 

1 

I 

1 

1 

massmm 

XR(I) 

XR(2) 

XR(3) 

XR(4) 

XR(5) 

XR(9)  1 

L. 

„J 

HELD 

FORMAT 

CON^fENl’S 

XR(1) 

FiO.O 

Radius  of  lower  torso  contact  .surface  ellipsoid 

XR(2) 

FI  0.0 

Radius  of  upper  torso  in  mid-saggital  plane. 

XR(3^ 

FIO.O 

Radius  of  head  in  mid-saggital  plane. 

XR(4) 

FIO.O 

Radius  of  upper  arm  contact  surface  cylinder . 

XR(5) 

FIO.O 

Radius  of  lower  arm  contact  surface  cylinder. 

XR(8) 

FIO.O 

Radius  of  upper  leg  contact  surface  cylinder. 

XR(9) 

FiO.O 

Radius  of  lower  leg  contact  surface  cylinder. 

ire  A  7.  Body  contact  surfaces  descrintion. 


TABLE  A-1.  SI’ANDARD  CONTACT  SURFACE  DIMEMSIQM 


Surface 

SyoiM 

Fraction 
of  Stature 

iRm _ 

Actual  Dimension 
for  50th- 

Percentile  Human 

MakiiiLl..  _ _ 

Pelvis 

Rl 

(1.0579 

4.00 

Chest 

Rl 

0.0689 

4.76 

Head 

R3 

0.0485 

3.35 

Ann 

R4,R6 

0.0263 

1.82 

Forearm 

Rs.R? 

0.0243 

1.68 

Tnigh 

Rs.Rio 

0.0466 

3.22 

Ixg 

R9.R11 

0.0344 

2.38 

Knee 

Ri2,Ri3 

0.0373 

2.58 

Foot 

R14.R15 

0.04C5 

3.10 

Hip 

R16,Ri7 

0.0515 

3.56 

Shoulder 

Rl8.f^l9 

0.0378 

2.61 

Elbow 

R20.R2  i 

O.O208 

1.85 

Hand 

R22.R23 

0.0339 

2.34 

LINE  22H:  Contact  Surface  Radii  (i!ontirnied  (only  if  IM AN  =  2  or  .1) 
DESCRIP110N:  Radii  of  contact  surfaces  12,  14,  16,  18,  20,  and  22  (in.). 
FORIrfAT  AND  EXAMPLE: 


2 

4 

5 

S 

7 

mmm 

1234567890 

1234567890 

1234567890 

1234567890 

1234567890 

1234567890 

XR(12) 

XR(14) 

XR(16) 

XR(18) 

XR(20) 

XR(22) 

FIELD 

FORMAT 

CONTENTS 

XR(12) 

FiO.O 

Radius  of  nec>  contact  surface  ellipsoid. 

XR(14) 

FIO.O 

Radius  of  f<x)t  contact  surface  sphere. 

XR(16) 

FIO.O 

Radius  of  hip  contact  surface  sphere. 

XR(18) 

FIO.O 

Radius  of  shoulder  contact  surface  sphere. 

XR(20) 

FiO.O 

Radius  of  elbow  contact  surface  sphere. 

XR(22) 

FIO.O 

Radius  of  hand  contact  surface  sphere. 

A  14 


LINE  221:  Spherical  Joint  and  Center  of  Mass  Offset  Distances  (only  if  IMAN  -=  2  or  3) 


Distances  that  spherical  joints  (shoulder  and  hip)  are  laterally  offset  from  the 
mid-saggital  plane,  and  the  anterior  offset  of  the  major  upper  body  segment 
(lower  torso,  upper  torso,  and  head)  center  of  masses  from  the  spine.  (See 
description  of  distances  in  Figure  A-6.)  Dimensions  are  in  indies. 


mRMAIANIIEMMPLE: 


4 

5 

1 

1 

1 

msmmsM 

XLH 

mmmM 

EM(1) 

EM(2) 

EM(3) 

'jjlllllHBIIHHHHBHi 

j 

FIELD 

TORMAJ' 

CONTE.NFS 

XLH 

FI  0.0 

Lateral  distance  of  center  of  hip  joint  from  mid-saggital  plane. 

XLS 

FIO.O 

Lateral  distance  of  shoulder  joint  from  mid-saggital  plane. 

EM(1) 

FIO.O 

Anterior  offset  distance  of  the  lower  torso  center  of  mass  from  the 
spine. 

EM(2) 

FIO.O 

Anterior  offset  distance  of  the  upper  torso  center  of  mass  from  the 
spine. 

EM(3) 

FIO.O 

Aiiterior  offset  distance  of  the  head  center  of  mass  from  the  spine. 

LINE  22J:  Abdomen  and  Chest  Compliance  (otily  if  IMfAN  =  2  or  3) 


DESCRII^TON:  Estimated  force-deflection  characteristics  (compliance)  ol'(x;cupant  chest 

and  abdomen  under  restraint  system  loads.  'Ehe  force,  I-,  is  computed  from 
cushion  deflection,  6,  according  to  F  =  C(e*^^  -  1 ). 

FORMAT  AND  EXAMPLE; 


1 

3 

4 

5 

6 

7 

1234567890 

1234567890 

1234567890 

1234567890 

(234567890 

1234567890 

1234567890; 

CARD 

BABD 

BCHE 

-  - i 

FIELD 

FORMAl’ 

CONTENTS 

CARD 

FIO.O 

C(5efficic.it  C  for  abdomen  compliance  (lb). 

BABD 

FI  0.0 

Coefficient  B  for  abdomen  compliance  (in. 

CCHE 

FIO.O 

Coefficient  C  for  chest  compliance  (lb). 

BCHE 

FIO.O 

Coefficient  B  for  chest  compliance  (in.’ V 

A  4( 


I.JNII  22 (C:  Axiid  SuiTne!;s  and  Damping  Properties  for  Spine  and  Ncc!;  (only  if  IMAN  -  2  or  3) 


DESCRIFHQN:  Axiai  forcc-tieflection  charactcrist'cs  for  tLe  S|iine  and  neck  beam  nicxicls 

and  associated  axial  damping.  'iTie  force,  F.  is  computed  from  deflection, 
5,  according  to  F  -  -  >.). 

FORMAT  AND  EXAMPU:: 


1 

21 

■HfllllllllllK 

4 

7 

TI345S785() 

T^456789(^ 

1234567890 

i2345678''{) 

T254567W 

CAXS 

BAXS 

DMPS 

CAXN 

BAXN 

DMPN 

J 

_ _ _ 1 

tJELL) 

EQEMAT 

CAXS 

FiO.O 

BAXS 

FIO.O 

DMPS 

FIO.O 

CAXN 

FIO.O 

BAXN 

FIO.O 

DMPN 

FI  0,0 

Coefficient  C  in  above  equation  for  axial  spinal  stiffness  (lb). 
Coefficient  B  in  above  equation  for  axial  spinal  stiffness  (in.  A). 
Axial  damping  in  spine  (lb  sec-in. '!). 

Coefficient  C  in  above  equation  for  axial  neck  stiffness  (lb). 
Coefficient  B  in  above  equation  for  axial  neck  stiffness  (in.‘  ^). 
Axial  damping  in  neck  (lb-sec  in." ^). 


I..1NH  221-.:  Rotational  StiiTncss  and  Dampnig  !Vofv?nies  for  Spine*  and  Neck  daily  if  IMAN  ^  2  or 
3) 

pn^tClllFTION:  Rotational  tnonicnt-angk".  characteristics  for  the  spiinc  and  neck  beam  mcxiels 

and  associated  rotational  damping.  The  moment,  M,  is  compared  from 
angular  deflectioti,  h,  according  <o  M  ■=  C(e'^^  -  1  k 


FORMAT  AND  EXAMPL.E: 


1 

. ■■■'2 

r  "  4 

5 

6 

7 

1 -234567890 

[T2T456789() 

1234567890 

h2345'678W 

1234567890 

1234567X90 

1234507890 

crTTt(IP 

BROT(l) 

XJ(1) 

CR0T(2) 

BR0T(21 

XJ(2)| 

FORMAi' 

CQMTGlsrrS 

CROT(l) 

FI  0.0 

Coefficient  C  in  above  et-juatioii  for  rotational  .spinal  stiffness 
(in.-'b.). 

EROT(l) 

FI  0.0 

Coefficient  B  in  above  equation  for  rotationai  spinal  stiffness 
(rad'p. 

XJ(!) 

FI  0,0 

Rotationai  damping  in  spine  (lb-sec). 

CROi(2) 

FI  0.0 

Coefficient  C  in  above  aiuation  for  rotational  neck  stiffness 
(in. -lb). 

BROT(2) 

FI  0.0 

Coefficient  B  in  abtjve  equation  for  rotational  neck  stiffness 
(radp. 

XJ(2) 

FI  0.0 

Rotational  damping  in  neck  (lb  sec). 

,A 
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PEgCRlPriON •  Din-^ensions  of  seac  model  as  shown  in  Figure  A-8. 
FX)RMA1^  AND  EXAMPI^: 


_ r 

_ 

4 

5 

J 

rz  1 

iifcgi3gHa;ifl 

12^5678901 

12M567890 

iXSEAT 

ZSEAT 

angSp 

ANGSB 

xlpan 

XWPAN 

SBHT 

i  10.0 

8.0 

_  iOO 

(■■KSiE 

r  20.(1 

FIELD 

FORMAT 

C(0N1EN1'S 

XSEAT 

ZSEAT 

2F10.0 

X-  and  Z.-coordinates  (in  aircraft-fixed  system)  of  intersection  of 
seat  pan  and  seat  back  planes  under  the  cushions  (in.). 

ANGSP 

ANGSB 

2F10.0 

Seat  pan  and  seat  back  angles  (in  aircraft-fixed  system),  directions 
as  defined  in  Figure  A-8  (deg). 

Xl.PAN 

FIO.O 

Seat  pan  length  (in.). 

XV/PAN 

FI  0.0 

Seat  pan  width  (in.). 

SBHT 

FIO.O 

Seat  back  height  (in.). 

A  1<J 


ANGSP 


A.l  ENERGY- ABSORBING  SHAT  INPUT 

Lines  24  through  26  provide  input  for  an  energy-absorbing  seat  option,  which  can  be  used  only  if 
NvSEAT  ~  0  on  IJne  3.  If  the  stroking  seat  weight,  SEATM,  on  Lini  ’4  is  zero  (or  blank),  this 
option  is  not  used  and  any  data  on  Lines  24  Uirough  26  are  ignored.  If  SEATM  is  nonzero,  the 
energy  absorber  force-deflection  data  illustrated  in  Figure  A-9(a)  must  be  provided  on  Line  23.  If 
the  mass  moment  of  ineitia  of  the  seat  with  respect  to  the  Y -axis,  YISEAT,  is  nonzero  on  Line  24, 
then  moment-rotation  data  must  be  provided  on  Line  26. 

If  NSEAT  =  1,  indicating  use  of  the  finite  element  seat  model,  Lines  24-26  are  ignored,  and  seat 
data  continue  with  Line  27. 

Note:  This  example  uses  the  finite  element  seat  model,  so  that  Lines  24-26  are  blank. 

However,  sample  case  no.  3  uses  the  energy-absorbing  seat  option. 


LINE  24:  Energy-Absorbing  Seat  Data 


DESCRIPTION:  Parameters  for  the  two-degree -of-freedom  (scat  stroke  and  rigid-body 

rotation)  energy- absorbing  seat  model.  (See  Figure  A-9  for  a  detailed 
description  of  the  parameters.) 


FORMAT  AND  EXAMPLE: 


1 

2 

3 

4 

5 

6 

7 

1234367890 

1234567890 

1234567890 

1234567890 

1234567890 

1234567890 

1234367890 

SEATM 

SDAMP 

YISEAT 

RUNl.OD 

RDAMP 

FIELD 


SEATM 


ANGEA 


SIINI..OD 
SDAMP 
VISE  AT 


KCNl.Ol) 

RDAMP 


FORMAT  CONTENTS 

FI  0.0  Weiglit  of  movable  par*  of  energy-absorbing  seat  (lb). 

F 1 0,0  Stroking  angle  for  guided  energy-absoibing  seat  (deg),  see  Figure 

A -9b. 

F !  0.0  Energy  absorber  unloading  slope  (Ib/in.). 

F  1 0,0  Damping  CiKdfieieni  tor  tite  energy  absor!x,‘r  (Ib-sec/in.). 

FI  0.0  Mass  moment  of  inertia  of  the  .seat  about  a  lateral  a.sis  through  point 

{,'  (it)  I-'igure  A  <^1  with  coordinates  X  --  XSIiAT,  /.  -  0  (lb-in, - 
see 

F  !().('  I ‘iiload.ng  slope  foi  roialionai  ileformalion  ol  seat  (in.  Ib/rad). 

FIO.O  Rotational  damping  eorlTieiv'iit  lor  the  seat  (in.  lb  see). 


\  ■■■.' 


i2  03 


Stioka  -  In. 


04 


a)  Force-stroke  characteristics 


b)  Geometry 


I  J.NE  25:  iMiet'gy  Absorber  Data 


DlvSCRbPriON:  Energy  absorber  force  versus  (.leflectioii  (iliustralci.!  in  b'igurc  /\  9a). 


1  ^ 

4 

5 

7 

1234567890 

tiaaERtiagii 

1234567890 

1234567890 

1234567890 

1 2.H567890 

FFEA(2) 

FFEA(3) 

FFEA(4) 

DDEA(2) 

DDEA(3) 

f)DHA(4) 

L_  .  ,  ! 

U.— 1.  1  ..1.1. 

FIEL.D  FORMAT  CONTENTS 

F2  3F10.0  Energy  ab.sorbcr  force  (lb). 

F3 

F4 

D2  3F1().0  Deflections  corresponding  to  above  forces  (in.);  see  Figure  A  9a. 

D3 

m 


.\  s.} 


LINE  26:  Rigid  Seat  Rotational  Stiffness  ParaJineters 


DESCRIPTION:  Applied  JBorncnt  versus  seat  rotiitional  angle  as  shown  in  Figure  A- 1 0. 

FORMAT  AND  EXAMPLE: 


r  1 

-  2 

■  '  4 

---  -- 

i234567F;a 

r 23456785011^567^ 

FFRT(2) 

FFRT(3) 

FFRT(4) 

DDRT(2) 

■BSHIl 

FIELD 

FORMAT 

CON'IENTS 

FFRT(2) 

FFRT(3) 

FTRT(4) 

3F10.0 

Applied  moment  on  rigid  seat  (in. -lb) 

DDRT(2) 

DDRT(3) 

DDRT(4) 

3F10.0 

Angulijr  seat  displacement  (rad). 

Moment  -  in.-io. 


Angular  deflection  -  rad 


f  igure  A  H).  Rigid  seat  nKKiei  rota'inrsai  sUfines 


\ 


A.2  NONRIGID  SEA7'  INPUT 

If  a  nonrigicl  seat  is  requested  by  setting  NSEAT  =  1  on  Line  3,  then  the  input  data  described  on 
the  following  lines  is  required  to  define  the  finite  elemcot  seat  model. 


LIMi  27:  Basic  Seat  Model  Data 


DESCRIPTION  :  Control  integers  for  finite  element  model  describing  the  number  of  nodes, 

elements,  materials,  and  cross  sections  in  the  model,  and  the  number  of 
plots. 


FORMAT  AND  EXAMPIT; 


i 

A 

L _ 2 

■nm 

4 

5 

6 

7 

1234567890 

mmm 

1234567890 

1234567890 

1234567890 

1234567890 

iglM 

NUMEL 

INUMAT 

p 

NSFXT 

KSPLT 

wm 

27 

2 

mm 

2 

_ 8 

1 

FIEU> 

FORMAT 

CONTENTS 

NUMNP 

15 

Number  of  real  nodes. 

NTJMEL 

15 

Number  of  elements. 

NUMA'F 

15 

Number  of  materials  (up  to  8). 

NUMDS 

15 

Number  of  displacement-sprecified  node  point:;  (at  which  the  aircraft 
displacement,  vckx:lty,  and  acceleration  aie  applied). 

NCORD 

15 

Numl-ter  of  inactive  beam  pointer  nodes,  which  are  used  to  orient  the 
y-axes  of  beam  cros.s  sections.  A  real  luxle  can  be  used  as  a  prtinter 
node.  Also,  a  single  node  can  be  used  as  a  pointer  node  for  more 
than  one  beam. 

NSECT 

15 

Number  of  different  beam  cross-section  types  (up  to  10). 

NSPLT 

15 

Number  of  requested  ;;eat  position  plots  (up  to  20). 
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LINE  2K:  Miscellaneous  Control  Flags 

DHSCRIPTTON:  Parameters  for  controlling  execution  of  finite  element  seat  simulation. 


EQKMAILAimJEMMPLE: 


i 

1 

5 

61  7] 

1254567^50 

1554567850 

1234567890 

1234567850! 1554567850 

KNTRL 

(0 

KNITU. 

(2) 

5 

5 

EELD.  EQRM^  CONTENTS 


KNTRL(l)  15 


Maximum  number  of  iterations  for  convergence  within  a  time  step 
(default  is  5). 


KNTRL(2)  15  Number  of  increments  to  enforce  the  floor  warping.  (See  Lines  43 

and  44.)  A  value  of  10  is  recommended  for  cases  where  the  floor 
warping  produces  plastic  and/or  large  deformations  of  the  seat 
structure. 


LfNl^  29:  Seat  Plot  Times  and  Viewing  Angles  (number  of  lines  required  -  NSIT  T  on  I  .ine  27) 


Times  when  seat  stmcture  plot  data  are  to  be  stored  on  unit  20,  which  must 
be  saved  as  a  permanent  file  for  subsequent  plotting.  'Hie  elevation  tmd 
azimuth  angles  corresponding  to  each  time  are  illustrated  in  Figure  A- 11. 


FORMAT  AND  EXAMPLE: 


1 

3 

4 

- 

u 

7 

1234567^ 

1234567890 

125^678901 

1234567890 

1234567‘8W 

p2345678'9() 

TSPLT 

THEPLT 

PHIPLT 

o.q 

45.(1 

0.025 

45.C 

20.0 

■■ItnH 

45.C 

20.0 

0.075 

45.0| 

20.0 

0.100 

45.0 

20.0 

0.125 

45.0i 

20.01 

0.150 

45.0 

20.0 

0J75 

45.0 

20.0 

EffiLG 

FORM/V[^ 

CONTF.NTS 

TSPLT 

FIO.O 

Plot  times  (sec). 

THEPLT 

FIO.O 

Azimuth  angle  for  viewing  seat  plot  (deg). 

PHIPLT 

FIO.O 

Elevation  angle  for  viewing  seat  plot  (deg). 

A  N) 


0  1=  Azimuth  angle  in  X-Y  plane  in  degrees  (-180°  <  9  <  +180°) 
(j)  =  Elevation  angle  in  degrees  (-90° <  (|)  <  +90°) 


No<.la!  Output  Selection  (used  only  if  lOU  r(8)  >  0) 

DESCRIPTION :  NtKle  numbers,  in  pairs,  to  specify  v/liich  X,  Y,  Z  dispia'.’ements  an;  to  be 

printed.  (The  siode  numbers  are  defined  on  Line  38.) 

E^EMA^LAMliLXAMilii: 


HHHb 

""“"7 

■“Tf - 

mmm 

KMX® 

KNt«E 

KNOWilKl'tOi:^. 

■  ijj  (0 

KTOOE 

(O 

KNOJ.*: 
(7)  (8) 

KNtXtl'jlCNOOG 
(9)  I  (10) 

! 

EOELD  B2RMAT  CQIffiEfcilS 

KNODE  5(215)  Nodtd  displacements  printed  for  nodes  beginning  with  KNODE(t) 

through  K.NODE(l-rd ),  inclusive.  Up  to  5  pairs  of  nodes  are 
permitted. 


LINS?.  3?  '  Beam  I-oati  and  Stress  Selection  (used  only  if  IOUl'(  1 0)  0) 


D£5.CRiErm: 


Element  numbers,  in  p)alr;>,  to  sprecify  whicli  stresses  <!re  to  be  pjrintcd. 
Maximuni  and  minimum  values  of  stress  ate  prrinted  at  botli  ends  of  selected 
beams. 


imMMlAmBXAMELE: 


cr 

1 

_  ■ 

4 

5 

(T  1 

12345678% 

mmm 

1 

I 

1 

i 

a 

KBEAM 

(L 

KBEAlvmBEAM 

HB 

iJmgwjjgalgSi 

laagsim 

HHi 

1 

1 

LJZ 

, 

_ 

J 

HELD.  EQRf^T  CQNTErfrS 

KBEAM  5(215)  Loads  and  stresses  printed  for  beam  elements  beginning  with 

KBEAM(l)  through  KBEAM(I+1),  inclusive.  Up  to  5  pairs  of 
elements  arc  permitted. 


Scat  Stnjcturc  Output  'Tnue  Intcival 


£Hi)S:Q?iEn.QiV;  Interval  iit  which  luxic  and  element  data  indic:tts<1  on  i.ities  30  and  3 !  are  to 
Ix'!  printed. 


r’'" . 

-—■j- . -“1 

4  5 

j 

1 

1 

i 

Lj 

15545071^)0'  !  534l67l‘5() 

1 2345678901 1 234.567890! 

j  DTSEAT' 

r  01515 

FaiFli)  FORMAT 

DTSEAT  FI  0.0 


CONTENTS 

Time  interval  in  seconds. 


LINE  33:  Material  Tyise  Number 


DESCRII^ION:  Material  type  designation  number.  Repeat  group  33  through  3.5  in  sa]uence 

NUMAT  times,  as  specified  on  Line  27,  one  sequence  for  each  material. 


FORMAT  AND  EXAMPLES: 


i 

1 

2 

3 

4 

5 

6 

7 

11234567890 

TI34567M) 

1234567890 

1234567890 

T234567890 

1234567890 

1254.567^1 

IMTYP 

MAT 

□ 

2024-T4  AL 

FIELD  FORMAT  CONTENTS 


MYTP 

15 

MAT 

AlO 

Material  type  designation  number.  The  element  data  on  Line  39 
specifies  die  material  type  by  referring  ro  this  nnmbtjr. 

Material  type  description  used  as  heading  for  material  property 
output. 


IJNE  34:  Material  Pioperties 


DESCRIPTION :  Material  physical  properties  as  described  in  Figure  A- 1 2. 

EmMAI.^ND_£MM£LE;* 


r  .  ■  i 

2|  3 

4 

6 

7| 

wsmmuwmM 

1 

1 

tmsmm 

12345^785^ 

E(l) 

E(2)  T  W) 

E(4) 

E(5) 

irfrnia 

E(7)  1 

W5m~  OTiii 

FIELI> 

K)RMAT 

CObriENTS  FOR  BEAM  ELEMENT 

E(l) 

FI  0.0 

Etensity  (Ib-sec^/in.^). 

E(2) 

FI  0.0 

M(xluius  ol  elasticity  ObAn.^). 

E(3) 

FI  0.0 

First  yield  stress;  Syi  Ob/in.2)  =  q  jf  elastic. 

E(4) 

FIO.O 

First  plastic  mcxlulus  (Ib/in.^)  =  0  if  elastic. 

E(5) 

FU).0 

Not  used. 

H(6) 

FlO.f) 

Ultimate  stress;  (Ib/in.^)  =  0  if  elastic. 

f:(7) 

F'lO.O 

Poissoii’s  ratio. 

UinBBKMRVISnVWtfWItMMWIvninHOi^  : 


LSNEi  35:  Material  !l-*rx>j:!'eri!KS  (ctxnsinued) 

PESCRIFFKJN :  Material  piiysical  ptoperties  as  described  in  I  ‘igurc  A  - 1 2. 


"T 

"  — 1 

6  7 

TIT^i367l9()iim567F;{l 

E(ir)“n 

E(^) 

EOS)  1 

Etll) 

E(12) 

r  "tsot: 

(515 

L 

FIELD 

E(8) 

E(9) 

E(l()) 

E(ll) 

E(12) 


POI^AT  CONTENTS  FOR  BEAM  ELEMENT 

FI 0,0  Second  yield  stress,  Sy2  (lb/in.2). 

F'10.0  Second  plastic  modulus  (lb/in.2). 

FI 0.0  Strain-rate  coefficient  =  0,  no  strain-rate  effect  considered. 

F 1 0.0  Strain-rate  exponent  ■-  0,  no  strain-rate  effect  considered. 

F 1 0.0  Explicit  moment  curvature  flag  1 ,  use  explicit  moment  curvature 

option  (plate)  0,  ignored  explicit  moment  curt'ature  option  (plate). 


v.impk-  for  ifu"  first  of  iwo  t/ronns  ilcicrniincd  bv  N*  )M  A  T  .!  ois  !  .h; 


A  6H 


LINE  36:  Beam  Cross-Secl ion  Djila 


DESCPJFnON :  Beam  ciement  cross-sectional  protK.rtics  as  described  in  Eignm  A  !  3. 

Repeat  gronp  36  and  37  NSF.CT  times,  as  sixecified  on  l  ine  71,  one 
sequence  lor  each  cross  section, 


ifmiipm 

r-  5 

4 

'  “  -  5 

fij  7 

■n3456’7^9’^ 

ra45678‘XtrOTU^7'?l9(') 

[ns  EG 

KLOS 

ABM 

FIXX 

nvY 

FIZZ 

mm 

BHi 

O.'^AV 

0.3028 

0.1514 

■KmE! 

FdHLD 

FORMAT 

CONTENl'S 

NSEG 

15 

Number  of  plate  segments  in  beam  cross  section. 

KIDS 

15 

Flag  for  closed-wall  sections 

KLOS-0;  closed  wall 

KLOS  =  I;  open  wall. 

ABM 

FI  0.0 

Cross-section  area  (in.^), 

FIX.X 

HYY 

FIZZ 

3F10.0 

Cross-section  moments  of  inertia  about  x,  y,  and  z.  principal  ax 
respectively  (111.“*).  The  cross  section  for  each  bea  n  element  is 
oriented  by  specification  of  a  pointer  node  on  the  y  -axis  in  the 
element  day  on  Line  39. 

'  (■  x.iinpk-  i( »i  iiU' 


two  pyonps.  ik'u  riiiinci!  (ly  NSlyi  I  on  i 


■\ 


a)  CraSvS-secticn  geometry 


Cross-Section 


ISESCi'RIFFlON'  Beam  eicmcnl  cross  scctiortil  dimensions  as  dcscrilxHl  in  Ingiirc  A- 1  ii. 

( I ^  Kc{>c:at  I., ire  37  NSEG  a-  Kf J3S  Utnes,  ftiilowing  Line  36, 

(2)  Repeat  the  setnicnce  of  Lines  Mt  and  37  MSfiCT  times,  as  defined  on  Line  27. 


r" 


pi  234.6?;7S9{y 


YC() 


Zd) 


.6 

I""" . ""“1^ 

.  _ 7 

iT3'4567fi90 

1  L234567}f9()!  123456780{)T1234567;390!  1234.^678901 

'Ti)"' 


'0.0 

- (n?54 

1' . oM;^! 

-0.590 

().590 

r  T).083 

■0P8^ 

“ini 

o.fiyipi 

-Ti.Wi 

1  oTiFOt 

0,0 

■Mi4 

110^3 

”1T.T9^ 

'(}Mi 

T5H34 

ITo 

r  ■  (r()g3 

FIELD 

FORMAT 

ii:d:iNTEEi’s 

Yd) 

FI  0,0 

Cios^-section  coordinates  of  point  a*  be^nnning  of  segment  I  (in.) 

Zd) 

FI  0.0 

(See  Figure  A-  I3a) 

T(I) 

FI  0.0 

Segment  diickness  tor  segmeiU.  between  poitus  I  and  T  +  \  (in.). 

* !  ;\ain[3c  tor  itic  ion  cross  sirction  hasod  on  jsSSrf  •  K  and  Ki  ( )S  ('  on  1  inn  hi 


LINE  38:  Ncxlal  Point  Data 


DESCRIF'nON:  Finite  element  ntxle  number  and  nodal  c(X)rdi nates  in  global  system. 

NOTCS:  Repeat  Line  38  NUMNP  +  NCORD  times. 

IimiAT-ANI3L£XAMPJLE:* 


“I 

3 

“7 

5 

6  7i 

12345678^0 

123456780(7 

T2545(j789S 

12342i67M) 

153456789(1  1234567890 

N 

~xC(Nr 

Y^CNT” 

ZC(N) 

‘mm 

8.0 

Oil 

FIELD 

TORMAT 

CONTENTS 

N 

15 

Node  number. 

5X 

XCCN) 

FI  0.0 

X 

YC(Nj 

FIO.O 

Y  coordinates  of  node  point  (in.). 

ZC(N) 

FIO.O 

Z 

i|i!c  lor  the  first  .'2  lines  !>.:sei.i  on  NLMNi’ 


21)  and  NLORi)  .  2  on  1  iiu-  27. 


LINE  39:  Element  Data 


DESCRIFHON :  Individual  element  property  descriptions. 

NOTE:  Repeat  Line  39  NUMEL  times. 

FORMAT  AND  EXAMPLE:* 


r— 1 

r  ^ 

r  4 

.5 

6 

1 

n 

Si 

T234367F5(T 

1234567890 

mmm 

M 

NCX5E 

(1) 

NODE 
(2)  , 

NODE 

(3) 

hfODE 

(4) 

NCXDE 
.  (5) 

■ssfiia 

NCX* 

(2) 

NODE 

NODE 

.  (2) 

liEn 

NODE 

(Ij) 

N-TO,  NOI3E 
(12)  (13) 

5 

2 

wm 

■IB] 

mm 

OlC 

000  oio 

FIELD  K)RM>^T 


M 

*  a  A 

15 

Element  number. 

NODE(l) 

NODE(2) 

215 

End  nodes. 

NODE(3) 

NODE(4) 

215 

Not  used. 

NODE(.5) 

15 

Stiffness  flag 

NODE(5)  =  0:  Use  plastic  beam  stiffness 

NODE(5)  =  1:  Use  elastic  beam  stiffness. 

N()DE(6) 

15 

Cross  section  type.  Tlie  first  set  of  Lines  36  and  37  is  assumed  to 
be  cross-section  type  no.  1 ;  the  second,  no.  2,  etc. 

NODE(7) 

15 

Pointer  node  for  orientation  of  initial  principal  beam  axis  y. 

NODE(8) 

15 

NODE(8)  -  2  for  beam  element. 

NODE(9) 

15 

Material  type  (assumes  1  if  left  blank  ). 

NODEflO) 

15 

Beam-end  conditions  (forces),  at  end  i.  Figure  A- 13(b) 

ABC  (packed  word,  right  justified,! 

A  -  Forc:c  rclca.se  in  .\  direction,  it  I 
H  [•orcc  release  in  y-diatetion,  if’  ) 
C  “  Force  relea.se  in  7.  diaa  tion,  if  1. 


N()I)H(1  1)  15  Beam  end  conditions  (mo:nerus).  at  end  i.  Figure  A  l.itb) 

Did'  (packed  word,  right  justified) 

D  -  Moment  a'iease  in  x dia'ction,  if  i 
F  Moment  lelease  in  y  iliaclior.,  it  I 
F  Moment  release  in  /.  direction,  d  1. 


'  i  saiiiple  toi  [lie  1  It  ill  oi  .’7  lines,  based  tui  N(  A1F1 .  7  on  I  me  .' 


N0DE(12)  15 


NODE(13)  15 


Bearn-end  conditions  (forces),  at  end  j,  I-igiire  A13(b) 
OIX?  (jjackeal  word,  right  justified) 

O  -  Force  release  in  x-direction,  if  1 
P  =  Force  release  in  y-direction,  if  1 
Q  =  Force  release  in  z-direction,  if  1 . 

Beam-end  conditions  (moments),  at  end  j.  Figure  A- 13(b) 
RST  (packed  word,  right  justified) 

R  =  Moment  release  in  x-direction,  if  1 
S  =  Moment  release  in  y-direction,  if  1 
T  =  Moment  release  in  z-direction,  if  1. 


A  ?■: 


UNE  40:  Seat  Pan  N«les 


Ntxies  on  which  seat  cushion  loads  will  be  applied,  and  which  are  used  to 
define  the  seat  pan  outline  for  the  (x;ciiparn  plots. 


EORMAIAND-EMMi^Ui: 


[ 

1 

mmmm 

mttm 

■PUP' 

- ^ 

”■"”7] 

mmm 

mmam 

mmm 

UiSMP 

NP,AN  1  NP.'\N 
(3)  (4) 

NPAN 
(<7  j 

NPAN  1  NPAN 
CO  (E) 

NPAN  j  NPAN 

NPAN 

NPAN 

on.. 

13|  14 

[ _ J 

'7 

1 

14^  15 

_ Z. _ «l _ [5 

_ 16 

_ 

BOem  tmtEMS 

NPAN  12If)  Nodes  on  which  seat  cushicn  loads  are  to  be  applied,  input  on  retir 

edge  first,  then  forward  edge,  and  from  right  to  left,  as  shown  in 
Figure  A- 14. 


Note  that  in  this  example  node  6  is  used  as  NPAN(2)  and  NPAN(5),  ncxle  14  as  NPAN(4)  and 
NPAN(7),  node  7  as  NPAN(6)  and  NPAN(9),  and  node  15  as  NPAN(8)  and  NPAN(1 1). 


LINE  41:  Seat  Back  Nwles 


DESCRIPnON  Nodes  on  wliich  back  cushion  loads  arc  to  be  applied,  and  which  are  used 
to  define  the  seat  back  outline  for  the  mxupanr  plots. 


i 

1 

r 

L 

. 4 

. ■■■■'5 

11234567890 

mmsm 

mmm 

mmm 

tlMMWfBgiM 

N.8AK  1 

NBAK 

NBAK 

r®AK  j  hba:< 

mAK 

NBAtt 

NBAK  j  NItAK 

[NBAK 

f'iBAK  1 

JH 

Ljzl- 

LjiL. 

Lj5j4_16i„ 

-JiLJ 

. (^0  J  (10) 

^11), 

(12) 

5 

_ 

iOj 

ill _ 1 

_ .ic 

Lrt^ 

19 

L..^ _ 

mm  Emm: 

NBAK  1215  Ncles  on  which  back  cushion  leads  are  to  be  applied,  input  on 

lower  edge  first,  iht  n  top,  and  fa:>m  right  to  left,  as  shown  m  Figure 

,Vi4. 

Note  tliat  in  this  example  .lode  6  is  used  as  NBAKClri  and  ?^EAK(5),  node  7  as  NEAK(6j  and 
NBAK(9),  etc. 


mSE  42-  Restraint  System  Anchor  Point  Nocles  (NOCC  lines) 

DESCRIFTTQN:  Nodal  points  on  seat  structure  to  which  restraint  system  is  attached  as 

shown  in  Figure  A- 14. 

EQEMATMDlMMJPI£: 


r  '  f 

“3  ■  '  ^ 

6 

_ z 

mmm 

wmBmfmmm 

"l'S4567^ 

123456789(1 

NIMA 

:H| 

NSHA 
-  0) 

NSi-TA 

1.  g) 

NTO 

10| 

j 

"TC 

j.i 

__  12 

1 

HEIJJ 

TORMAT 

CONTENTS 

NLBA 

215 

Seat  structure  nodes  at  which  lap  belt  is  attachvjd,  right  side  first, 
then  left,  as  shown  in  Figurc  .A- 14.  Not  usext  if  lap  belt  is  attached 
to  aiicraft  floor  ratlier  than  to  the  seat. 

NSFIA 

215 

Seat  stnicture  nodes  at  which  shoulder  harness  load  is  to  be  applied 
(one  node),  or  distributexl  (two  n«xles).  l-eave  blank  if  shoulder 
harness  is  not  used  or  not  attaclied  to  seat. 

NTD 

15 

Seat  structure  node  at  which  lap  belt  tiedow'n  strap  load  is  to  be 
applied.  Ix;ave  blank  if  tiedown  .stiap  m  nor.  us<;d. 

Note  that  for  this  exantple  node  10  is  used  a.s  both  NrLBA(2)  for  passengct  1  and  NliBA(  l)  for 
passenger  2,  and  node  1 1  is  used  as  both  NLBA(2)  for  passenger  2  and  NLEA(i)  for  passenger  3, 
as  common  points  of  attachment  for  lap-belts  are  usually  found  on  transpoit  seats,  Also,  in  this 
example  no  s'loulder  harness  is  used. 


LINE  43:  Node  Constraint  Data 


DEvSCRIFTTON:  Packed  (encoded)  word  for  each  nodal  point  that  is  constrained  in  at  least 

ont5  degree  of  freedom. 


NOTE:  ( 1 )  Repeat  Line  43  NTJNIDS  times.  Omit  if  NUMItS  =  0. 

(2)  If  any  of  the  displacement/rotation  codes  is  set  to  2  to  enforce  floor  warping, 
include  Line  44  immediately  following  the  conesponding  Line  43  data. 


FORMAT  AND  EXAMPLE:* 


1 

2| 

4 

5 

6 

7 

1234567890 

1234567890 

1234567890 

1234567890 

1234567890 

1234567890 

1234567890 

NODDIS 

■um 

FIELD  FORMAT  CONTENTS 


NODDIS  1 1 0  Packed  word  -  N  ABCDF.F  (right  j  ustified) 

N  :=  Node  number 

A  -  Displacement  code  in  X  direction 

B  -  Displacement  code  in  Y  direction 

C  =  Displacement  code  in  Z  direction 

D  =  Rotation  code  in  X  direction 

F,  -  Rotation  cotle  in  Y  direction 

!'  -  Rotation  ctxle  in  Z  direction 

A.  B,  C,  D,  E,  or  F  =  0,  no  constraint 

=  1,  constrained  for  zero 

displacement/rotation 
=  2,  constrained  for  n(H)r  warp 

displacenient/rotation. 


ii\e  MIK'S, 
li'  '  ;  iMe  It 


1 1  n  \  i  MI'S 


\ 


LIN  E  44:  Floor  Warp  Data 


DESCRIPTION:  Floor  warp  displaccment/rotation. 


NOTE:  (1)  Repeat  Line  44  lor  each  ci'splacement/rotation  ccxle  set  to  2  in  Line  43,  in  the 

order  from  displacement  in  X-direction  to  rotation  in  Z-direction. 

(2)  Rotations  are  input  in  radians;  displacements,  in  inches. 


FORMAT  AND  EXAMPLE:* 


1 

- 1 

“I  4 

5 

6 

7 

1 12^56789(1 

1234567890 

ITB4'567890|T234W890| 

12345678‘^d 

12345678^0 

1234567890 

FWARP 

-0.5 

FIELD  FORMAT  CONTENTS 


FWARP 


FI  0.0 


l-kxrr  warp  displacement/iotation. 


If  ccniact  with  the  seat  back  is  to  be  simulated  (by  IOllT(4)  ~  1  or  2),  the  following  li?ies  of  input 
data  are  required  to  describe  the  surfaces  on  die  seat  brick  (only  witli  NSEAT  =  0  on  Line  3). 
ITtese  lines  would  direedy  follow  Line  If-  (with  the  finite  element  seat  model  data  omitted). 


LIME  45:  Seat  Back  Contact  Surface  i)imcnsioi;K 
DESCRII^'tOM:  Dimensions  of  contact  surfaces  on  sec. 


I  -Mck,  as  illustrated  in  Figure  A- 15, 


FORMAT  AND  EXAMPLE: 


1 

SSHHP 

- 

6 

7 

!234567.?90 

123456789()|  1234567890 

1234567890 

1234567890 

TTI' 

\,VTT 

HTT  1  TAR 

WAR 

HAR 

XI.AR 

s 

_ 

_ 

_ 

FIELD 

FORMAT 

CONTE N IN 

ITL 

FT  0.0 

Distance  from  top  of  seat  back  to  top  edge  of  stowed  tray  tabic  (in.). 

WIT 

FT  0.0 

Width  of  tray  tal  (in.). 

HTT 

^TO  0 

Height  of  tray  table  fin,). 

TAR 

FT  0.0 

Distance  from  top  of  seal  back  to  top  of  annrest  (in.). 

WAR 

FT  0.0 

Width  of  armre.si  (tn.). 

HAR 

F10.''i 

Height  of  arrniest  (in.). 

XI.,AR 

FT  0.0 

Length  of  anttrest  (in.). 

LINE  46-48;  Seat  Back  Force-Defonnati(>p  Prop<Mlies 


DESCP-IEnC/N :  Force-deflection  characteristic:;  and  daioping  for  st'at  back  suifaces,  'I'he 

force  is  computed  from  deflection,  6,  accttrdiiig  to  I-  =  ■  I ). 


FQFfMAT  AND  EX/  MPLR: 


r" - 1 

"  -  3 

4  5| .  “  bT""  "  “T 

1234567890 

1234.567890 

1234567890 

1234567890  1234567890!  1234567890!  i2;^4567890 

CCON(l) 

BCON(l) 

DCONd) 

CC6N(2) 

BCON(2) 

DCON(2) 

(:CON(3) 

BCO.Nf3) 

DCON(3) 

'  i 

FIELD 

FORMA!' 

CONTENT'S 

(ICON 

FI  0.0 

Coefficient  C  in  above  equation  (lb). 

BCON 

FI  0.0 

Coefficient  B  (in."*). 

DCXDN 

FI  0.0 

Damping  ctxjfficient  at  zero  load  (!b-sec/in.) 

Three  lines  of  data  are  input  in  the  above  format.  The  first  (46)  apjtlics  to  the  cushion  surfaces  (I, 
3,  and  4,  in  Figure  A-J  5).  The  secottd  line  (47)  applies  to  the  tray  table  (2  in  Figure  .A-l.o).  The 
third  refers  to  the  imnnesi  surfaces  (5-8  ut  Figure  A- 1 5). 


LINE  49:  Seat  Back  Weight  and  Row  Pitch 


DESC’RIPnON :  Weight  of  the  movable  seal  back  for  use  in  the  breakover  mode!,  damping 

coefficient  for  seat  back  breakover,  and  seat  row  pitch. 


BQEM4T.Al!fllEXAMELE: 


T 

4 

5 

- —B 

7 

r234567890 

1234567890 

hQ34567890 

1234567890 

T534537890 

BKWT 

dpbo 

SpTtTm 

HELD 

FORMAT 

CONTENTS 

BKWT 

FIO.O 

Seat  back  weight  (lb). 

DPBO 

FIO.O 

Damping  coefficient  (in.-lb-sec). 

SPITCH 

FIO.O 

Seat  row  pitch  (in.). 

A 


LINE  50:  Seat  Back  Breakover  Resistance 


DESCRIPTION:  Seat  back  breakover  moment  versus  rotatioin  angle,  similar  to  that  shown  in 

Figure  A -10. 


FORMAT  AND  EXAMPLE: 


1 


11345671}^ 


ffb(5(:>) 


TU4567§9?) 


ffbOgT 


3! 


1234567890 


1234567890 


5 


T254567M 


1234567890 


1234567890 


“FFB0(4) 


DDBO(2) 


DDB0(3r 


DDBO(4) 


FIELD 

FORMAT 

CON-’IENTS 

F’FJ30(2) 

FFB0(3) 

FFBO(4) 

3F10.0 

Resisting  moment  of  seat  back  (in. -lb). 

DDBO(2) 

DDBO(3) 

DDBO(4) 

3F10.0 

Angular  displacem.ent  of  seat  back  (rad.). 

A  S(. 


APPENDIX  B 


EXAMPLES  OF  OCCLTPANl’  CHARACTERISTICS  AND  MATERIAL  PROPERTIES 

A  significant  problem  encountered  in  mathematically  modeling  a  physical  system  lies  in 
determination  of  system  characteristics  and  properties.  In  tliis  appendix  are  presented  examples  of 
the  following: 

•  Occupant  dimensions  and  characteristics. 

•  Restraint  system  webbing  load-elongation  characteristics. 

•  Qishion  load-deflection  characteristics. 

•  Stmctural  material  stress-strain  characteristics. 


The  characteristics  and  properties  contained  in  this  appendix  are,  of  course,  not  intended  to  be  all 
inclusive,  but  rather  are  intended  to  provide  the  program  user  with  examples  that  may  aid  in  setting 
up  new  input  cases. 

B.l  OCCUPANT  MODELING  CHARACTI-.RTSTICS 


As  described  in  Chapter  2,  dimensions  and  inertial  properties  for  two  standard  occupants,  a  50th- 
percentile  civilian  male  and  a  50th-percentile  anthropomorphic  (Part  572)  dummy,  are  included 
within  the  program.  If  a  nonstandard  occupant  is  desired,  additional  data  must  be  provided  on 
Lines  22A  tlirough  22L.  The  fonnat  for  nonstandard  occupant  data  is  displayed  in  Figure  B-1,  and 
parameters  are  defined  on  pages  A-34  through  A-49.  In  Figure  B-2  are  presented  the  properties 
that  ai'e  used  in  the  program  for  the  standard  (Part  572  50th-percentile)  dummy  occupant.  Figure 
B-3  presents  a  set  of  data  for  a  95th-percentile  dummy,  w'hich  have  simply  been  scaled  from  the 
50th-percentile  data.  The  use  of  this  scaling  method  is  not  suggested  if  measured  properties  can  be 
obtained;  however,  to  complete  a  partial  set  of  properties  or  obtain  a  quick  estimate  of  the  solution, 
use  of  the  scaling  approach  can  be  justified. 

The  scaling  method  is  based  on  multiplying  the  50th -percentile  properties  by  the  appropriate 
nondimensional  scaling  factor.  All  properties  with  length  dimensions  are  multiplied  by  the  ratio  of 
nonstandard  occupant  sitting  height  to  .50th- percentile  sitting  height.  In  this  example: 


Irength  Factor 


95th  %  Sitting  Height 
5()th  %  Sifting  Height 


37.8  in  ^ 
35.7  in 


1.06 


Similarly,  (Kxupant  properties  based  on  weight  are  scaled  by  the  (xcupant  weight  ratio,  i.e.; 


Weight  Factor  = 


95th  %  We^ht  212  lb 
5()th  %  Weight  IM  lb 


1.29 


The  factor  for  .scaling  moments  of  inertia  was  deriveri  from  a  dimensional  analysis  for  the  vtu'iables 
involved.  Tie  resulting  scaling  factor  iv; 


,,,  .  (95lli%Weight)2(95!h%SiltingHeight)- 

Moment  of  Inertia  l  acior  ^  '  c  /  ' .  -  / 

(5(hh  %  Weight)^  (.5()th  %  Sitting  I  leighi)’ 


1,1,5 


it  I 


vSince  there  is  no  valid  basis  for  scaling  stiffnesses,  the  50th-percentiie  spine  and  neck  stiffness 
properties  were  retained. 


B.2  WF.BBTNG  I.OAD-ELONGATIQN  CHARACTERISTICS 

Figures  B-4  and  B-5  present  static  load-elongation  characteristics  for  several  types  of  nylon  and 
polyester  restraint  system  webbing,  respectively.  Very  little  dynamic  data  for  webbing 
deformation  exist;  however.  Figures  B-6  and  B-7  present  some  dynamic  results  taken  from 
reference  B.l. 

The  damping  coefficients  for  the  restraint  components  are  based  on  three  assumptions:  that  the 
webbing  damping  coefficient  is  not  a  function  of  strain  condition,  that  it  is  independent  of  strain 
rate,  and  that  the  Voigt-Kelvin  ntodel  (shown  in  Figure  B-8)  can  be  used  to  represent  the  webbing. 

Tile  first  assumption  allows  the  use  of  a  linear  approximation  to  the  static  and  dynamic  load-sPain 
curves  for  the  webbing  material.  The  single  slope  approximation  should  be  the  best  estimate  for 
the  expected  range  of  webbing  loads,  and  not  for  the  entire  curve.  The  second  assumption 
indicates  that  the  damping  coefficient  will  be  applicable  to  all  possible  strain  rates  encountered  in 
the  simulation.  The  accuracy  of  the  damping  coefficient  can  be  maximized  by  basing  tJie  calculated 
value  on  dynamic  webbing  test  data  measured  at  an  applicable  strain  rate.  The  procedure  for 
calculating  the  damping  coefficient  for  nylon  webbing  (MIL-W-4088  TYPE  VII)  is  given  below. 

The  static  load-elongation  curve  for  the  nylon  webbing  sample  is  shown  in  Figure  B-9.  A.  linear 
approximation  tc  this  curve  is  1 1,000  lb/in./in.  over  the  expected  load  range  of  0  to  2000  lb.  The 
slope  of  the  dynamic  test  data,  measured  at  a  strain  rate  of  40.9  in./in./sec,  is  approximated  as 
26,000  lb/in./in.  Based  on  the  assumption  of  a  parallel  spring-damper  model,  the  dynamic  load  at 
any  elongation  value  must  be  equ.il  to  the  static  load  plus  the  damper  force,  i.e.. 


PoYNAMIC  =  PsTATTC  +  PdAMPING 
=  Ke  O 


(B-1) 


Where; 


K  is  sloj>e  of  the  .-•''ad-strain  curve  (Ib/in./in.) 
C  is  the  damping  coefficient  (Ib-secfin./in.) 

e  is  the  strain  (in./in.) 

fc  is  the  stnain  rate  (in./in./sec) 


Thtrefore,  the  dan:.ping  coefficient  can  be  calculated  using 


Pffi'NAMIC:  Kf  (B  2) 

Using  as  a  representative  point  a  liyitatiiic  load  -ol  2000  ih  and  0.0K2S  m  /in.  strain,  the  damping 
coefficient  for  the  nylon  weblung  is  caicuiaied  its 


H  2 


c  -  (n/XX)  lh/in./in)(().0825  in. /in.) 

40.9  in./in./sec 

™  26.7 

in./xn. 


):DEELEaaOIi-.CmRMri^^ 


The  seat  cushion  represented  in  Program  SOM-LA/SOM  TA  accounts  for  the  stiffness  and 
damping  properties  of  the  cushion  combined  witli  tlie  occupant  butttxrks.  Tlxis  naxielirig  appixxach 
is  desirable  in  order  to  avoid  the  numerical  problems  associated  with  springs  ixi  a  seiies 
configuration.  An  experiment  was  performed  to  develop  load-deflection  properties  for 
representative  cushions.  The  experiment  consisted  of  applying  a  known  static  load  in  the 
downward  direction  to  the  lower  torso  segment  of  an  Alderson  VIP-95  dummy.  Tins  downward 
load,  Vi'hich  was  applied  at  the  spine  base  plate,  caused  both  the  buttock.;  and  cushion  to  deform. 
The  deflections  of  the  combined  system,  buttocks  and  cushion,  and  the  buttocks  separately  were 
measured  for  each  applied  load. 


A  description  of  the  cushions  used  in  load-deflection  tests  is  given  in  Table  B- 1.  'Hie  cushions 
were  selected  to  provide  a  spectnim  of  tlie  possible  cushion  configurations  that  the  user  may  select. 
Combined  load-deflection  curves  for  the  VIP-95  buttocks  and  cushions  are  presented  in  Figurcs  B- 
10  through  B-14.  The  fonri  that  the  load-deflection  cums  tiike  is  a  lineiw  slope  followed  by  an 
exponential  stiffening  as  the  cushion  and  occupant  "bottom  out."  These  curves  can  be 
approximated  by  an  expression  of  the  form: 

F  =  C(e»8-l)  {B3) 

Representing  the  load-deflection  curves  with  a  smooth  function  alleviates  a  convergence  problem 
encountered  previously  with  the  numerical  integrauon  around  the  slope-change  points  of  a 
piecewise,  linear  representation.  The  exponvcntial  representation  of  the  five  load-deflection  curves, 
developed  with  a  leasSt-squares  approximation  .routine,  is  presented  as  the  dashed  line  in  each 
figure.  Also  presented  in  this  section  are  the  sepai’ate  ioad-deflection  curves  (Figure  B- 15)  for  the. 
Alderson  VlP-95  dummy  buttocks  when  tested  with  each  of  the  five  cushion  types.  This  i.s 
piesented  for  the  user  wdio  may  want  to  .s}'nthesize  a  combined  load  deflection  curve  by  adding  the 
<iesiR‘.d  cushion  properties  determined  undex-  a  rigixl  indenter  to  an  average  deflection  curve  for  tlie 
dummy  buttocks.  'ITie  indenter  .should  be  confignied  like  the  dummy. 

B.4  SIRUCrURAL  MAIFRIAL.STRBSti- STRAIN  CURVES 

Figures  B  16  through  B-20  present  approximated  stress -sira in  cun'es  for  three  steels  and  two 
aluminum  alloys.  From  each  of  these  cntvss,  six  characteristics  are  pi-ovided  as  input  to  the  finite 
element  seat 


T/vB.[.,E  B-1.  1'Yl‘E  DHSlGNATtDN  AND  DI'SCSIFriON  OF'  CUSHIONS  {•OR 
_ _ _ _ _ _ _ _ _ _ _ 


Typz 

Nauiffi£i„  _ _ _ I><;i!.aaj£U.aa _ 


1  Coritoui'Cil,  miiit'’!ay(;;.red  cusiucn  designed  to  minimize  occupant  rebound  in  a  crash 
situation. 

2  Contoured,  rigid  tbanr  cushion  designed  for  negligible  deflection. 

21  Contoured  funiiture  foam  cushion  approximately  1.5  in.  thick  (undefonned)  over 

buttock  contact  .«ea. 

4  Furniture  foam  slab,  1.2  ib/ft-^  density,  approximately  3.0  in.  thick  (undefonned). 

a  Furniture  foam  slab,  1.4  Ib/ft^  density,  appioximately  3.0  in,  thick  (undefonned). 
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oad-elongation  characteristics  for  polyester  v/ebbing. 


K  (!b/in.in.) 

C  (ib-sec/in./in.) 


B-8.  Voigt-Kelvin  model  of  restrain  +  system  webbing, 


bOOO 


Exponential  representation:  F  =  73.6  S.-jj 
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lire  B-iO.  Combined  load-deflection  curve  and  expcnendal  representation 
for  Type  1  cushion  and  VIP-95  dummy  pelvis  and  buttocks. 
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Figure  B-i  1.  Combined  ioad-dcflection  curve  and  exponential  representation 
for  Type  2  cushion  and  VIP-95  dummy  pelvis  and  buttocks. 


Exponential  representation 


Combined  load-deflection  curve  and  exponential  representation 
for  Type  3  cushion  and  VIP-95  dummv  pelvis  and  buttocks. 


Exponential  representation 


Combined  load-deflection  curve  and  exponential  representation 
for  Type  4  cushion  and  VIP-95  dummy  pelvis  and  buttocks. 


Exponential  representation: 
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i  B-14.  Combined  ioad-defiection  curve  and  exponennal  represeritation 
for  Type  b  cushion  and  ViP-95  dummy  pelvis  and  buttocks. 


Lead '  lb. 


Deflectilon  -  in. 


F^'ure  B- 15.  Load-defkxtion  curves  for  AkiersoD 
Vil’  95  dummy  pelvis  and  bint<K‘k;i 
tested  witli  five  differeiit  cushions, 
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ApproKsmiai"-  Material  Properties 

Modulus  of  elasticity,  E(2)  30  x  U/^psi 

First  yield  stress,  H(3)  ==  58,700  psi 
First  plastic  modulus,  E(4)  2.9  x  lO^’psi 

Ultimate  stress,  E(6)  ”  67,000  psi 
Second  yield  stress,  E('8)  =•■  63,500  psi 
Second  plastic  modulu-s,  E{9)  --  75,000  psi 
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Figure  B -  i  6.  t’iecewisc,  liueai  approxiniaOon  o!  stress- 
sirain  curve  for  ;0i0  coFt  drawn  siccl. 
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Approximate  Material  Properties 
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Stress  -  ksl 


Approx^ate  Material  Properties 

Modulus  of  elasticity,  E(2)  =  29.1  x  10  ^  psi 
First  yield  stress,  E(3)  =  160,000  psi 
First  plastic  modulus,  E(4)  =  7.8  x  10  psi 
Ultimate  stress,  E(6)  =  180,000  psi 
Second  yield  stress,  E(8)  =  170,000  psi 
Second  plastic  modulus,  E(9)  =-  8.85  x  10  ^  psi 


Strain  •  percant 


Figure  0-18.  Typicai  len.s'le  strcs.s-<ilnun  curve  for  AI.S!  4340  steel, 
heat  treated  to  180  ksi  ullitnaie  stress  and  pi'-ce wise, 
iistear  aji|>n)xirnatiori  to  curve 


Stress  -  ksJ 


Approximate  Material  Properties 

Modulus  of  elasticity,  E(2)  =  10.5  x  10  psi 
First  yield  stress,  E(31  =  44,000  psi 
First  plastic  modulus,  £(4)  =  4.9  x  10  ^  psi 
Ultimate  stres.s,  £(6)  =  62,000  psi 
Secoiid  yield  stress,  E(8)  =  58,(X)0  psi 
Second  plastic  modulus,  E(9)  "  6.2  x  10  ^  psi 


Figure  B  19.  Typical  tensile  stress  sirain  curve  ior 

2024  T4  alum!  urn  aHoy  and  picv  c  ss  i  sc, 
lt,riear  ajrpro.vj  iiatiuii  to  curve 


Approximate  Material  Properties 

Modulus  of  elasticity,  E(2)  =  9.9  x  10  ^  psi 
First  yield  stress,  E(3)  =  36,200  p.si 
First  plastic  modulus,  E(4)  =  1.1  x  10  ^  psi 
Ultimate  stress,  E(6)  =  42,000  psi 
Second  yield  stress,  E(8)  =  40,200  psi 
Second  plastic  modulus,  E(9)  =  3.0  x  10  psi 
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APPENDIX  C 


PROGRAM  STRUCTURE 


Til*?  0’(-srall  orgartiAidon  of  Program  SOM-IA/SOM-TA  i?  illustf?.ted  in  Figuif^  C-1.  The  maitr 
progrem  controis  Lhs  overall  solution  piocedare  by  calling  tw  o  indivlaual  sets  of  snbrootines,  on?; 
for  the  occupant  stignisuts  of  the  program,  and  the  other  for  the  sear  segment  of  the  program. 
Detailed  descriptions  of  the  occupant  subroutines  are  presented  in  Section  C.i,  and  the  seat 
subroutines,  in  Section  C.2. 


At  the  start  of  execution,  the  main  program  calls  subroutine  INFF  to  read  input  data  for  the 
occupant  vnodel  and  subroutine  READIN  for  seat  input  data,  Subroutines  CONST  and  II-ilTIL 
calculate  constants  and  detemiine  initial  values  of  generalized  coordinates  for  the  occupant,  and 
subroutine  ASSBLE  performs  prelimina'^'  calculation.?  for  tlse  seat  model.  ITien,  a  solution  loc^p  is 
entered  at  initial  time  and  passed  through  for  each  time  step.  During  each  pass  through  the  solution 
loop,  subroutine  RKAM  advances  the  solution  for  the  occupant  equations  of  motion  one  time  step 
and  provides  forces  to  be  applied  to  the  seat  model  by  die  occupant.  After  the  call  to  RKAM,  if  the 
finite  element  seat  model  is  being  used,  subroutine  SOLVE  advances  the  solution  by  ihe  seat 
.stnictural  analysis  to  the  same  point  in  time  that  has  been  attained  by  RKAM.  At  lime  intervals 
aek'cted  by  user  input,  subroutine  ANSWER  stores,  in  arrays,  user-selected  items  of  output  data. 
Daia  for  po.  t-processing  plot  programs  are  written  on  external  files  14  and  20  for  die  occupant  and 
seat,  respectively.  Additional  data  are  written  on  unit  26,  as  de.scribed  iu  Section  3.4.  These  files 
must  be  saveri  if  plots  are  desired 


C.I 


The  relationship  among  the  subroutines  .in  the  occupant  segment  of  the  program  is  illustrated  in 
Figure  C-2.  Indiv  dual  subroutines  are  described  below. 

C.  i .  1  Subroutine  AMATRX.  Called  by  EQUATF';  calculates  elenietus  of  the  inertia  matrix  ( A] 
for  the  three -dhriensional  (x;cup;mt  model. 

C.1.2  Subroutine  AMATX2.  Called  by  EQUAT?,;  calculate.?  eleme.nts  of  the  inertia  matrix  for  the 
plane-motion  txrcupant  model. 

C.  1.3  Subrourinc  ANSWER.  Called  by  M.AJN;  calciihues  accelerations  Ar(I,  J)  of  body 
.segme.Tis  in  the  inerliai  ccvrrdinate  system  'i.nd  iTari.sf(>m!s  the  accelerations  to  segment  fixed 
c(K>rdinate  sy;stC'.n;s  Ca’rulates  severity  indicet'i  juid  organizes  position,  velocity,  and  foR-e  data  for 
output  Wrtes  plot  data  oo  units  bt  ana  26.  if  dalu  fihcring  is  atqiiested  by  user  iiquit,  ANSWER 
w,d:es  (K'cup.i'U  acedriaiums  on  unit  9  and  sea:  acceiera.tioris  on  unit  10  for  subsequent  filtering  by 
sabrnuune  OU TPT  Also  <-aIied  by  EQUATE,  HQCA'n,  !  OBODY,  E1NV3I-,  or  SOLVE  iiuhe 
event  of  abnormal  tenriina^K!!! 


c  E4  Casitd 

accrie ration  co'tifKsnenis  tu  thr  aircraft 
iKXC;vn'U.i-')t:  to  determic.'  ve 
meitiai  sysieiu;.  When  iuTM;  r;  greiUer 
i.iiai  veiocitv  thftn  Rroams  constain. 


umialiv'  .iNE'E,  then  by  lAlS'f'ON;  calculates  current 
IUku'  (A(  ’C(,1),  I  =  1, 3)  based  cm  input  iicceleration  pulses, 
iloviry  (in  aircrafi  cooniinatc  .system)  ami  displaccmeni  (in 
that!  the  ripui  jsuise  diuotiun,  aecelcialion  i.s  set  to  zero,  so 


(  1 


Figure  C-l.  Overall  organization  of  Program  SOM-TA. 


igure  C.2  SOM-LA^SOM-TA  Program  Stmcture;  Occupant  Segtnent 


i'‘.  1,5  Subroutine  A.UXSUB.  Called  by  RK,AM  {also  ittitially  by  MAJlNj;  calculates  clori natives 
atid  forms  two  1  k  2N  anays  of  ^^'^l■iables  aiiddeiivarive:;-. 
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where  N  is  the  nurnber  of  vlegrces  of  freedom,  either  12  for  ilic  plane-moti0.ri  model  or  29  for  the 
three-dimen, sional  model  The  velocity  and  acceleration  of  the  DR[  rnotlel  are  assigtiCd  to 
DER(2N-rl)  and  DER(2NTi),  res;pecdvely.  If  the  two- deiP'ee- of- freedom  energy-absorbing  seat 
model  is  used,  it,s  velochies  md  accelerationj'.  are  assigned  to  DE'R(2N-t-3';  through  DEiri(2N-r6). 

EQUATE  is  called  to  provide  va?i>es  of  the  derivatives  (the  gci.erali^ed  velocities,  QD(J),  and 
accelerations,  QDD(j)).  RKAM  ilten  integrates  the  two  sy.stem®  of  fust-order  equations. 

C.1.6  ouoroutine  BMATR^-  Celled  by  IKHIATE;  calculates  elements  of  velocity-dependent 
vector  {Bj  forthe  three  dimensional  nrtcsiel 

C.1.7  Subroutine  BrA'lATX2.  Called  by  EQUAT2;  calcr’ates  elements  of  vector  (B)  for  the 
plane-motion  model 

C.1.8  Subroutine  BUC.iILE.  Called  by  FORCES;  deiormint.s  (wsiaon  of  the  fKrint  of  intersection 
between  abdominal  centre?  surface  and  thigh  contact  surfaces  <  projected  on  X-Z  plane). 

C.i.9  Subroutine  CABIN.  Calle  J  by  INPT;  computes  coordinates  fo.f  the  image  of  the  seat  in 
front  of  that  being  modeled,  Tbe.se  coordinate.s  arc  fhen  available  for  use  by  the  passenger  plot 
program.  If  seat  back  breakover  is  being  modeled,  CABJK  c,dit(i  ly  Al'XoUB  to  conipaie 
derivatives  for  seat  back  motion. 


C.l.iO  CaUed  by  INPl,  once  Um  eacn  passenger;  h'-st-d  cn  input  data, 

‘  4  I  /•'  I  »  t  '  ‘  f  fa  O  1  '  '■»  I  »  l  t’  <  «  fr-  ^  r.\  <  •  J  «•  '»-j.  .vv-vr-.. 


alculate.s  values  of  parameters  that  remain  ctuistant  ihroiJgl'o.c  pitnuam  evecctitin.  'The:- 
onstants  include  functions  of  voccupant  diinenMOf's  use.d  in  the  equations  ..af  motion  and  jon 


constants 
reaustance  paran'eters 


se 
joiin 


CONST  also  w-rites  (xteupaiU  difiicfcsiops  on  mat  14  tur  piiHtiug, 

*  l  Id  1  ‘hibroutiii;,  EflUAlT  Ct-llea  bv  AUXSuB  for  tlu;  dn  v  dirnensmna!  occuparu,  mscs  ihc 
'atcsi  values  of  gesie<'.d)2.cd  cootdinaics  ana  veloc-Ues  sOcab'uKOc  sc’u.',  in  r'crupaut  cijuaco  is  ,4 
moiioii,  Solves  equatioi’s  cif  mot, on  Tsi  ac  cc'enituuis  QU!  )(il. 

C'llls  AMA  IKS',  HMA'ITX,  yfd.\’!  TX,  RMAd’HX,  am!  Q'AA'i  px  a;  '-ettnu'  iqi  ‘hr-  cqualu-ns  of 

uKiiiou  uudealiv  !.  Tr  picir  saiuuon 


C.i.i;2  Siib’.'outmc  ECHJ '4,T2.  Called  by  AUXSUB  for  the  piane-rnoiion  occapant;  u:ies  latest 
valiiC-s  of  generalis'.ed  ccK'n  iinates  and  velticilies  to  calculate  terms  in  (Kxuj/ant  equations  of  motion  . 
S:dves  equations  of  motion  for  i.ccelctaiions  QDD(J). 

Calls  AMATX2,  BMAIX?.,  VMATX2,  RMATX2,  end  QMATX2  in  setting  up  the  equations  of 
motion  and  calls  1.21^  V 3F  fot  their  solution. 


C.L13  SubfOfttinp  F'BHLTS.  Cabled  by  FORCES  to  coinpute  restraint  system  forces. 


C.C14  Subrr)utine  FCUSHiV.  Called  by  FORCES  to  compute  scai  ciishicn  forces.  The  cushion 
forces  include  frictional  contponetus  whose  directions  oppose  tlte  current  velocity  of  the  occupant 
with  respect  to  the  cushion. 

C.1.1.3  Subroutine,  I'OBQDV.  Called  by  FORCES  if  IRSYS  -  0  for  lap  belt  only;  checks  for 
head/leg  and  chest/leg  contact.  If  contact  occurs,  the  force  is  ca'culated  based  on  an.  exponential 
function.  The  components  of  the  force  are  then  added  to  the  F  array  already  computed  bv 
FORt2ES. 

C.1.16  Submutine  FORCES.  Called  by  QMATRX  or  QMA1X2;  calls  FBELTS,  FQJSHN,  and 
FCBCDY  to  compute  extemai  fotces,  calculates  forces  exerted  on  the  occupant  by  the  floor.  Sums 
forces  and  transforms  them  to  inertial  coordinate  systc  >  ^or  aiuations  of  motion. 

The  forces  arc  placed  in  an  array  (F(I,  J), !  -- 1,  1 1,  J  -  1,  3)  for  use  in  QMATRX  or  0MATX2. 

If  the  two  -degree-of-fieedom  energy-absorbing  seal  ntodel  is  used,  the  a-ansiational  and  rotational 
ao.telcradons  are  calculated. 


C.  Li?  Subroutine  IMFACT.  Called  by  FORCES;  computes  the  point  of  closest  proximity 
bf  i.'ween  each  contact  surface  on  the  occupant  and  the  seat  in  front.  DELIMP(N,J)  is  the 

pc,  netration  of  occupant  surface  N  into  surface  J  on  the  seat  back  in  front.  If  DEL1MP(N,J)  S  0  the 
irapact  velocity  V.ELIMP(N,J)  is  calculated.  Calculates  forces  due  a)  seat  back  contact. 

C.  1  18  Subroutine  INITIL.  Called  by  MAIN,  once  for  each  passenger;  calculates  initial  values  of 
tl  i:  generalized  coordinates  and  '/elocities  for  the  occupant  and  the  initial  defiections  of  the  seat. 
O'^fTTlL  first  uses  input  values  of  GAM(J)  to  determine  the  angular  position  of  the  lKx.ly  segments  1 
through  7.  Based  on  the  air.raft  orientation,  the  occupant's  weight  ss  applied  u>  the  seat  and 
restraint  system,  and  the  po.sition  of  the  lower  torso  seg.nient  (Xj,  Vj,  Z^)  iS  determined.  From 

the  X  and  Z  cixirdinates  of  segment  1  (computed  hcfc)  and  of  the  ocrupan’  s  heels  (from  INPT) 
t.hc  position  of  the  leg  set,menis  is  calculated.  Thioughout  *hesc  con'ipu!.a».io.us,  the  body  is 
a.sstime.d  to  be  symmetric  with  respect  to  the  aircraf;  (X-Z)  plane. 

If  the  sear  is  iii.iriaily  waqied  so  thai  the  seat  back  angle  is  changed,  the  values  of  GAM(l), 
G.4M(2),  and  GA.Mi  3)  are  adjusted  accordingly 

In  the  event  that  the  input  initial  conditions  imptise  umeasonable  leqnirements  on  occupant 
geometry,  a  diagnosiic  sne-ssage  is  provided  and  execution  is  stopiicd, 

C,1.19  Subrouiiiic  INF'',  f ’ailed  by  MA IN,  reads  (.HX'upam  input  data  Detailed  (iescri(nioiis  of 
input  are  presented  in  C’hapter  2  and  A[;pef)dix  A. 


i-20  SuhrxniOuc  ItriNTS.  ('ailed  by  KMATRX  or  KM 
uans/tion  region  ol  liie  jmiu  stopping  moinee;  ,. 


rX2;  tCs  a  t  iihie  curve  into  the 


C.L21  SubroutitiR  I.DFN.  Called  by  FBELTS  and  by  rORCES;  uses  linear  interpolation  in  a 
table  of  force  (Y)  versus  deflection  (X)  values.  A  description  of  the  parameters  in  the  calling 
sequence  follows: 

X  A  table  of  the  independent  variable,  xj,  such  that  xi+ 1  >  xj  (if  ICf  IK  ==  0). 

Y  llie  tabic  of  the  dependent  variable,  yj  =  yfxj)  (if  IQdK  -  0). 

N  The  number  of  entries  in  each  of  the  above  tables;  i  =  1,...N. 

XA  The  independent  variable,  x,  for  which  interpolation  is  requested. 

XLAST  Previous  values  of  XA  and  YA. 

YI^ST 


ICHK  Index  which  is  0  or  1  depending  on  call  during  loading  or  unloading. 

XCURV  Point  on  loading  curve  from  which  unloading  started. 

YCURV 

C  Unloading  slope,  used  only  if  lUNLD  =  2. 

YA  The  dependent  variable,  y  =  y(x),  being  determined. 

lUNLD  Index  which  is  2  if  unloading  slope,  C,  is  to  be  used,  1  if  unloading 
proceeds  along  basic  loading  function. 

C.1.22  Subroutine  LINV3F.  Called  by  EQUATE  and  EQUAT2;  performs  linear  equation 
solution; 


A 

B 

IJOB 


N 

lA 

Di,D2 


Input/output  rtiatrix  of  dimensions  N  x  N.  See  parameter  IJOB. 

Input/output  vector  of  length  N.  On  input,  B  contains  the  right-hand  side  of 
the  equation  AX  =  B.  On  output,  the  solution  X  replaces  B. 

Input  option  parameter.  IJOB  =  2  implies  solve  the-  equation  AX=B.  A  is 
replaced  by  the  LU  decomposition  of  a  rowwise  permutation  of  A,  where  IJ 
is  upper  triangular  and  L  is  lower  triangular  with  unit  diagonal.  ITie  unit 
diagonal  of  L.  is  not  stored. 

Order  of  A.  (input). 

Row  dimension  of  A  as  specified  in  the  calling  program.  lA  must  Ihj  greater 
than  or  equal  to  N.  (input). 

If  D1  is  non-negative  on  input,  then  1)1  and  D2  will  be  components  of  the 
detenninant  on  output  such  that  detemiinant  (A)  -  D1  *2**!  )2. 


WKARliA  Work  area  iength  at  least  N  when  IJOB  2. 


IfiR 


I'lTor  parameter.  I  erminai  error  -  r.^K  fN,  wlieieN 
A  is  algoritlimiraliy  singular. 


{ 


(t 


2  iiuiieates  dial  iiiairix 


C.L23  Subroutine  LUDATF.  Called  by  L1NV3F;  pcrfonns  L-U  decomposition  by  the  Croul 
algorithm  with  optional  acctjracy  test. 

A  Input  matrix  of  dimension  N  x  N  containing  the  matrix  to  be  decomposed. 

LU  Real  output  matrix  of  dimension  N  x  N  containing  the  L-U  decomposition 

of  a  rowwise  permutation  of  the  input  matrix. 

N  Input  scalar  containing  the  order  of  the  matrix  A. 

lA  Input  scalar  containing  the  row  dimension  of  matrices  A  and  LU  in  the 

calling  program. 

IDGT  Input  option. 

If  IDGT  is  greater  than  zero,  the  non-zero  elements  of  A  are 
assumed  to  be  correct  to  IDGT  decimal  places.  LUDATF  performs 
an  accaracy  test  to  determine  if  the  computed  decomposition  is  the 
exact  decomposition  of  a  matrix  which  differs  from  the  given  one  by 
less  than  its  uncertainty. 

If  IDGT  is  equal  to  zero,  the  accuracy  test  is  bypassed. 

D1  Output  scalar  containing  one  of  the  two  components  of  the  determinant.  See 

description  of  paraircter  D2,  below. 

D2  Output  scalar  containing  one  of  the  two  components  of  the  determinant. 

The  determinant  may  be  evaluated  as  (Dl)  (2**D2). 

IPVT  Output  vector  of  length  N  containing  the  permutation  indices. 

EQUIL  Output  vector  of  length  N  containing  rexiprocals  of  the  absolute  values  of 
the  largest  (in  absolute  value)  element  in  each  row. 

WA  Accuracy  test  paran*eter,  output  only  if  IDGT  is  greater  than  zero. 

EER  Error  paranK;ter, 

Terminal  error  =  1 28-rN. 

N  =  1  indicates  that  matrix  A  is  algorithmically  singular. 
Warning  error  =  32-fN. 

N  =-  2  indicates  that  the  accuracy  test  failed.  The  computed 
solution  may  be  in  error  by  more  than  can  be  accounted  for 
by  the  uncertainty  of  the  data.  This  warning  can  be 
prcxJuced  only  if  IE)GT  is  greater  than  0  on  input. 

C.l  .24  Subroutine  LUELMF.  Called  by  LINV31';  jrerforms  the  elimination  part  of  the  solution  of 
AX  -  B,  in  full-storage  rnixle. 

A  'Fhe  result,  I.U,  computed  in  the  subroutine  LL'DA  ri',  where  I.  is  a  lower 

tiiangular  matrix  with  ones  on  the  main  diagonal.  IJ  i.s  upper  u  iangular.  ! . 
and  1)  are  sto.ad  as  a  single  matrix  A,  and  the  unit  diagonal  of  I  is  not 
stored. 

B  is  a  vet  tor  of  iergth  N  on  the  right  hant!  side  of  the  eipianoii  AX  =  B 


( '  7 


B 


IPVl’  'Fhc  fjcm-sutaliosi  rnamK.  nitunied  ftoin  the  subroutine  LT  IDAl’F,  storcd  as  an 

N-iength  vector. 

N  Order  of  A  and  number  of  rows  in  B. 

TA  Nuniter  of  txtws  in  the  ilimt'-nsion  statement  for  A  in  the  calling  program, 

X  IT'^e  result  X. 

C.1.25  Subroutine  QUTPT.  Called  by  MAIN;  writes  output  data,  along  with  headings,  on  the 
output  file.  The  paratneter  lOUT(J)  lirotn  input  deteiTnines  whether  the  output  file  receives  data  of 
Type  J.  For  example,  ouiput  category  No.  1  is  occupant  segment  position  information.  If 
lOUT(l)  =  1,  these  data  go  to  output;  if  lOUT(l)  =  0,  they  do  not.  Also  performs  filtering  of 
acceleration  data  if  requested  in  input.  Called  by  EQUATE,  F:QUAT2,  FOBODY,  or  SOLVE  in 
the  event  of  abnormal  termination. 

C.1.26  Subroutine  PLOT!'.  Called  by  OUTFf;  pit)vides  printer  plots  for  up  to  three  dependent, 
continuous,  single-valued  functions  (Yl,  Y2,  Y3)  against  an  even-incremental  independent 
variable  (XT 

M  The  number  of  dependent  variables  ( 1 , 2,  or  3). 

N P  The  number  of  points  to  be  plotted  for  each  dependent  variable. 

X  The  independent  variable. 

Yl  llie  dependent  variables. 

Y2 

Y3 

C.i.27  Subroutine  POSTON  Called  by  QMATRX  or  QMATX2;  uses  equations  of  the  fc»rm 


to  compute  absolute  positions  of  29  points  on  bt>dy  (XC,  YC,  ZC).  Computes  positions  of  the 
same  29  points  in  aircraft  coordinate  system  (XCA,  YCA,  ZCA).  Calculates  velocities  IXCUA, 
YCDA,  ZCDA)  for  output  and  for  use  in  vekx:ity  deixjndent  form  cominitation. 

C.  1.28  Subroutine  QMATRX.  Called  by  EQUATE  for  the  three  dimensional  model;  calculates 
elements  of  generalizea  force  vector  [Qf ).  Calls  FORCES  for  computtition  of  external  forces 
acting  on  occupant. 

C.1.29  Subroutine  OMATX2  Called  by  LQIJAT2  for  the  plane-niotion  model;  calculates 
elements  of  generalized  force  vector  (Qf).  Calls  FORCES  for  computation  of  ex  ternal  forces 
acting  on  occupant. 

(’.’,30  Subroutine  K liST'R'T.  ('ailed  by  MAIN  at  unput  S[H't:ified  imeiN'a!  .  fo  tn-Mi-rale  da.ta  tiles 
on  unit  25  for  resiariing  solution  at  time  later  tluin  /ero.  Called  by  siduiiou  i;-, 

resiarted  'o  reat-  data  previously  stored  <>(>  uiiit  3.^ 


( ■  K 


C.1.'3l  Subroutiric  RKAM.  Called  by  MAIN;  salves  a  set  of  N  simultaneous,  rirst-orde!', 
OKiinai7  differential  equations.  Because  of  the  iiniwrlance  of  t)ie  integration  scheme  to  the  success 
of  any  dynamic  analysis  prograni,  a  detailed  discussion  of  the  method  is  provided  along  with  the 
description  of  the  FOR'FR  AN  subiDutine. 

Method  -  Tlie  user  is  allowed  an  option  of  using  either  the  Runge-Kutta  classical  fourth -order 
method  or  the  Adaivis-Moulton  pmdictor-corrcctor  method  using  the  Range  Kipta  meti/od  for 
starting  the  prrtcess. 

'Fite  system  of  equation, s  to  be  solved  is: 


yi’  =  fi(x.  yi.  y2-.  yp^) 

yi(><o)  =  yio 


i  =  1, N  (c.i) 


Let  yin  be  the  value  of  y;  at  x  =  Xn  and  fin  the  derivative  of  yi  at  x  ~  Xj,.  and  let  h  be  the  increment 
(step  size)  of  the  independent  variable  x.  The  classical  Runge-Kutta  fourth-order  method  uses  the 
fonnulas 


kii  =  hf;(x„,  yin), 

ki2  -  hfi(xn  +  V2  h,  yin  +  V2  kil)> 

ki3  =  hfi(xn  +  V2  h,  yin  +  V2  ki2),  (C.2) 

ki4  hfi(xn  +  H.  yin  f  ki3), 

yi.n-H  =  yn  +  V6  (kil  +  2ki2  +  2ki3  +  ki4) 


The  nonnai  option  is  to  continue  the  integration  with  Adams- Moultou  predictor-corrector  forrmilas 
once  enough  back  values  have  been  gcnenitcd  by  the  Range  Kutta  metlKKi. 

The  Adams-Moulton  predictor-comxtor  formulas  for  the  system  ((  ’  1)  arc 


^  ^  bdf,.,  ,  a  -  9f\,  3) 


V..,  t  {yf’’'  ,  f  ICf  ,,  ■  bt,.,  i  a  f,  „  \ 

J  j  {,  ^  j  .7  ^  j  5,  .  j  I,  J..i  J  l.n  .  J 


The  cora'Ctor  formula  (C.4)  is  applied  ooiy  once  per  s?c 
are  na^ded  for  each  Adams-MouJf.om  iotegratioj)  step, 
initially  obtaiiied  using  the  Rnnge-Kutta  iinethod,, 


p  so  tliat;  <fnh'  ovo  derivative  evahiadons 
The  slatting  valiu-ts  needed  irt  (C.3)  are 


The  AdamS'Moulton  methfti  m&y  i)e  used  with  eidicr  a  fixed  step  size  or  a  variable  step  size.  Tlie 
step  size  to  te  used  in  the  vaj'iable  mode  is  deienTjineci  fn>ni  the  diffeicttce  berovecn  tiie  predicted 
and  corrected  values.  The  integi'ation  step  .size  is  thus  cor.troUctt  dynasTticaliy  between  piescribed 
error  bounds  so  that  execudon  speed  and  accuracy  can  b&  optimized. 


Restrictioiis  -  An  auxiliary  loutine  must  be  provided  for  evaluation  of  the  fxist-ovder  derivatives. 
(See  AUXSUB  under  C:a]iing  Sequence.) 


Initial  conditions  for  both  variables  and  derivatives  must  be  stored  in  their  respective  locations  prior 
to  entering  RKAM. 


CalLitjg  Sequence 
XIJP 
HDP 
V.AR 
DER 

AUXSUB 

N 

OPT 

EU 

e. 

UMAX 

HMIN 

iMX'T 

TFMl-'.S 

\'i-i 


X,  tlie  independent  variable 
h,  the  integration  step  size 

N-diinensional  vector  of  dependent  variables  (yi,  yi,—^  Vn) 

N-dimensional  vector  of  derivatives  (yj',  y2’i."-.yn') 

Name  of  the  atixiliarj'  routine  that  computes  derivatives  and  stores 
them  in  DFiR(l)  to  DER(N).  The  main  program,  which  calls 
RKAM,  must  contain  an  EXTERNAL  statement.  No  items  are 
allowed  in  the  calling  sequence. 

Numlxer  of  etiuaiions 

Option  indicator,  zero  for  AM,  non  .lero  for  RK  only 
N-dimensional  vector  of  upjier  bound.i  from  m<dn  pi-ogram 
N-dimensiona!  vector  of  lov»^er  Ixtunds  froio  main  program 
Absolute  value  of  ntaximuni  aiiowable  step  si;® 

Absolute  value  of  minimum  allowable  step  size  (IIMIN  ;>  0) 

Interna!  counter,  .set  K>  zero  iniiiaby  in  MAIN 


A  two  dimens'otiai,  (h.N)  stoiaj-.e  region,  ITMP  f 

inns!  be  set  ti)  zero  initially  or  when  irstaiting 


1.  N 
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FunctiOiTial  Ikiscdplica  -  The  subroutine  employs  the  fourth-order  Adams-Moulton  predictor- 
cfH'recfor  trrevfsod  usinji  the  classioal  foutiJi-ortler  Runge-Kutta  method  to  obtain  starting  values. 

AM  has  the  fodowing  advjintage.s  witli  respect  to  RK: 

1 .  Gv'tly  half  a.s  many  derivative  ev.aiaations  per  integration  step  are  required  to 
attain  the  same  order  of  accuracy. 

2.  The  local  truncation  error  may  be  estimated  at  the  conclusion  of  each 
integration  step  thereby  providing  a  means  for  stip  size  control. 

For  each  vaiiable,  the  h>cal  truncation  error  is  approximately  one-fourteenth  the  difference  between 
tire  predicted  aiid  coirected  values,  that  is 


In  RKAM,  die  differences  Dj  =  y^T^  are  formed  and  cornpaied  with  positive  numbers  EUi 
and  ELj.  If  Dj  ^  EUi  for  ariv  i,  the  sic  p  size  is  halved  provided  lh/2!'  >  HMIN.  If  Dj  <  ELj  for  gll  i 
and  lor  three  successive  steps,  the  step  size  is  doubled  provided  i2hl  <  HMAX.  (Note  that  h  may 
be  held  fi.xcd  either  by  setting  HMIN  -  HMAX  or  by  making  EUj  and  ELi  prohibitively  large  and 
small,  respectively.)  If  halving 's  ccJled  for  during  the  first  AM  step  following  the  three  initial  RK 
steps,  the  .step  size  is  halved,  Ae  independent  \'ariable  is  set  back  to  its  initial  value,  and  the  three 
RK  steps  are  repeated.  This  v/U!  continue  until  the  first  AM  step  »s  successfully  taken.  From  this 
point  on,  halving  is  effected  by  interpolation  of  past  data  whereas  doubling  is  accomplished  by 
alternate  selection  of  past  data. 

Ill  selecting  F,U  and  EL,  one  should  note  the  following: 

1 .  Fhe  test  is  an  absolute  test.  To  control  relative  error  EUj  and  ELj  should  be 
computed  as  functions  of  y^  prior  to  each  integration  step. 

2.  Although  tlie  Ux:al  truncation  eiror  in  is  not  allowed  to  exceed  EUi, 
dex^s  not  imply  that  die  '.'umulative  error  will  not  e,cceed  FlJi,  Therefore, 

FiUj  and  EL.i  should  depend  upon  the  maximum  aiiowable  cumulative  error 
and  the  numbcj  of  integ  iition  steps. 

.1,  Since  doubling  h  will  multiply  the  truncation  error  by  a  factor  ol  2.*’,  El.; 

sl’cruhi  tie  chosec.  less  than  EU,/32  if  the  advantages  of  Joubling  a’c  not  to 
bit  short  lived. 


C.  1..12  Subroutine  KMA I'RX  vrilicd  by  EQUATE  tor  the  tham-dimfisional  nKtciel,  calculates 
elen'icnts  of  joint  resistance  vector  IK),  ln(nit  paraiiKder  IMAN  (klenniucs  whether  hiiniaii  (IMAN 
Oi  or  dummy  (  IM.3N  -  1)  rnociel  is  iiseii. 

Shb.'OuUric  RMATX2.  ('ailed  by  EtyLlA'i'?  Uu' the  p!am;  motion  n.iKid  caicuiaic'^ 
clciocats  of  join:  noastatu.c  o'Ctor  (R  ' 


U. !  S.ul'ifUiUUt.b.F.AXii'i'  ('aiicd  b>  iNl’T;  reads  innut  d  iM  required  lor  rigid  sc  u  model  atid 
loicrgv  atworbing  op'^jir 


C.  1,35  Subroutine  VMATRX.  Called  by  EiQUATB  for  the  three  dimensional  model;  calculates 
elements  of  force  vech^r  (Fp}  (jerived  from  system  potential  energy. 


C.  1.36  Subroutine  VMATX2,  Called  by  EQUAT?  for  the.  plane-motion  ntodcl;  calculates 
elements  of  force  vector  {Fp '  derived  from  system  potential  energy. 


c.: 


.miiom 


'Fhe  relationships  among  die  subroutines  in  the  seat  segment  of  the  program  arc  illustrated  iri  Idgure 
C-3.  Individual  subroutines  are  described  below. 


C.2. 1  Subroutine  ASSBLE.  Called  by  MAIN;  initializes  the  elemeni  data  storage.  The  mass 
matrix  and  the  initial  transformations  Bji  for  the  tiodaJ  coordinate  systems  are  assembled,  and  the 

initial  values  of  the  pointing  vectors  ri,  "n,  and  ^  and  the  normal  components  of  the  rigid  links  A  are 
generated. 

C.2.2  Subroutine  ASSMBL.  Called  by  PLSTF  and  BMSTF;  assembles  the  master  stiffness 
matrix  in  a  banded  symmetric  form.  Tltis  subroutine  calls  subroutine  KADE,  which  adds  a 
particular  element  of  the  square  element  stiffness  matrix  to  the  banded  master  stiffness  matrix. 


C.2.3  Subroutine  BASME.  Called  by  ASSBLE;  forms  the  initial  element  coordinate  system  E  for 
beam  and  spring  elements. 


C.2.4  Subroutine  BFRCIN.  Called  by  FRCIN;  calculates  tlie  beam  and  spring  element 
deformations  and  nodal  forces  in  the  element  coordinate  system.  Performs  the  operations 
asscjciated  with  the  master-slave  lelaiions  and  transfomss  the  forces  to  the  nodal  ctxtrdinate  system. 

C.2.5  Subroutine  BGEQM.  Called  by  READ.IN,.  reads  the  data  describing  the  cross-section 
projteilies  of  beams  mid  springs.  Generates  certain  additional  data,  such  as  segment  lengths  and 
torsional  constants. 


C.2.6  Subroutine  BMEND.  Called  by  BMSTF;  calculates  the  reduced  stiffness  matrix  due  to 
axial  force,  shear,  and  morrtent  discontinuities 


C.2.7  Subroutine  BMSTF.  Called  by  SOLVE;  calculates  beam  or  spring  element  stiffness  matrix. 
The  principal  subroutines  called  include  BMSTFl  for  clastic  material,  BMSTF2  for  inelastic 
material,  BMEND  for  the  irKdif  cation  of  the  stiffness  due  to  special  end  conditions,  and  ASSMBL 
for  assembly  of  element  stiffness. 


C.2.8  Subroutine  BMS'lTl.  Called  by  BMS'i'l  ';  calculates  the  cacuic  stiffness 
or  spring  eien^ent. 


inatri.x  for  a  beam 


C.2.9  Subroutine  BMS  rF2-  Called  by  BMSTF,  calcnbifc's  liic  uiiigentKi!  stiffness  for  a 

iHttun  or  spnng  element. 

('’.2. 10  Suhroutine  BCi)lJND  Called  by  .SOLVii;  applies  the  spccihcd  bo!.mda''Y  coiuiitions  fu  the 
assembh  (i  master  stiffness  matrix. 


C',2. !  i  Snhroiiiiiic  (’R(..)SS.  Utility  subniiitine;  ciilculatcs  the  (irndiit  t  rf  two  ui.Uncc:,. 
C.2,i2  .S,u!}roUiUlLiTILV-.T  (’ailed  l\v  f.i'T.S'l'b:  coiiunn.s  nliitc  iKnulim:  cii;  \  atait:  tithcas 


ft/ 


iTVi'iiittimmimrwr  inriiT 


C.2.13  Subroutine  CUR  V  AT.  Called  by  'fT-'RCIN;  provides  algebraic  expressions  for  curvature 
components  at  the  midpoints  of  the  three  sides  of  the  plate  elements  rs  functions  of  !he  nodal 
rotations. 

C.2.14  Subroutine  DECOD.  Utility  subroutine,  deco^les  a  packed  word. 

C.2.15  Subroutine  DGMPRD.  Double  [srecision  version  of  G.M  PRD. 

C.2.16  Subroutine  DGTPRD.  Double  precision  version  of  G  ITRD, 

C.2.17  Subroutine  EPTSTF.  Called  by  PLSTF;  calculates  plate  element  stiffness  matrix.  This 
subroutine  calls  TRIANG  for  the  in-plane  stiain-displacement  relationship,  CRVTBL  for 
curvature,  and  FORMK  for  calculation  of  appropriate  elements  in  the  stiffness  matrix. 

C.2.18  Subroutine  F/"OG.  Called  by  BMSTF  and  PLSTF;  transfonns  appropriate  variables  from 
the  element  coordinate  system  to  the  global  coordinate  system, 

C.2.19  Subroutine  FORMK.  Called  by  EPTSTF;  calculates  the  products  oi  three  different 
matrices. 

C.2.20  Subroutine  FRCIN.  Called  by  SOLVE;  calruiates  internal  ncxial  forces.  The  program 
updates  tire  nodal  coordinate  transformations  B.  and  calls  subroutines  Tl  'RCIN  for  plate  forces  and 
BFRCIN  for  beam  forces. 

C.2.21  Siibroudne  FREEFD.  Called  by  SOLVE;  calculates  the  external  forces  including  restraint 
system  forces  and  forces  exerted  by  the  (X-’ciipants  on  the  see.t  [ran  ar,  J  seat  back. 

C.2.22  Subroutine  FWARP.  Called  by  SOLVE'  modifies  forces  to  acc«'iint  for  specified  HcKrr 
warp,  displacements  and  rotations. 

C.2.23  Subroutine  GENldfM.  Called  by  LOCFRC,  numericallv  integrates  .stresses  over  the  cross 
section  of  the  Iream  to  (rbtaio  internal  fotves  and  moments. 

(’.2.24  Subroutine  GENr01/GFNF()2.  ('alUd  by  'i}'K(’’N;  coinpiues  moments  and  .forces  at  a 
cross  section  of  an  clastic-plastic  plate  (wiih/widKnit')  integrating  tiirough  the  thickness, 

(’.2.25  Subfouiine  GMPRD.  Utility  sebrc-riin'';  iXTUirms  general  mairix  nniltiplii  ation 

(.'.2.26  Subroutine  GTPRLV  Udiity  suhu)iitine:  calcul  no,  tire  [iroinct  of  ilie  tiar)S[rosc  of  a  matns 
wiili  another  matrix. 

(’  2.2?  Suhroiiliiie  KADD  ('ailed  by  A.S.SMHI  .,  add-  a  parliciihu  elenient  ot  ilie  vvjuare  mauix  to 
the  banded  mairix. 

Subroutine  i  .tX’hRC’.  ('ailed  by  eaf'uiaics  mnlplane  stra.ins,  yi,rv:itv,ires,  ruxlal 

lorccs,  ,ind  iiioments  in  ihe  he.un  cleinent  cisiniinatc  sy.su  ui,  1  Join  ai.on  ainl  iindai  fotoc'.  .at-  also 
.iH  i,ii;ncd  tor  the.  spring  in  the  r'eriKUii  c(M>rdin;i!c  svsu-sn. 

(',2  ,.’d  Subroutine  1,a..>CS B.  t  'alleti  t’y  HMS'!'!'.  lR)UNi),  and  K.AUft,  eompuu's  lue  I.K'aoon  sg  a 
p.iiiiculat  element  ot  a  stinau'  matrix  ''Vtu-ii  asseinbletl  inlo  th.e  bandcvi  ss  tn.mcnu'  tiMin, 


C.2.30  Sabrootine  MCHB.  Called  by  SOLVE;  solves  the  linear  system  of  equations  H  x  -  F  for 
K  (displacements).  The  master  stiffness  matrix  K  .is  assumed  to  be  symmetric  positive  definite  and 
stored  in  the  compressed  form,  that  is.  main  diagonal  and  upper  codiagonals  rowwise  in 
successive  storage  locaticns.  F  ts  the  applied  force  vector.  This  is  a  two-step  equation  solver  that 
uses  Choksky's  method.  In  the  first  step  the  master  stiffness  matrix  K  is  factored  into  an  upper 
diagonal  rikatrix  U  and  a  lower  diagonal  matrix  L. 


.K  =  Ln 


(C.6) 


and  let 


Ux  =  V 


(C.7) 


so  that  the  linear  system  of  equations  Kx  =  F  is  equivalent  to 
Lv  =  F 


(C.8) 


In  the  second  step,  equation  (C.8)  is  solved  by  forward  reduction  for  v  and  finally  equation  (C.7) 
is  solverl  for  x  by  back  substituting  foi'  v. 

C.2,31  Subroutines  MISES1/MISES2.  Called  by  'FFRCIN;  computes  biaxial  elastic-plastic 
stress-strain  relations  using  Von  Mises  yield  criterion. 

C.2.31  Subroutine  NODALF.  Called  by  ITRCIN;  calculates  nodal  forces  and  moments  for  the 
plate  element. 

C.2.33  Subroutine  NTOCi.  Called  PLSTF;  transforms  appropriate  variable, s  from  the  global 
coordinate  system  to  the  nodal  coordinate  system. 

C.2.34  Subroutine  OUTPTS.  Called  by  SOL’^^'E;  organizes  and  tabulates  output  for  those 
quantities  selected  for  output.  Deformed  seat  mrdel  plot  data  are  written  onto  file  20  at  user- 
selected  tiines. 

C.2.35  Subroutine  PLSTF.  Called  by  SOLVE;  calls  liP'rS  l’F  to  fonn  the  plate  clement  stiffness 
matrix  and  then  uses  ASSMBI  to  assemble  the  element  stiffness. 

r.3.36  Subroutine  READIN.  Called  by  MAIN;  a^ads  all  input  data  and,  if  required,  initializes  the 
data  files.  Undefomied  seat  rntxle!  data  and  mtKlel  parameters  aa*  written  onto  file  20  if  recjiic'.ted 

('.2. .17  Subroutine  .SOIA'Ii,  ('ailed  by  MAIN;  performs  the  main  solution  [irocedure.  The 
principal  subroutines  called  are  HMSTl  for  beam  or  spring  stiffness,  P!  STl  tor  plate  sliftiiess, 
I'RFFFD  for  ajiplied  forces,  MC'llB  for  the  solutioi  of  displacements,  ami  OUTi^TS  for  printed 
output  and  plot  of  selectetl  paraiiieters. 

('.2,3b  Subrouiinc  SPRi.Nt  i  ('ailed  from  I .(,)( 'K [■( ',  ealeniates  ihe  element  forces  for  i  spiing 
elemcin. 


('.2.3d  .Suhioutine  .SJ'RI'.S  ('ailed  tiom  ('diN'f  ltM  and  SPKiNdi;  pno  ules  .ui  a^  'luiliini  foi 
uniaKiai  aiess  strain  relationshij). 

(,',e  •}('  .Siihroniine  L'VSMl.  (  alha!  trorn  .A.SSH!  h,  :,  .d(  fa  enfUnt^  ooii'  lo  fiic  lumptsi 

mass  matrix  Im  iilate  eleni!  in.  i  oinrs  tlie  inmal  eleniiiu  ef.i|i!jn  ne  s's.lsin  i  toi  (he  piaie 
elemen^.^, 


C.2.4T  SubrolUiiiC  TFRCIN.  Caikii  iron'.  FKCIM;  calculiarcs  plave  eJen'eni  def<rnv4tion'j  ra\{ 
forces  in  the  elernern  cwirdinate  system  'lliC  fovi  .-s  are  then  iiansfoirned  h?  ’.he  'lOtiai  ctK'uiinat;; 

system. 

C.2.42  Subroutine  TFvMNG  Ca’led  frein  EFFSFT;  coniains  aigebraic  ^.xaress’ons  of  tisr-  strain- 
dispiaccment  relationship  for  a  plate  element. 

C.2.4.?  Subroutine  V^CITRl.  Called  from  TP'RCiN;  calculates  the  defonned  lengths  of  plate 
eiemeni  sides. 


C.2.44  Sttltrouttne  VfiCrR2.  Called  by  TFRCIN;  determines  th.e  component  s  of  .a  vector  rionijai 
to  a  reference  .siviface. 


APPENDIX  D 
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